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Abstract 

Traditional  signal  processing  algorithms  rely  heavily  on  models  that  are  inherently 
linear.  Such  models  are  attractive  both  for  their  mathematical  tractability  and  their 
applicability  to  the  rich  class  of  signals  that  can  be  represented  with  Fourier  methods. 
Nonlinear  systems  that  support  soliton  solutions  share  many  of  the  properties  that 
make  linear  systems  attractive  from  an  engineering  standpoint.  Although  nonlinear, 
these  systems  are  solvable  through  inverse  scattering,  a  technique  analogous  to  the 
Fourier  transform  for  linear  systems.  Solitons  are  eigenfunctions  of  these  systems 
which  satisfy  a  nonlinear  form  of  superposition  and  display  rich  signal  dynamics  as 
they  interact.  By  using  solitons  for  signal  synthesis,  the  corresponding  nonlinear 
systems  become  specialized  signal  processors  which  are  naturally  suited  to  a  number 
of  complex  signal  processing  tasks.  Specific  analog  circuits  can  generate  soliton  signals 
and  can  be  used  as  natural  multiplexers  and  demultiplexers  in  a  number  of  potential 
soliton-based  wireless  communication  applications.  These  circuits  play  an  important 
role  in  investigating  the  effects  of  noise  on  soliton  behavior.  Finally,  the  soliton  signal 
dynamics  also  provide  a  mechanism  for  decreasing  transmitted  signal  energy  while 
enhancing  signal  detection  and  parameter  estimation  performance. 
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Chapter  1 


Introduction 


Many  traditional  signal  processing  applications  rely  on  models  that  are  inherently  lin¬ 
ear  and  time-invariant  (LTI).  Much  of  the  success  of  such  methods  can  be  attributed 
to  their  being  mathematically  tractable,  often  leading  to  efficient  signal  representa¬ 
tions  and  fast  algorithms.  Linear  techniques  have  also  proven  effective  for  modeling 
a  variety  of  signals  of  practical  interest  such  as  speech  or  financial  time-series  and 
systems  of  interest  such  as  the  telephone  or  radio  broadcast  channels.  However,  we 
increasingly  turn  to  nonlinear  models  to  capture  some  of  the  more  salient  behavior  of 
physical  or  natural  systems  that  cannot  be  expressed  by  linear  means,  such  as  thresh¬ 
old  phenomena,  amplitude-dependence,  or  chaotic  behavior.  Nonlinear  systems  also 
hold  the  potential  to  produce  more  efficient  algorithms  or  models  for  a  variety  of  signal 
processing  and  communication  problems  where  linear  techniques  are  suboptimal. 

Systems  that  support  solitons  may  be  a  natural  choice  for  a  class  of  nonlinear 
systems  to  explore  since  they  share  many  of  the  properties  that  make  LTI  systems 
attractive  from  an  engineering  standpoint.  Although  nonlinear,  these  systems  are 
solvable  through  inverse  scattering,  a  technique  analogous  to  the  Fourier  transform 
for  linear  systems  [1] .  Solitons  are  eigenfunctions  of  these  systems  which  satisfy  a  non¬ 
linear  form  of  superposition.  We  can  therefore  decompose  complex  solutions  in  terms 
of  a  class  of  signals  with  simple  dynamical  structure.  Solitons  have  been  observed  in  a 
variety  of  natural  phenomena  from  water  and  plasma  waves  [31,  57]  to  crystal  lattice 
vibrations  [14]  and  energy  transport  in  proteins  [31].  Solitons  can  also  be  found  in 
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a  number  of  man-made  media  including  super-conducting  transmission  lines  [56]  and 
nonlinear  circuits  [29,  60].  Recently,  solitons  have  become  of  significant  interest  for 
optical  telecommunications,  where  optical  pulses  have  been  shown  to  propagate  as 
solitons  for  tremendous  distances  without  significant  loss  or  dispersion  [23]. 

In  this  thesis,  we  view  solitons  from  a  decidedly  different  perspective.  Rather  than 
focusing  on  the  propagation  of  solitons  over  nonlinear  channels,  we  consider  using 
these  nonlinear  systems  to  both  generate  and  process  signals  for  transmission  over 
traditional  linear  channels.  By  using  solitons  for  signal  synthesis,  the  corresponding 
nonlinear  systems  become  specialized  signal  processors  which  are  naturally  suited 
to  a  number  of  complex  signal  processing  tasks.  This  thesis  can  be  viewed  as  an 
exploration  of  the  properties  of  solitons  as  signals.  In  the  process,  we  explore  the 
possibility  of  using  these  signals  in  a  potential  multi-user  wireless  communication 
context. 

For  example,  we  consider  the  problem  of  multiplexing  a  number  of  users  onto  a 
single  carrier  for  transmission  over  a  lossy  channel.  A  variety  of  linear  multiplexing 
schemes  have  been  proposed  such  as  time-,  frequency-,  or  code-division  multiple  access 
systems.  However,  the  number  of  nonlinear  multiplexing  strategies  available  is  rather 
limited.  One  potential  benefit  to  using  solitons  as  carrier  signals  and  the  nonlinear 
systems  as  multiplexors,  is  that  the  soliton  signal  dynamics  provide  a  mechanism  for 
simultaneously  decreasing  transmitted  signal  energy  and  enhancing  communication 
performance. 

1.1  Outline  of  the  Thesis 

There  is  a  large  body  of  literature  on  soliton  theory  spanning  over  a  century  of  re¬ 
search.  Rather  than  attempting  to  provide  a  self-contained  summary,  in  Chapter  2 
we  present  a  brief  overview  of  the  components  of  soliton  theory  that  we  will  draw 
upon  throughout  the  thesis.  The  chapter  includes  an  introduction  to  solitons  and 
some  of  the  nonlinear  systems  that  support  them.  We  focus  our  discussion  on  the 
Toda  lattice,  a  particularly  simple  soliton  system  that  forms  the  basis  of  many  of  our 
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examples  throughout  the  thesis.  In  order  to  exploit  the  rich  mathematical  structure 
of  these  systems,  we  review  some  elements  from  inverse  scattering  theory.  In  addi¬ 
tion  to  providing  an  efficient  mechanism  for  the  construction  of  soliton  signals  and 
the  solutions  to  their  dynamics,  the  inverse  scattering  framework  can  also  be  used 
to  synthesize  families  of  soliton  systems.  Throughout  the  thesis,  the  inverse  scatter¬ 
ing  transform  plays  an  analogous  role  to  the  Fourier  transform  in  the  analysis  and 
processing  of  soliton  signals  in  the  presence  of  noise. 

To  facilitate  real-time  generation  and  processing  of  soliton  signals,  it  will  be  im¬ 
portant  to  explore  implementations  of  these  nonlinear  systems.  In  Chapter  3  we 
develop  new  circuit  models  for  two  soliton  systems.  The  first  is  a  diode  ladder  imple¬ 
mentation  of  the  Toda  lattice  which  appears  to  be  the  first  circuit  implementation  to 
display  true  soliton  behavior.  We  also  develop  a  lattice-circuit  implementation  of  the 
discrete-KdV  equation,  which  will  be  important  for  processing  discrete-time  soliton 
signals. 

These  circuit  models  then  form  the  basis  for  a  communication  paradigm  presented 
in  Chapter  4,  where  multiple  signals  can  be  multiplexed  onto  soliton  carriers  using 
such  circuits  as  tuned  modulators  and  demodulators.  As  we  will  see,  the  nonlinear  in¬ 
teraction  of  multiple  solitons  can  be  exploited  as  a  means  for  reducing  the  transmitted 
signal  energy  in  a  multi-user  communication  context.  Such  low  power  transmission 
techniques  may  be  applicable  to  a  variety  of  portable  or  power-limited  communication 
applications. 

Before  soliton  systems  can  be  used  in  a  practical  communication  context,  we  need 
accurate  models  for  the  effects  of  disturbances  on  the  dynamics  of  these  systems.  The 
robustness  of  such  systems  to  additive  corruption  will  have  a  direct  impact  on  the 
demodulation  performance  of  the  nonlinear  receivers.  In  Chapter  5,  we  analyze  the 
effects  of  small  amplitude  corruption  on  the  dynamics  of  solitons  in  the  Toda  lattice 
and  characterize  the  statistics  of  the  noise  as  it  is  propagates  through  the  system. 

In  order  to  develop  effective  communication  strategies  by  modulating  the  param¬ 
eters  of  soliton  carriers,  we  need  to  have  accurate  models  for  the  ability  of  a  receiver 
to  resolve  these  parameters.  In  Chapter  6,  we  compute  Cramer-Rao  bounds  for  their 
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estimation  error  variance.  We  show  that  in  addition  to  reducing  signal  energy,  the 
nonlinear  interaction  of  multiple  solitons  can  also  enhance  parameter  estimation  per¬ 
formance.  Based  on  our  characterization  of  the  noise,  we  develop  a  set  of  parameter 
estimation  algorithms  in  which  maximum-likelihood  (ML)  estimates  can  be  obtained 
from  corrupted  measurements.  In  Chapter  7,  we  demonstrate  how  soliton  circuits  can 
be  used  to  enhance  the  detection  of  multiple  solitons  in  noise. 

Finally,  Chapter  8  summarizes  the  main  contributions  of  the  thesis  and  indicates 
some  interesting  and  potentially  important  directions  for  future  study. 
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Chapter  2 


Soliton  Systems 


Solitons  play  an  important  role  in  the  study  of  a  large  class  of  nonlinear  evolution 
equations.  As  we  shall  see,  the  ability  to  describe  their  long-term  behavior  analytically 
makes  this  class  of  systems  attractive  for  modeling  a  variety  of  nonlinear  phenomena 
in  a  number  of  diverse  areas  in  mathematics,  physics,  biology,  and  engineering.  Ex¬ 
amples  include  topics  from  surface  or  internal  water  waves,  to  energy  transport  along 
long  protein  chains,  and  the  propagation  of  optical  pulses  along  nonlinear  fibers.  Fur¬ 
thermore,  the  development  of  optical  and  electrical  analogs  for  many  of  these  systems 
makes  soliton  signals  of  great  practical  interest. 

The  theory  of  solitons  dates  back  to  the  work  of  Korteweg  and  de- Vries,  and 
was  motivated  by  attempts  to  explain  the  unusual  water  wave  observations  of  Scott 
Russell  in  1834.  The  method  of  solution  for  soliton  systems  began  with  the  work  of 
Gardner,  Greene,  Kruskal,  and  Miura  on  the  Korteweg-de  Vries  (KdV)  equation  [18]. 
These  techniques  were  generalized  by  Lax  [42]  and  then  by  Ablowitz  et  al.  [1]  and 
used  to  solve  what  is  now  a  large  class  of  solvable  nonlinear  evolution  equations. 

This  chapter  is  designed  as  a  brief  overview  of  the  components  of  soliton  theory 
that  we  will  draw  upon  throughout  the  thesis.  Although  no  new  results  will  be 
presented,  this  chapter  serves  several  purposes.  First,  we  present  an  introduction  to 
solitons  and  to  some  of  the  nonlinear  systems  that  support  them.  We  will  focus  our 
attention  on  the  Toda  lattice,  a  simple  mechanical  system  that  will  form  the  basis 
of  many  of  the  examples  throughout  the  thesis.  Second,  we  summarize  some  of  the 
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Figure  2-1:  The  Toda  Lattice. 

main  results  of  inverse  scattering  theory,  the  method  by  which  soliton  systems  can  be 
solved  analytically.  Rather  than  attempting  to  fully  develop  this  topic,  we  concentrate 
on  portions  of  the  theory  that  we  will  exploit  in  the  applications  developed  in  later 
chapters  of  the  thesis.  Finally,  we  mention  a  few  of  the  techniques  that  can  be  used  to 
construct  families  of  soliton  systems.  As  we  will  see  in  later  chapters,  these  systems 
are  potentially  useful  for  a  number  of  communication  applications.  Therefore,  there 
may  be  a  practical  use  for  the  generation  of  increasingly  complex  soliton  systems. 

2.1  The  Toda  Lattice 

The  Toda  lattice  is  a  conceptually  simple  mechanical  example  of  a  nonlinear  system 
with  soliton  solutions.  A  comprehensive  treatment  of  the  lattice  and  its  associated 
soliton  theory  can  be  found  in  the  monograph  by  Toda  [69].  As  shown  in  Fig.  2-1, 
the  Toda  lattice  equations  describe  the  displacements  of  an  infinite  chain  of  masses 
connected  with  nonlinear  springs.  Each  of  the  springs  satisfies  the  nonlinear  force  law 


fn  =  o(e 


-1), 


(2.1) 


where  fn  is  the  force  on  the  spring  between  masses  with  displacements  yn  and  yn~i 
from  their  rest  positions.  The  equations  of  motion  for  the  lattice  are  given  by 


myn  =  a  (e  b^Vn  yn~^  —  e  b(yn+1  yn^  , 


(2.2) 
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Figure  2-2:  A  propagating  wave  solution  to  the  Toda  lattice  equations.  Each  trace 
corresponds  to  the  force  fn(t)  stored  in  the  spring  between  mass  n  and  n  —  1. 

where  yn  is  the  displacement  of  the  n-th  mass  from  its  rest  position,  m  is  the  mass, 
and  a  and  b  are  constants.  This  equation  admits  pulse-like  solutions  of  the  form 

fn(t)  =  /?2sech2 (sinh~ 1  ( \Jm/ab  /3)n  -  fit),  (2.3) 

which  propagate  as  compressional  waves  stored  as  forces  in  the  nonlinear  springs.  A 
single  right-traveling  wave  fn(t )  is  shown  in  Fig.  2-2.  The  bottom  trace  in  the  figure 
corresponds  to  the  force  in  the  spring  between  masses  “zero”  and  “one”  of  an  infinite- 
length  Toda  chain.  This  compressional  wave  is  localized  in  time,  and  propagates  along 
the  chain  maintaining  constant  shape  and  velocity.  Since,  for  example,  the  wave 
fn(t)  appears  on  the  thirtieth  mass  at  a  later  point  in  time,  this  wave  is  therefore 
propagating  to  the  right  along  the  lattice  as  viewed  in  Fig.  2-1.  The  parameter  (3 
appears  in  both  the  amplitude  and  in  the  temporal-  and  spatial-scales  of  this  one- 
parameter  family  of  solutions  giving  rise  to  tall,  narrow  pulses  which  propagate  faster 
then  small,  wide  pulses.  This  type  of  localized  pulse-like  solution  is  what  is  often 
referred  to  as  a  solitary  wave. 


17 


Definition  1:  A  solitary  wave  solution  to  a  partial  differential  equation  for  de¬ 
pendent  variable,  y,  with  temporal  and  spatial  variables,  t  and  n,  is  a  traveling-wave 
solution  of  the  form, 

y(n,  t )  =  f(n  -  ct )  =  f(z),  (2.4) 


where  c  is  a  fixed  constant,  and  the  energy  of  f(z)  is  localized. 

The  history  of  solitary  waves  dates  back  to  the  work  of  John  Scott  Russell  in  1834 
and  perhaps  the  first  recorded  sighting  of  a  solitary  wave.  Here  is  his  original  account 
of  his  sighting: 

I  was  observing  the  motion  of  a  boat  which  was  rapidly  drawn  along 
a  narrow  channel  by  a  pair  of  horses,  when  the  boat  suddenly  stopped  — 
not  so  the  mass  of  water  in  the  channel  which  it  had  put  in  motion;  it  ac¬ 
cumulated  round  the  prow  of  the  vessel  in  a  state  of  violent  agitation,  then 
suddenly  leaving  it  behind,  rolled  forward  with  great  velocity,  assuming 
the  form  of  a  large  solitary  elevation,  a  rounded,  smooth  and  well  defined 
heap  of  water,  which  continued  its  course  along  the  channel  apparently 
without  change  of  form  or  diminution  of  speed.  I  followed  it  on  horseback 
and  overtook  it  still  rolling  on  at  a  rate  of  speed  some  eight  or  nine  miles 
an  hour,  preserving  its  original  figure  some  thirty  feet  long  and  a  foot  to 
a  foot  and  a  half  in  height.  Its  height  gradually  diminished  and  after  a 
chase  of  one  or  two  miles  I  lost  it  in  the  windings  of  the  channel.  Such  in 
the  month  of  August  1834  was  my  first  encounter  with  that  singular  and 
beautiful  phenomenon...  [58]. 

What  Scott  Russell  actually  observed  was  a  solitary  wave  solution  to  what  is 
now  known  as  the  KdV  equation  [57],  The  Toda  lattice  is  in  many  ways  analogous 
to  KdV;  in  fact,  KdV  can  be  derived  as  a  continuum  limit.  A  detailed  discussion  of 
linear  and  nonlinear  wave  theory  including  KdV  can  be  found  in  [75].  In  a  1965  paper, 
Zabusky  and  Kruskal  performed  numerical  experiments  with  KdV  and  noticed  that 
these  solitary  wave  solutions  retained  their  identity  upon  collision  with  other  solitary 
waves.  Since  the  velocities  of  KdV  solitary  waves  are  proportional  to  their  amplitudes, 
a  collision  of  solitary  waves  will  occur  for  any  solution  with  a  taller  pulse  to  the  left  of 
a  shorter  pulse.  As  the  individual  solitary  waves  approach  one-another,  they  begin  to 
interact  nonlinearly.  However,  after  passing  through  one-another,  they  regain  their 
shape  and  speed  with  only  a  slight  positional  shift  as  evidence  of  their  interaction  [78]. 
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This  process  indicates  that  the  nonlinear  dynamic  system  has  a  generalized  coordinate 
transformation,  or  underlying  dynamic  structure  which  produces  a  set  of  “normal 
modes”  that  are  almost  completely  uncoupled. 

Zabusky  and  Kruskal’s  discovery  that  these  solitary  wave  solutions  behaved  much 
like  particles  in  that  they  essentially  retained  their  identity  after  collisions  prompted 
them  to  coin  the  term  soliton  implying  that  these  solitary  waves  possess  a  particle¬ 
like  nature.  This  is  the  type  of  behavior  we  would  expect  for  linear  wave  equations, 
that  is,  the  ability  to  form  solutions  from  a  superposition  of  simpler  waves.  However, 
that  a  nonlinear  equation  such  as  KdV  or  the  Toda  lattice  permits  such  a  form  of 
superposition  is  an  indication  that  it  belongs  to  a  rather  remarkable  class  of  nonlinear 
systems.  This  property  distinguishes  solitons  from  other  solitary  waves.  We  will  adopt 
the  following  working  definition  of  solitons: 


Definition  2:  A  soliton  solution  to  a  partial  differential  equation  is  a  solitary  wave 
which  asymptotically  retains  its  shape  and  velocity  upon  collision  with  another  soliton, 


i.e. 


y{n ,  t) 


f{n  —  ct),  t  — >  — oo, 
f(n  —  ct  +  6),  t  — y  +oo, 


(2.5) 


for  (n  —  ct)  fixed  and  5  an  arbitrary  constant. 


As  we  shall  see,  the  term  soliton  often  implies  additional  mathematical  structure 
involving  the  inverse  scattering  transform. 

Returning  to  the  Toda  lattice,  Fig.  2-3  illustrates  soliton  behavior  for  two  solutions 
of  the  form  of  Eq.  (2.3).  The  bottom  trace  in  the  figure  corresponds  to  the  force  in  the 
spring  between  masses  with  indices  zero  and  one  in  the  lattice.  Note  that  as  a  function 
of  time,  a  smaller,  wider  soliton  appears  before  a  taller,  narrower  one.  However,  as 
viewed  by,  e.g.,  the  thirtieth  mass  in  the  lattice,  the  larger  soliton  appears  first  as  a 
function  of  time,  i.e.  has  traveled  faster.  Note  that  when  the  larger  soliton  catches  up 
to  the  smaller  soliton  as  viewed  on  the  fifteenth  node,  the  combined  amplitude  of  the 
two  solitons  is  actually  less  than  would  be  expected  for  a  linear  system,  which  would 
display  a  linear  superposition  of  the  two  amplitudes.  Also,  the  signal  shape  changes 
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time 

Figure  2-3:  Two  solitary  wave  solutions  to  the  Toda  lattice. 


significantly  during  this  nonlinear  interaction.  Both  of  these  characteristics  of  soliton 
interaction  will  have  useful  implications  in  the  context  developed  in  Chap.  4. 

An  analytic  expression  for  the  two-soliton  solution  for  Pi  >  0  >  0  is  given  by  [29] 

f  rt\  _  ^^sech^r/x)  +/g|sech2(772)  +  Asech2(77i)sech2(7/2) 
n  ab  (cosh(<£/2)  +  sinh(^/2)  tanh(?7i)  tanh(j?2))2  ’  '  7 


where 


A  =  sinh(<£/2)  [[01  +  00s)  sinh(^/2)  +  20  p2  cosh (<£/2))  , 


<f>  =  In 


sinh((px  -p2)/2)\ 
sinh((pi  +p2)/2)J  ’ 


and  0  =  Jab/m  sinh(pj),  and  rji  =  Pin  —  0 (t  —  0). 

Although  Eq.  (2.6)  appears  rather  complex,  Figure  2-3  illustrates  that  for  large 
separations,  |5x  —  <52|,  /„(£)  essentially  reduces  to  the  linear  superposition  of  two 
solitons  with  parameters  0  and  0.  As  the  relative  separation  decreases,  the  multi¬ 
plicative  cross  term  becomes  significant,  and  the  solitons  interact  nonlinearly.  This 
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asymptotic  behavior  can  also  be  evidenced  analytically, 


77? 

fn{t )  =  — /?isec h2(pin- pi{t-61)±<l)/2)  (2.9) 

+  ^/3fsech2(p2™  -  fait  -  82)  T  0/2),  t  ±00, 

ao 

where  each  component  soliton  experiences  a  net  displacement  4>  from  the  nonlinear 
interaction. 

The  Toda  lattice  also  admits  periodic  solutions  which  can  be  written  in  terms  of 
Jacobian  elliptic  functions  dn(-)  and  sn(-).  These  solutions  can  be  expressed 

U(t)  =  {dn2  [2  (s  ±  vt)  *]-§}-  (2-10) 

for  wavelength  A  and  frequency  u,  with 

fnh  /  F\  — 

~  \j~  \  sn~2(2K/X)  —  1  +  — )  .  (2.11) 

For  further  description  of  Jacobian  elliptic  functions,  the  reader  is  referred  to  [2] 
and  [69]. 

An  interesting  observation  can  be  made  when  the  Toda  lattice  equations  are  writ¬ 
ten  in  terms  of  the  forces, 

^1u{1+t]  =  i.{u+i  ~ 2/" + (2-12) 

If  the  substitution 

fn(t)  =  (2-13) 

is  made  into  Eq.  (2.12),  then  the  lattice  equations  become, 

^  (^n  =  fin  4*n-l$n+l-  (2-14) 

In  view  of  the  Teager  energy  operator  introduced  by  Kaiser  in  [35],  the  left  hand  side 
of  Eq.  (2.14)  is  the  Teager  instantaneous- time  energy  at  the  node  n,  and  the  right 
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hand  side  is  the  Teager  instantaneous-space  energy  at  time  t.  In  this  form,  we  may 
view  solutions  to  Eq.  (2.14)  as  propagating  waveforms  that  have  equal  Teager  energy 
as  calculated  in  time  and  in  space,  a  relationship  also  observed  by  Kaiser  [36]. 


2.2  The  Inverse  Scattering  Transform 

Perhaps  the  most  significant  discovery  in  soliton  theory  was  that  under  a  rather  gen¬ 
eral  set  of  conditions,  certain  nonlinear  evolution  equations  such  as  KdV  or  the  Toda 
lattice  could  be  solved  analytically.  That  is,  given  an  initial  condition  of  the  system, 
the  solution  can  be  explicitly  determined  for  all  time.  This  discovery  eventually  led  to 
a  theory  of  solvable  nonlinear  dynamic  systems  to  which  KdV  and  Toda  belong  along 
with  many  other  nonlinear  systems  which  exhibit  soliton  behavior.  Since  much  of 
inverse  scattering  theory  is  beyond  the  scope  of  this  thesis,  we  will  only  present  some 
of  the  basic  elements  of  the  theory  which  we  will  exploit  in  later  chapters  and  refer 
the  interested  reader  to  the  comprehensive  treatment  given  in  the  text  by  Ablowitz 
and  Clarkson  [1]. 

The  nonlinear  systems  which  have  been  solved  by  inverse  scattering  belong  to  a 
class  of  systems  called  conservative  Hamiltonian  systems.  These  are  state-variable 
systems  where  the  dynamics  can  be  written  explicitly  in  terms  of  the  gradient  of  an 
energy  function,  or  the  Hamiltonian  of  the  system,  H(x).  For  example,  if  the  state 
of  the  system  were  given  by  x,  and  the  Hamiltonian  were  'H(x),  then  the  dynamics 
could  be  written  in  the  form 


0  I 
-I  0 


Vft(x). 


(2.15) 


For  conservative  Hamiltonian  systems,  'H(x)  will  not  be  a  function  of  time.  A  rigorous 
treatment  of  Hamiltonian  systems  can  be  found  in  the  text  [41].  Conditions  under 
which  conservative  Hamiltonian  systems  become  solvable  can  be  found  in  [3].  For  the 
nonlinear  systems  we  discuss  in  this  thesis,  an  integral  component  of  their  solution  lies 
in  the  ability  to  write  the  dynamics  of  the  system  implicitly  in  terms  of  an  operator 
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differential  equation  of  the  form, 


=  B(t)m  -  L(t)B(t), 


(2.16) 


where  L{t)  is  a  symmetric  linear  operator,  B(t)  is  an  anti-symmetric  linear  operator, 
and  both  L(t )  and  B(t)  depend  explicitly  on  the  state  of  the  system  x(£).  Lax  has 
shown  [42]  that  when  the  dynamics  can  be  written  in  the  form  of  Eq.  (2.16),  then 
the  eigenvalues  of  the  operator  L(t)  are  time-invariant,  i.e.  A  =  0.  The  ability  to 
write  the  dynamics  implicitly  in  the  form  of  Eq.  (2.16),  coupled  with  the  structure  of 
conservative  Hamiltonian  systems  is  sufficient  to  guarantee  the  solvability  of  each  of 
the  soliton  systems  that  will  be  discussed  in  this  thesis. 

Using  the  Toda  lattice  as  an  example,  the  operator  L  would  be  the  symmetric 

.  (2-17) 

K  =  j,  (2.18) 


matrix 


L  = 


•  Q-n—l 
®n— 1 

On 


where 


n  —  -p(Vn~yn+ 1)/2 
un  —  £  ’ 


for  mass  positions  yn  in  a  solution  to  Eq.  (2.2).  When  B  is  given  by  the  anti-symmetric 
matrix 


B  = 


•  on_  i 

On—  1  0  CZn 

On 


(2.19) 


then  Eq.  (2.16)  implicitly  contains  the  Toda  lattice  equations.  Although  each  of  the 
entries  of  L(t)  evolve  with  the  state  of  a  solution  to  the  Toda  lattice,  the  eigenvalues 
of  L(t )  remain  constant. 

If  we  assume  that  the  motion  on  the  lattice  is  confined  to  lie  within  a  finite  region 
of  the  lattice,  i.e.  the  lattice  is  at  rest  for  |n|  — >  oo,  then  the  spectrum  of  eigenvalues 
for  the  matrix  L{t)  can  be  separated  into  two  sets.  As  depicted  in  Fig.  2-4,  there 
is  a  continuum  of  eigenvalues  A  e  [—1,1]  and  a  discrete  set  of  eigenvalues  for  which 
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Figure  2-4:  Spectrum  of  eigenvalues  of  the  matrix  L(t)  for  the  Toda  lattice. 

|Afc|  >  1.  For  any  solution  to  the  infinite-length  Toda  lattice,  the  continuous  set  of 
eigenvalues  will  all  be  present.  However  the  discrete  eigenvalues  will  only  exist  if  there 
are  solitons  present  in  the  lattice  —  one  discrete  eigenvalue  for  each  soliton  excited. 
This  separation  of  eigenvalues  of  L(t)  into  discrete  and  continuous  components  is 
common  to  all  of  the  nonlinear  systems  solved  with  inverse  scattering. 

The  inverse  scattering  method  of  solution  for  soliton  systems  is  schematically 
similar  to  methods  used  to  solve  linear  evolution  equations.  For  example,  consider  a 
linear  evolution  equation  of  the  form 

^y(x,t)  =  £(t,x,y),  (2.20) 

where  £  is  a  linear  function  of  time  t,  space  x,  the  state  y  and  its  derivatives.  Given 
an  initial  condition  of  the  system,  y(x,  0),  a  standard  technique  for  solving  for  y(x,  t ) 
employs  Fourier  methods.  By  decomposing  the  initial  condition  into  a  superposition 
of  simple  harmonic  waves  of  the  form, 

1  r°° 

y(x,  0)  =  —  /  Y(k,0)ejkxdk,  (2.21) 

each  of  the  component  harmonic  waves  can  be  independently  propagated  using  the 
dispersion  relation  defined  by  £.  Given  the  Fourier  decomposition  of  the  state  at  time 
t,  the  harmonic  waves  can  then  be  recombined  to  produce  the  state  of  the  system 
y(x,t).  This  process  is  depicted  schematically  in  Fig.  2-5. 

An  outline  of  the  inverse  scattering  method  for  soliton  systems  is  similar.  Given 
an  initial  condition  for  the  nonlinear  system,  y(x,  0),  the  eigenvalues  A  and  eigen¬ 
functions  tp(x,  0)  of  the  linear  operator  L(0)  can  be  obtained.  This  step  is  often 
called  forward  scattering  by  analogy  to  quantum  mechanical  scattering.  To  obtain 
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»(*,0)  — >  F.T. 


k,  Y(k,  0) 


Jj.  Evolve  simply  in  time 
jj,  e-Mk)t 

y(x,t )  <—I.F.T.±—  k,Y(k,  t) 

Figure  2-5:  Schematic  solution  to  linear  evolution  equations. 
y(x,  0)  — >  F.S. — X,ip(x,0) 

1J.  Evolve  simply  in  time 

y{x,t)  < — I.S.  < —  A ,ip(x,t) 

Figure  2-6:  Schematic  solution  to  soliton  equations. 

the  eigenvalues  and  eigenfunctions  at  a  point  in  time  t,  all  that  is  needed  is  the 
time  evolution  of  the  eigenfunctions,  since  the  eigenvalues  do  not  change  with  time. 
For  these  soliton  systems,  the  eigenfunctions  evolve  simply  in  time,  according  to  lin¬ 
ear  differential  equations.  Given  the  eigenvalue-eigenfunction  decomposition  of  L(t), 
through  a  process  called  inverse  scattering ,  the  state  of  the  system  y(x,t)  can  be 
completely  reconstructed.  This  process  is  depicted  in  Fig.  2-6  in  a  similar  fashion  to 
the  linear  solution  process. 

Once  again,  we  turn  to  the  Toda  lattice  as  an  example  to  illustrate  the  inverse 
scattering  method  in  slightly  more  detail.  Although  the  discussion  that  follows  is 
specific  to  the  Toda  lattice,  the  general  approach  and  much  of  the  theory  applies 
directly  to  all  soliton  systems. 

Once  a  pair  of  linear  operators  for  the  nonlinear  evolution  equation  is  known, 
the  dynamics  of  the  nonlinear  system  can  be  inferred  through  the  behavior  of  the 
eigenvalues  and  eigenfunctions  of  L(t).  Here  we  consider  the  eigenvalue-eigenfunction 
decomposition 

L(t)ip{t)  =  A^(t),  (2.22) 


25 


the  n-th  row  of  which  corresponds  to 


(L{t)ip(t))n  =  an-i(t)ip(n  -l,t)  +  bn(t)%l>{n,t)  +  an{t)ip(n  4-  M)  =  \ip(n,t).  (2.23) 
The  eigenfunctions  satisfy, 

*?&$■  =  B(t)1,(n,t),  (2.24) 

the  n-th  row  of  which  yields, 

=  a -  1,  t)  -  On(t)ip(n  +  1,  i).  (2.25) 

at 

Assuming  that  the  lattice  is  at  rest  (an  =  1/2  and  bn  =  0)  for  |n|  1,  the 

eigenvalues  of  L(t)  consist  of  a  continuum  A  €  [-1,1],  and  a  finite  set  of  discrete 
eigenvalues  | A^j  >  1.  The  corresponding  eigenfunctions  can  be  specified  in  terms  of 
their  asymptotic  behavior  in  the  regions  where  the  motion  on  the  lattice  vanishes. 
Specifically,  the  eigenfunctions  (k(n,  t)  corresponding  to  the  discrete  eigenvalues  A* 
are  specified  as 

Ck(n,t)  «  ck(t)z£,  n— >oo,  (2.26) 

where  \zk\  <  1  and  Afc  =  {zk  +  z^1)/ 2,  with 

OO 

E  (C*(M))2  =  1-  (2.27) 

n=— oo 

The  eigenfunctions  t )  corresponding  to  the  eigenvalues  A  in  the  continuum  are 
specified  as 

ip(n,  t)  «  z~n  +  R(z,  t)zn,  n  — »  oo,  (2.28) 

z~n 

il)(n,t)  «  — - — —  n  — >  — oo,  (2.29) 

for  A  =  (z  4-  z-1)/2. 

We  can  now  solve  the  initial  value  problem  for  the  Toda  lattice.  Given  y„(0), 
the  initial  eigenvalue-eigenfunction  decomposition  can  be  found  for  L(0).  The  collec- 
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tion  of  discrete  eigenvalues  and  the  asymptotic  description  of  both  the  discrete  and 
continuous  eigenfunctions  are  often  referred  to  as  the  nonlinear  spectrum  for  soliton 
systems.  Here  our  nonlinear  spectrum  consists  of 


S{ A,  0)  =  [{(zk,  ck( 0))  :  1  <  k  <  N}  ,  R(z,  0),  a(z,  0)] .  (2.30) 

From  (2.25),  the  time  evolution  of  the  nonlinear  spectrum  can  be  expressed  as 

cfe(t)  =  c*(  0)e^\  (3k  =  (z^-zk)/2  (2.31) 

a(z,t )  =  a(z,  0),  (2.32) 

R(z,t )  =  R(z,  0)e2/3t,  P  =  (z~1-z)/2  (2.33) 


which  yields  the  nonlinear  spectrum  at  time  t,  S(X,t).  The  inverse  scattering  problem 
is  then  to  reconstruct  the  state  of  the  system,  an(t),  and  bn(t)  at  time  t,  given  S( A,  t). 
From  the  nonlinear  spectrum,  we  can  define  the  following  function, 

F(n ,  t)  =  Xj  4(1)4  +  ^  f  R(z,  t)zn~ldz.  (2.34) 

We  then  seek  a  solution  K(n,  m,  t )  to  the  Gel’fand-Levitan-Marchenko  (GLM)  equa¬ 
tion, 


K,(n,m,t)  +  F(n  +  m,t)  +  X]  K,(n,n',t)F(n' +  m,t)  =  0,  m  >  n,  (2.35) 

n'  =n+ 1 


where, 


and 


K(n,n,t)  2  =  1  +  F(2n,t)+  XI  K(n,n',t)F(n'  +  m,t),  (2.36) 

n'=n+l 


«(n,  m,  t )  = 


i^(n,  m,  t) 


(2.37) 


K(n,n,t )  ’ 

The  GLM  equation  form  a  linear  discrete- integral  equation  in  the  terms  «(n,  m,  t). 
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The  solution  yn(t )  can  then  be  found  from 


fln(t) 

bn(t) 


K(n  +  1  ,n  +  l,i) 

2  K(n,n,t) 

K(n,n  +  l,t )  K(n  —  l,n,  t) 

2  K(n,n,t)  2K(n  —  l,n  —  l,t) 


(2.38) 


The  forces  on  the  springs,  /„(t),  can  be  obtained  directly 


/«(*)  = 


AT(n,  n,  t) 
K(n  —  l,n  —  l,t) 


(2.39) 


For  a  large  class  of  soliton  systems,  the  inverse  scattering  method  generally  in¬ 
volves  solving  either  a  linear  integral  equation  or  a  linear  discrete-integral  equation. 
The  general  form  of  the  GLM  equation  follows  that  for  the  Toda  lattice.  Although 
the  equation  is  linear,  finding  the  solution  of  the  GLM  is  often  very  difficult  in  prac¬ 
tice.  However,  when  the  reflection  coefficient,  R(z ),  is  zero,  then  Eq.  (2.35)  reduces  a 
set  of  simultaneous  linear  equations.  This  corresponds  to  a  solution  made  up  of  pure 
solitons,  one  soliton  for  each  discrete  eigenvalue.  Examples  of  this  method  for  the 
Toda  lattice  are  carried  out  for  a  single  soliton  and  a  multi-soliton  solution  as  well  as 
a  solution  with  a  single  pole  reflection  coefficient  in  Appendices  2.A-2.C. 


2.3  Other  Systems  Exhibiting  Solitons 

Since  the  discovery  of  the  inverse  scattering  method  for  the  solution  to  KdV,  there 
has  been  a  large  class  of  nonlinear  wave  equations,  both  continuous  and  discrete,  for 
which  similar  solution  methods  have  been  obtained.  In  most  cases,  solutions  to  these 
equations  can  be  constructed  from  a  nonlinear  superposition  of  soliton  solutions.  For 
a  comprehensive  study  of  inverse  scattering  and  equations  solvable  by  this  method, 
the  reader  is  referred  to  the  text  by  Ablowitz  and  Clarkson  [1].  We  briefly  mention  a 
few  to  indicate  the  rich  class  of  such  equations. 
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1.  The  sine-Gordon  equation 


ipxx  -  fpu  =  Sin(^).  (2.40) 

The  sine-Gordon  (SG)  equation  has  been  used  as  a  model  for  a  variety  of  phys¬ 
ical  phenomena  including  the  propagation  of  crystal  dislocation  [57]  and  the 
propagation  of  magnetic  flux  along  a  Josephson  strip  line  [56].  A  simple  phys¬ 
ical  model  can  be  constructed  by  hanging  pendula  from  an  elastic  line  -  the 
angular  displacement  of  the  pendula  are  approximately  governed  by  the  SG 
equation.  A  soliton  solution  is  given  by 

ip{x,t)  =  4  tan-1  e±('/i-u2)?  (2-41) 

which  is  often  called  a  “kink”  soliton  since  it  corresponds  to  a  change  in  phase. 
An  inverse  scattering  solution  for  the  SG  equation  has  been  developed  by 
Ablowitz,  Kaup,  Newell,  and  Segur  [1]. 

2.  The  Nonlinear  Schrodinger  Equation 

iut  +  uxx  ±  k\u\2u  =  0.  (2.42) 

The  most  widely-known  application  of  the  nonlinear  Schrodinger  (NLS)  equa¬ 
tion  has  been  to  describe  the  propagation  of  lightwave  pulses  along  nonlinear 
fiber-optic  cables  [23].  Other  applications  include  the  propagation  of  heat  in 
solids,  and  Langmuir  waves  in  plasma  [57] .  The  NLS  equation  admits  envelope 
solitons  of  the  form 

u  =  u0sech  {y/k2u0(x  —  vet)j  eI(^/2)(i -vct),  (2.43) 

where  ve  and  vc  are  the  envelope  and  carrier  velocities  respectively  [57].  The 
inverse  scattering  solution  to  the  NLS  equation  was  given  by  Zakharov  and 
Shabat  [79]. 
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3.  The  Discrete-KdV  Equation 


un  =  eUn_1  —  eVn+1 .  (2.44) 

The  discrete-KdV  (dKdV)  equation  was  discussed  by  Kac  and  Van  Moerbeke 
in  relation  to  the  Toda  lattice.  Manakov  also  studied  the  dKdV  lattice  for  its 
relation  to  a  model  for  Langmuir  oscillations  in  plasma.  An  inverse  scattering 
theory  for  dKdV  was  developed  by  each  independently  in  [33]  and  [45].  Manakov 
presents  the  single  soliton  solution  for  Nn  =  eUn , 

_  cosh (77(71  -  sp  -  2))  cosh(r/(n  -  zp  +  1)) 
n{)  cosh(r](n-xo-l))cosh(rj(n-xo))  ’  J 

where 

*.(()=  *„(0)+* BMMt  (2.46) 

4.  Cellular  Automata 

A  number  of  fully-discrete  dynamic  systems,  or  cellular  automata  (CA),  have 
been  shown  to  possess  solutions  with  many  of  the  properties  of  solitons.  Ablowitz 
et  al.  have  shown  several  automata  that  possess  soliton-like  solutions  which  are 
remarkably  similar  to  the  solitons  of  KdV  or  the  Toda  lattice  [1,  52]  as  shown  in 
Fig.  2-7.  Park,  Steiglitz  and  Thurston  have  also  observed  soliton-like  behavior 
in  automata,  and  even  suggested  possible  computational  applications  of  such 
systems  [53,  63,  62].  Takahashi  has  presented  a  CA  similar  to  the  automata 
of  Ablowitz  et  al.,  which  also  not  only  possesses  soliton-like  solutions,  but  also 
has  an  infinite  number  of  conserved  quantities,  a  property  shared  by  solvable 
nonlinear  systems  [66]. 

2.4  Soliton  Hierarchies 

In  this  section  we  mention  a  few  of  the  techniques  that  can  be  used  to  construct 
families  of  solvable  nonlinear  evolution  equations.  Since  we  shall  see  that  the  soliton 
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Figure  2-7:  Soliton-like  solutions  to  the  automata  of  Ablowitz  et  al. 

solutions  of  such  systems  are  potentially  useful  for  communication,  there  may  be  a 
practical  use  for  the  generation  of  increasingly  complex  soliton  systems. 

Many  solvable  nonlinear  evolution  equations,  including  the  Korteweg-de  Vries  and 
nonlinear  Schrodinger  equations,  can  be  written  in  the  form 

^  =  *(«),  (2-47) 

where  K(u)  may  be  a  function  of  u  and  its  ^-derivatives.  In  [42],  Lax  demonstrated 
conditions  under  which  the  eigenvalues  of  a  self-adjoint  linear  operator  L,  parameter¬ 
ized  by  u,  would  remain  invariant  under  time  evolution  of  u(x,t)  and  hence,  of  L(t). 
This  is  the  so-called  Lax  equation,  L  =  BL  —  LB,  where  B  is  an  anti-symmetric 
operator. 

In  many  cases  this  leads  directly  to  an  inverse  scattering  framework,  from  which 
the  nonlinear  evolution  equations  can  be  solved.  Unfortunately,  given  a  nonlinear 
evolution  equation,  there  is  no  constructive  means  known  to  generate  a  Lax  pair 
for  the  system  or  even  to  tell  if  one  exists.  Lax  does  indicate  in  [42]  how  given  an 
operator,  L ,  an  associated  operator  B  may  be  constructed  such  that  the  resulting 
nonlinear  evolution  equation  is  nontrivial.  Although  such  evolution  equations  have 
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often  been  discounted  for  having  no  physical  analog,  such  equations  clearly  can  be 
used  in  the  signal  processing  contexts  discussed  in  this  thesis  either  through  numerical 
integration  or  with  the  aid  of  high-speed  nonlinear  circuitry. 

A  similar  technique  was  applied  to  the  matrix  L(t)  for  the  Toda  lattice,  viz., 

L[n,  ?7i]  —  0,n—l^n,m+l  4-  bn8n>m  "h  ®n^n+l,mi  (2.48) 

by  Moser  in  [48].  Here,  the  evolution  equation  will  be  defined  in  terms  of  L\  hence 
we  require  BL  —  LB  to  be  a  tridiagonal  symmetric  matrix  yielding  the  evolution  of 
the  parameters,  an(t )  and  bn(t).  A  Toda  lattice  hierarchy  can  be  found  by  specifying 
an  algorithm  for  selecting  the  antisymmetric  matrix,  B,  though  the  calculations  are 
cumbersome.  The  first  of  such  matrices  would  be  the  simple  matrix, 

5i[tI,  m]  Cn—lbn,m+l  Cra^n+1,7715  (2.49) 

which  leads  to  the  following  set  of  evolution  equations, 

On.  =  On(6„  ^n+l)  (2.50) 

bn  =  2c(a£_1-  al). 

With  c  =  1,  Eqs.  (2.50)  are  the  Toda  lattice  equations. 

The  next  in  the  Toda  hierarchy  can  be  found  by  allowing  the  matrix  B  to  be  a 
penta-diagonal  matrix, 

B2[n,  m]  —  dn—2&n—2,m  Cn—\bn—l,m  Cn^n+ l,m  dn8n+2,m-  (2.51) 

This  leads  to  the  following  set  of  evolution  equations, 

On  =  On(on_j  —  an+1  +  b^  —  5^+1)  (2.52) 

bn  i(bn— x  H-  bn )  an(bn  -I-  fen+i))- 
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This  higher  order  lattice  equation  contains  the  discrete-KdV  equation  as  a  special 
case,  when  bn  =  0,  reduce  to 

On  —  on(on_1  (2.53) 

Defining  an  as 

an  =  \eUn/2,  (2.54) 

leads  to  the  discrete-KdV  equation, 

un  =  (e""-1  -  eUn+1) .  (2.55) 

This  system  has  been  studied  by  Moser  in  [48].  The  inverse  scattering  solution  for 

this  system  was  given  by  Manakov  for  the  infinite  line  case  in  [45],  and  Kac  and  van 
Moerbeke  in  [33]  and  [34]  for  both  the  semi-infinite  and  periodic  cases.  Though  little 
emphasis  has  been  placed  on  this  system,  since  it  has  no  physical  analog  (though 
Manakov  suggests  a  relation  to  the  spectra  of  Langmuir  oscillations  in  plasma  [45]), 
realizations  of  the  discrete-KdV  system  may  implemented  either  numerically  or  with 
analog  circuits  making  the  system  of  potential  use  in  the  signal  processing  contexts 
discussed  in  this  thesis.  In  Chapter  3,  we  show  both  a  single  circuit  implementation 
of  this  system,  along  with  an  implementation  using  two  Toda  lattice  circuits  based 
on  a  relationship  suggested  by  Kac  and  van  Moerbeke  [33]  and  [68]. 

As  the  number  of  nonlinear  evolution  equations  solved  using  the  inverse  scattering 
method  grew,  more  methods  for  generating  hierarchies  of  nonlinear  equations  solvable 
by  inverse  scattering  were  developed.  Ablowitz,  Kaup,  Newell,  and  Segur,  developed 
a  framework  for  generating  large  classes  of  nonlinear  evolution  equations  solvable 
by  inverse  scattering  [1].  This  method  is  a  generalization  of  the  ideas  of  Lax,  and 
centers  around  a  more  general  form  of  eigenvalue  problem.  For  further  information, 
the  reader  is  referred  to  [1],  and  [49]  and  the  list  of  references  therein. 
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2. A  A  Soliton  Solution 


The  discrete  integral  equation,  or  GLM,  is  often  difficult  to  solve,  except  under  certain 
interesting  conditions.  Specifically,  for  the  case  when  the  reflection  coefficient  vanishes 
identically,  R(z,  t)  =  0,  the  GLM  reduces  to  a  set  of  linear  equations  which  can  be 
solved  by  simple  methods.  The  only  eigenfunctions  that  result  in  this  case  are  those 
corresponding  to  discrete  eigenvalues. 

We  consider  first  the  case  when  the  solution  has  no  contribution  from  the  contin¬ 
uous  eigenvalues,  R(z )  =  0,  and  has  only  a  single  term  from  the  discrete  spectrum, 
which  we  write, 

20  =  ±e-7-  (2.56) 

That  zQ  can  be  taken  as  real  comes  from  the  symmetric  assumption  on  L(t).  We  com¬ 
plete  the  specification  of  the  solution  by  selecting  the  initial  value  of  the  normalization 
coefficient,  Co(0).  By  forming  the  function, 

F(m,t)  =  clzg,  (2.57) 

where  c0  =  co(0)eA<,  and  (3q  =  ±sinh(7). 

The  GLM  equation  can  now  be  written  as 

°0 

k(ti,  m)  +  c$z£+m  +  c%z™  ^2  k(ti,  ti')zq  =  0.  (2.58) 

re'=n+l 

Following  [15],  (2.58)  can  be  solved  by  assuming  a  solution  of  the  form, 


leading  to 


where 


k{u,  m)  —  AnCoZ™, 

A  -Co*# 

n  1  +  e^02(n+1)  ’ 

es  =  c0(t)  _  es0+p0t 


(2.59) 


(2.60) 


(2.61) 
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Substituting  (2.59)  into  (2.38),  after  some  manipulations,  we  have 

e-(9n-«n- 1)  —  i  —  ffiseda2  (jn  —  0q  t  —  6 o).  (2.62) 

This  corresponds  to  a  single  soliton  solution  whose  direction  of  propagation  is  deter¬ 
mined  by  the  sign  of  /30,  and  hence  the  sign  of  z0. 


2.B  Multi-soliton  Solutions 


Once  again,  we  consider  the  case  where  R(z )  =  0,  and  for  the  discrete  spectrum,  we 
allow  N  eigenvalues  corresponding  to  z\, . . . ,  z^.  In  this  case  the  function  F(m)  in 
the  GLM  is  given  by, 

N 

-  (2.63) 


F(m)=E<& 

J=1 


This  results  in  the  GLM 


N  oo  JV 

K(n,m)  =  '£c2jz]+m+  £  «(^')E^,+m  =  0. 

j= 1  n'=n+l  3= 1 


(2.64) 


By  assuming  a  solution  of  the  form 


N 


K(n,  m)  =  ]T  Aj,nCjZ™, 

j=l 


(2.65) 


and  substituting  into  (2.64),  for  each  n,  after  some  algebra,  we  have  a  set  of  N  linear 
equations  for  the  N  unknowns  Aj,n, 


£  +  ^TZ~)  =  -<**?,  0  =  1.  • .  . .  JV). 

i= 1  \  1  / 

This  can  be  written  in  matrix  form, 

B(n) 


(2.66) 


A\yn 

~CiZ1 

-cNz% 

(2.67) 
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which  lead  to  a  direct  formula  for  the  soliton  solution, 


.(*-*-1)  _  i  =  |B(n)llg(n-2)l  _ 
\B(n  —  1)|2 


2.C  Single  Pole  Reflection  Coefficient 


(2.68) 


We  now  consider  an  example  when  the  reflection  coefficient  contains  a  single  pole  in 
the  z-plane,  at  z  =  7,  e.g, 

R(Z)  =  T^F V  (2-69) 

where  7  >  0,  and  c  are  constants.  Since  the  discrete  spectrum  contains  zeros  of  a(z), 
there  must  be  at  least  one  discrete  eigenvalue  corresponding  to  the  eigenfunction 


(j>(z)  ~  07",  n  — »  00.  (2.70) 

From  the  definition  of  the  function  F(m ),  we  have 

F(m)  =  ^  +  ±ijT^_z^dz.  (2.7!) 

By  the  Cauchy  residue  theorem,  the  integral  on  the  right  hand  side  evaluates  to 

=  (2-72) 

where  u[m\  is  the  discrete  unit-step  function.  This  results  in 

F(m)  =  cymu[— m  —  1].  (2.73) 

The  resulting  GLM  for  this  nonlinear  spectrum  is  given  by 

OO 

K(n,m)+C7n+mu[-(n+m)-l]+  ^  K(n,n')cyn'+mu{-(n' +m)-l\  =  0,  m  >  n. 

n’  —n+1 

(2.74) 
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For  the  diagonal  terms,  n  =  m,  we  have 


OO 

K(n,n)~ 2  =  1  +  C72nu[—  (2n)  —  1]  +  Y  /c(n,  n')c7n,+nu[— (n'  +  n)  —  1].  (2.75) 

n'=n+l 


From  (2.74)  we  see  immediately  that 


k(u,  m)  =  0,  n  +  m>0. 


(2.76) 


For  n  +  m  <  0,  we  have 


— TO—  1 


k{u,  m)  +  C7n+m  +  Y  K(n>  n')c')n  +m  =  0,  m  >  n.  (2.77) 

n'=n+ 1 


One  method  of  solving  (2.77),  is  to  take  the  first  difference  with  respect  to  m, 
resulting  in 


k(ti,  m  +  1)  —  K(n ,  m)  +  c'yn+m(/y  —  1)  + 

-(m+l)— 1  —  m— 1 

Y  K{n,n')cjn'+m+1  -  Y  <n,ri)c 7n'+m  )  =  0, 

n'=n+l  n'=n+ 1 


which  leads  to 


K(n,m  +  1)  —  K(n,  m)  +  c'yn+m(j  —  1)  + 


-TO— 1 


Y  K(ni  n,)cryn'+m  )  (7  -  1)  -  >c(n,  -m  —  l)c  =  0. 

Vn'=n+1 


By  substitution,  we  arrive  at 


K(n,m+l)-«;(n,m)+C7n+m(7-l)-(-(-«;(n, m)-c7n+m)(7-l)-K(n,  -m-l)c  =  0, 

(2.78) 

which  is  solved  if  c  =  1  -  7  and  /c(n,m)  =  (7  —  l)/7-  Using  these  values  in  (2.75), 
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yields 


K(n ,  n) 
K(n,  m) 


0,  n  >  0 

\Z^T>  n  <  ^ 


0,  n  +  m  >  0 

n  +  m  <  0. 


(2.79) 

(2.80) 


The  initial  condition  that  led  to  this  one  pole  nonlinear  spectrum  is  given  by 


1  ^7(27-1)- 7 

an  =  2  +  ^ - ^ - <5[n+l] 

,  n  ,  (l-7)(27-l)„_1 
On  -  0  H - — - o|n], 


(2.81) 

(2.82) 


which  is  a  localized  perturbation  to  the  lattice  at  rest.  Even  for  small  disturbances, 
|7|  <C  1,  we  see  that  the  localized  perturbation  gives  rise  to  both  continuous  and  dis¬ 
crete  components  of  the  spectrum.  In  [16],  Flaschka  demonstrates  that  any  localized 
disturbance  must  give  rise  to  both  continuous  and  discrete  components. 
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Chapter  3 


New  Electrical  Analogs  for  Soliton 
Systems 

Since  soliton  theory  has  its  roots  in  mathematical  physics,  most  of  the  systems  stud¬ 
ied  in  the  literature  have  at  least  some  foundation  in  physical  systems  in  nature.  For 
example,  the  KdV  equation  was  developed  to  describe  the  dynamics  of  small  ampli¬ 
tude  surface  water  waves,  and  a  large  class  of  weakly  nonlinear  systems  have  been 
demonstrated  to  reduce  to  KdV  in  a  variety  of  physical  applications.  Such  systems 
range  from  ion-acoustic  waves  in  plasma  [78],  to  pressure  waves  in  liquid  gas  bubble 
mixtures  [57].  As  a  result,  the  predominant  purpose  of  soliton  research  has  been  to 
explain  physical  properties  of  natural  systems.  In  addition,  there  are  several  exam¬ 
ples  of  man-made  media  that  have  been  designed  to  support  soliton  solutions  and 
thus  exploit  their  robust  propagation.  The  use  of  optical  fiber  solitons  for  telecom¬ 
munications  and  of  Josephson  junctions  for  volatile  memory  cells  are  two  practical 
examples  [56,  57]. 

Whether  its  goal  has  been  to  explain  natural  phenomena  or  to  support  propa¬ 
gating  solitons,  this  research  has  largely  focussed  on  the  properties  of  propagating 
solitons  through  these  nonlinear  systems.  In  this  thesis,  we  take  a  decidedly  different 
perspective.  We  view  solitons  as  signals  and  consider  exploiting  some  of  their  rich 
signal  properties  in  a  signal  processing  or  communication  context.  This  perspective 
is  illustrated  graphically  in  Fig.  3-1,  where  a  signal  containing  two  solitons  is  shown 
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Figure  3-1:  Two-soli  ton  signal  processing  by  a  soliton  system. 

as  an  input  to  a  soliton  system  which  can  either  combine  or  separate  the  compo¬ 
nent  solitons  according  to  the  evolution  equations.  From  the  “solitons-as-signals” 
perspective,  the  corresponding  nonlinear  evolution  equations  can  be  viewed  as  spe¬ 
cial  purpose  signal  processors  that  are  naturally  suited  to  performing  complex  signal 
processing  tasks  such  as  signal  separation  or  sorting.  As  we  shall  see,  these  systems 
also  form  an  effective  means  of  generating  soliton  signals. 

Although  many  soliton  systems  are  described  by  partial  differential  equations,  the 
concept  of  performing  signal  processing  operations  with  partial  differential  equations, 
or  wave  equations,  is  not  new.  For  example,  signal  generation  and  detection  in  pulse 
compression  radar  systems  can  be  accomplished  using  a  dispersive  delay  line  [61]. 
Such  systems  can  be  constructed  either  using  bulk  material,  like  a  surface  acoustic 
wave  (SAW)  device,  or  with  a  lumped  element  line.  In  either  case,  a  signal  is  induced 
into  the  device  and  allowed  to  propagate  according  to  the  wave  equations  that  describe 
the  device.  Once  the  signal  has  propagated  along  the  device  and  the  appropriate  delay 
or  other  processing  has  been  accomplished,  the  signal  can  be  extracted  from  the 
medium.  SAW  devices  can  also  be  used  for  a  variety  of  signal  processing  applications 
including  chirp  Z-transforms  or  digital  filtering  [11]. 

In  this  chapter,  we  study  circuit  implementations  of  two  of  the  nonlinear  evolution 
equations  discussed  in  the  previous  chapter.  The  first  of  these  is  a  nonlinear  LC 
ladder  network  developed  by  Hirota  and  Suzuki  [29].  Although  this  circuit  serves 
as  a  useful  electrical  model  for  the  Toda  lattice,  it  is  difficult  to  implement  using 
standard  components.  We  present  a  new  circuit  model  for  the  Toda  lattice  based  on 
a  precise  electrical  analog  of  the  exponential  spring  mass  system.  This  circuit  has 
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Figure  3-2:  Nonlinear  LC  ladder  circuit  of  Hirota  and  Suzuki. 

been  implemented  in  hardware  using  standard  components  and  appears  to  be  the  first 
such  circuit  to  display  true  soliton  behavior.  We  also  present  a  new  circuit  model  for 
the  discrete-KdV  equation,  for  which  there  is  no  prior  electrical  or  mechanical  analog. 

3.1  Toda  Circuit  Model  of  Hirota  and  Suzuki 

Motivated  by  the  work  of  Toda  on  the  exponential  lattice,  the  nonlinear  LC  ladder 
network  implementation  shown  in  Fig.  3-2  was  given  by  Hirota  and  Suzuki  in  [29]. 
Rather  than  a  direct  analogy  to  the  Toda  lattice,  the  authors  derived  the  functional 
form  of  the  capacitance  required  for  the  LC  line  to  be  equivalent.  The  resulting 
network  equations  are  given  by 

i1  +  Tr)  =  Zck(K-lW  -  2Ut)  +  Wt))’  (3'1) 

which  is  equivalent  to  the  Toda  lattice  equation  for  the  forces  on  the  nonlinear  springs 
given  in  Eq.  (2.12).  This  amounts  to  an  implicit  mapping  from  force  to  voltage, 
fn(t)  —>•  Vn(t).  The  capacitance  required  in  the  nonlinear  LC  ladder  is  of  the  form 

GWo 

Vq  +  V' 

where  Vo  and  Co  are  constants  representing  the  bias  voltage  and  the  nominal  capac¬ 
itance,  respectively.  In  their  implementation,  varactor  diodes  with  nonlinear  capaci¬ 
tance 

C{V)  «  27{V  -  H)-48pF,  (3.3) 
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where  V&  is  a  bias  voltage,  were  used  to  approximate  the  required  capacitance  of 
Eq.  (3.2). 

Although  the  varactor  diode  capacitance  can  be  biased  to  yield  a  fairly  good 
match  for  small  voltages,  for  larger  voltages,  the  deviations  from  the  ideal  capacitance 
becomes  apparent.  Moreover,  as  the  length  of  the  lattice  increases,  the  effects  on  any 
propagating  solitons  accumulate.  The  net  result  is  that  interaction  between  solitary 
waves  of  appreciable  amplitude  will  not  result  in  soliton  collisions;  rather  such  a 
collision  will  also  produce  a  nontrivial  amount  of  ripple  [29].  Also,  since  the  circuit 
is  only  accurate  for  small  voltages,  where  the  velocity  difference  between  solitons  is 
small,  large  numbers  of  nodes  are  required  to  bring  about  collisions. 

After  publication  of  their  circuit  [29]  and  subsequent  publication  of  modulation 
experiments  using  the  circuit  [64, 65],  several  papers  have  appeared  in  the  literature  on 
a  variety  of  related  topics.  In  [39],  Kolosick  et  al.  analyze  a  similar  nonlinear  network. 
In  [32]  Islam,  Singh  and  Steiglitz  studied  the  effects  of  dissipation  on  the  propagation 
of  individual  solitons  as  well  as  the  interaction  of  pairs  of  solitons.  They  found 
that  dissipative  effects  led  to  a  decrease  in  amplitude  and  an  increase  in  the  width 
of  solitons  as  they  propagate  through  the  lattice.  These  findings  are  in  agreement 
with  the  numerical  work  of  Okada,  Watanabe  and  Tanaca  in  [50],  whose  studies 
showed  similar  effects  due  to  parameter  fluctuation  in  the  periodic  Toda  lattice.  In  [4] 
Ballantyne  et  al.  observed  the  Jacobian  elliptic  function  solutions  in  a  periodic  version 
of  the  nonlinear  LC  line.  Toda  also  demonstrated  properties  of  the  nonlinear  line  and 
illustrated  the  existence  of  modulated  solitons,  by  relating  the  lattice  to  the  nonlinear 
Schrodinger  equation  in  [68].  Finally,  Cho,  Wakita  and  Miyigawa  developed  a  similar 
nonlinear  network  as  an  equivalent  circuit  model  for  the  propagation  of  nonlinear 
surface  acoustic  waves  in  thin-bar  and  broad-plate  vibrations.  They  also  have  shown 
that  the  nonlinear  LC  network  is  an  accurate  model  for  a  metallic  grating  waveguide 
and  use  this  circuit  model  to  explain  certain  nonlinearities  observed  in  SAW  devices 
including  the  generation  of  acoustic  phase  conjugate  waves  [10]. 
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Figure  3-3:  Diode  ladder  network. 


3.2  Diode  Ladder  Circuit  Model  for  Toda  Lattice 

Although  the  nonlinear  ladder  network  realizations  of  the  Toda  lattice  retain  many 
of  the  properties  of  the  ideal  lattice,  as  suggested  in  Section.  3.1,  the  dynamics  of 
these  circuits  are  limited  to  a  small  range  of  voltages  and  therefore  their  applicability 
is  inherently  limited.  In  this  section,  we  present  a  new  circuit  model  that  accurately 
matches  the  Toda  lattice  and  is  a  direct  electrical  analog  of  the  nonlinear  spring  mass 
system.  If  voltages  un_i  and  vn  are  applied  to  the  terminals  of  a  junction  diode,  then 
the  current  through  the  device  is  accurately  modeled  by 

in  =  Is  (eiVn-l~Vn)/vt  -  l),  (3.4) 

where  Is  is  the  saturation  current  and  vt  is  the  thermal  voltage.  If  we  place  the 
diodes  in  a  ladder  configuration  as  shown  in  Fig.  3-3,  then  the  current  through  the 
n-th  shunt  impedance  is  given  by 


*n+ 1  —  Is 


(Vn-1-Vn)/Vt  _  e(Vn~Vn+l)/vt  ^ 


(3.5) 


In  analogy  to  (2.2),  we  see  that  if  the  shunt  impedance  has  the  following  voltage- 
current  relation 

^^-=a(in(t)-in+1(t)),  (3.6) 

then  the  governing  equations  become 


djVji(f)_  _  J  (  (v„-i(t)-vn(t))/vt  _  Jvn(t)-vn+l(t))/vt 
dt2  ~  al*  e 


(3.7) 
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Figure  3-4:  Double  capacitor  circuit  diagram, 
or, 

+  =^(i„-1(i)-2i„(t)+i„+1(«)),  (3.8) 

where  *i(£)  =  ijn(£).  These  are  equivalent  to  the  Toda  lattice  equations  with  a/m  = 
als  and  b  =  l/vt.  The  required  shunt  impedance  is  often  referred  to  as  a  dou¬ 
ble  capacitor,  which  can  be  realized  using  ideal  operational  amplifiers  in  the  gy- 
rator  circuit  shown  in  Fig.  3-4,  yielding  the  required  impedance  of  Zn  =  ajs 2  = 
R3/RiR2C2s2  [30,  59]. 

When  «in(£)  in  Fig.  3-3  is  of  the  form 


kn  (t)  =  Istt2  sech2(7t),  (3.9) 

7  =  SlyJl'Ctfvu 

a  single  soliton  is  induced  in  the  line  resulting  in 

in{t )  =  IsQ2sech2(pn  -  yt),  (3.10) 

where  Q  =  sinh(p).  Note  that  the  saturation  current  Is  may  be  absorbed  into  the 
parameter  Q,  yielding 

in{t)  =  ^2sech2(pn  -  /?r),  (3-11) 

where  /3  =  vX  sinh(p),  and  r  =  t^Ja/vt.  Since  Is  is  generally  on  the  order  of  pico- 
amps,  the  operating  range  of  the  circuit  can  be  on  the  order  of  milliamps  over  a  wide 


44 


range  of  values  of  the  soliton  wavenumber  p.  As  a  result,  the  diode  ladder  circuit 
model  is  very  accurate  over  a  large  range  of  soliton  wavenumbers,  and  is  significantly 
more  accurate  than  the  LC  circuit  of  Hirota  and  Suzuki. 

Solitons  of  the  form  of  Eq.  (3.11)  are  solutions  of  the  infinite-length  Toda  lattice 
equations.  In  practice,  a  finite-length  lattice  can  be  constructed  to  yield  soliton  solu¬ 
tions  if  the  diode  ladder  circuit  can  be  appropriately  terminated  to  limit  reflections. 
As  a  starting  point,  we  consider  the  termination  that  would  yield  no  reflections  for 
the  small  signal  model.  This  can  be  obtained  from  the  impedance  of  the  line  when 
the  diodes  are  replaced  with  their  equivalent  linearized  resistance  q  =  vtj id,  where 
id  is  the  current  in  the  linearized  diode.  This  results  in  an  impedance 


Z\n  — 


(3.12) 


For  typical  component  values,  a  «  10u.  If  R?q  is  taken  to  be  ut/25mA  =  1,  then  for 
frequencies  below  1  MHz,  a  load  impedance  consisting  of  a  lfi  resistor  and  a  0.3/iF 
capacitor  approximate  Eq.  (3.12)  well  and  yield  negligible  reflections  in  practice. 


3.3  Circuit  Simulation 

The  diode  lattice  has  been  simulated  using  realistic  component  models  in  the  circuit 
simulation  package  HSPICE  [47].  The  diodes  used  are  model  ln4148  with  a  saturation 
current  of  Is  «  .OlpA.  Setting  the  operation  range  of  the  circuit  to  produce  solitons 
on  the  order  of  10mA  yields  a  value  of  p  «  14.  To  fix  the  time  scale  of  the  circuit,  we 
set  the  pulse  width  of  a  soliton  to  approximately  bps,  which  leads  to 

JTs’  (313} 

or  a  «  10n.  The  resistor  values  in  the  double  capacitor  circuits  can  now  be  chosen  to 
prevent  saturation  of  the  operational  amplifiers.  By  calculating  the  transfer  function 
from  the  driving  point  of  the  double  capacitor  to  each  of  the  operational  amplifier 
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output  voltages,  we  obtain 


(3.14) 


R.2  +  R3 

Gi  =  ~rT’ 

G2=  1  +  (3.15) 

K3 

where  G\  and  G2  are  the  transfer  characteristics  from  the  voltage  vn  to  the  outputs 
of  the  top  and  bottom  amplifiers,  respectively.  In  order  to  select  a  valid  set  of  resistor 
values,  the  range  of  voltages  seen  at  the  top  of  the  double  capacitor  is  needed.  For  a 
single  soliton  solution,  the  closed  form  solution  for  the  voltage  is 


vn(t)  =  vt  In  jcosh  (p{n)  —  (3t\Ja/v^j  j  —  vt  In  jcosh  (p(n  +  1)  —  fit\Ja/v^J  j  +  const. 

(3.16) 

The  limiting  voltage  in  Eq.  (3.16)  is  given  by 


lim  vn(t) 

t->  00  v  ' 


vtp  +  constant, 


and 


lim  vn(t)  =  —vtp  +  constant. 

t-»— 00 


Selecting  the  constants  such  that  vn(— 00)  =  0,  gives 


(3.17) 


(3.18) 


fim  vn(t)  =  2 vtp.  (3.19) 

For  p  «  14,  this  leads  to  a  final  voltage  amplitude  on  the  order  of  vn  «  0.75  volts.  For 
each  soliton  that  passes  through  a  given  node,  the  voltage  on  the  double  capacitor 
will  increase  by  2 vtp.  In  order  to  keep  the  amplifiers  from  saturating  due  to  a  single 
soliton,  Gi  and  G2  must  remain  less  than  about  a  factor  of  10.  Since  each  soliton  that 
passes  through  a  given  node  will  result  in  a  similar  voltage  increase,  we  would  like 
these  gains  to  be  as  small  as  possible  to  avoid  saturation.  However,  it  is  also  important 
to  maintain  voltages  levels  in  the  double  capacitors  that  are  above  the  noise  level  for 
the  circuit,  which  implies  that  these  gains  cannot  be  too  small.  A  reasonable  balance 
can  be  obtained  by  setting  i?2  ~  R3,  and  R\  <C  l/C  which  can  be  met  by  selecting 
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Figure  3-5:  HSPICE  simulation  of  the  evolution  of  a  two-soliton  signal  through  the 
diode  lattice.  Each  horizontal  trace  shows  the  current  through  one  of  the  diodes  1, 3, 4 
and  5. 


Ri  =  R.2  =  R$  =  Ikfi  and  C  —  .OlyxF.  These  values  permit  soliton  pulse  widths 
of  about  5/iS  with  amplitudes  of  about  10mA  and  with  voltages  at  the  amplifier 
outputs  within  the  double  capacitors  on  the  order  of  1  Volt.  The  double  capacitors 
use  precision  LT1028A  operational  amplifiers  with  a  gain  bandwidth  product  of  about 
65  MHz.  Shown  in  Fig.  3-5  is  an  HSPICE  simulation  with  two  solitons  propagating 
down  a  length  10  Toda  chain. 

A  significant  difference  between  soliton  solutions  to  this  circuit  and  those  of  the 
nonlinear  LC  line  lies  in  the  scale  of  operation.  Due  to  biasing  constraints  for  the 
LC  line,  solitons  were  generally  restricted  to  a  small  range  of  wavenumbers  in  the 
neighborhood  of  p  «  1.  Over  this  range,  the  propagation  velocity  of  the  solitons, 
which  is  proportional  to  sinh(p)/p  does  not  vary  greatly  between  solitons  of  different 
wavenumbers.  This  led  to  the  use  of  chains  with  hundreds  of  nodes  in  order  to  demon¬ 
strate  soliton  collisions.  The  diode  ladder  circuit,  however,  can  operate  in  the  range 
p  «  14  for  solitons  with  amplitudes  in  the  mA  range.  Due  to  the  exponential  nature 
of  the  sinh  function,  the  velocities  of  solitons  with  slightly  different  amplitudes  for 
currents  in  the  mA  range  yield  exponentially  different  velocities.  This  enables  soliton 
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collisions  to  take  place  with  far  fewer  nodes  than  with  the  nonlinear  LC  network. 

As  illustrated  in  the  bottom  trace  of  Fig.  3-5,  a  soliton  can  be  generated  by 
driving  the  circuit  with  a  square  pulse  of  approximately  the  same  area  as  the  desired 
soliton.  As  seen  on  the  third  node  in  the  lattice,  once  the  soliton  is  excited,  the  non- 
soliton  components  are  quickly  stripped  away.  For  the  example  in  the  figure,  a  small 
pulse  followed  by  a  larger  pulse  are  used  to  drive  the  circuit  giving  rise  to  a  small 
soliton  followed  by  a  larger  amplitude  soliton.  This  property  has  been  demonstrated 
experimentally  for  a  number  of  soliton  systems,  c.f.  [23]  for  the  nonlinear  Schrodinger 
equation  and  [29]  for  the  Toda  lattice.  It  has  been  shown  theoretically  for  KdV,  c.f. 
[1]  and  [12],  that  practically  any  smooth,  localized  disturbance  of  the  proper  area  will 
result  in  a  soliton  with  that  area,  if  such  a  solution  exists.  Multi-soliton  signals  could 
also  be  generated  by  using  inverse  scattering  techniques  to  determine  a  drive  signal 
iin(t)  that  would  give  rise  to  the  desired  solitons.  However,  this  method  requires  more 
complex  drive  circuitry  than  simple  pulse  generators. 

Note  that  as  the  faster  soliton  overtakes  the  slower  as  viewed  on  the  fourth  node 
in  Fig.  3-5,  the  joint  signal  amplitude  is  significantly  less  than  the  sum  of  the  indi¬ 
vidual  amplitudes.  Also,  the  signal  shape  changes  significantly  during  the  nonlinear 
interaction.  These  two  effects  will  impact  both  the  energy  of  multi-soliton  signals  and 
the  ability  to  recover  their  signal  parameters. 

3.4  Circuit  Implementation 

To  perform  real-time  experimentation  and  to  verify  the  operation  of  the  model  using 
standard  components,  the  diode  ladder  circuit  has  been  implemented  in  hardware. 
Real-time  implementation  also  enables  rapid  testing  of  soliton  processing  techniques 
and  enables  measurements  of  actual  circuit  noise  levels.  Such  noise  measurements 
permit  experimental  verification  of  some  of  the  theoretical  results  we  obtain  in  Chap.  5 
concerning  system  noise. 

In  the  construction  of  the  circuit,  there  were  several  practical  matters  to  be  dealt 
with.  First,  the  diode  ladder  is  driven  by  a  current  source.  In  our  implementation 
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Figure  3-6:  Precision  bipolar  current  source. 

the  precision  bipolarity  current  source  shown  in  Fig.  3-6  taken  from  [30]  was  used. 
When  implemented  with  the  low  noise  LT1028A  operational  amplifiers,  this  circuit 
provides  a  reliable,  accurate  current  source,  with  low  leakage.  In  practice,  leakage 
current  turns  out  to  be  a  problem,  since  the  double  capacitor  circuits  are  marginally 
stable.  The  node  voltages  are  double  integrators  of  their  current  and  therefore  any 
excess  current  will  lead  to  large  deviations  in  the  node  voltages  and  corrupt  soliton 
propagation. 

In  addition  to  the  voltage  deviations  from  leakage  current,  each  soliton  that  passes 
through  a  node  on  the  ladder  contributes  a  net  voltage  increase  of  2 vtp  or  approxi¬ 
mately  1  Volt.  Therefore  a  signal  containing  three  solitons  will  leave  the  node  with 
a  net  voltage  increase  of  nearly  3  volts.  If  several  such  signals  are  processed  by  this 
circuit,  the  operational  amplifiers  in  the  double  capacitors  will  eventually  saturate. 
This  problem  can  be  overcome  by  resetting  the  node  voltages  after  each  signal  has 
been  processed  by  the  circuit  using  analog  switches  as  shown  in  Fig.  3-7. 

Finally,  since  the  solitons  are  present  in  the  diode  ladder  circuit  as  current  wave¬ 
forms,  there  must  be  an  adequate  means  of  measuring  the  current  through  the  diodes 
without  significantly  affecting  the  dynamics  of  the  circuit.  This  can  be  accomplished 
by  placing  a  small  resistance  in  series  with  each  diode  in  the  lattice  a  shown  in  Fig.  3- 
7.  The  current  through  the  diodes  can  then  be  observed  by  measuring  the  voltage 
drop  across  each  of  the  resistors  with  a  differential  amplifier. 

A  hardware  implementation  of  the  diode  ladder  circuit  with  twelve  nodes  is  shown 
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Figure  3-7:  Diagram  for  the  double  capacitors  used  in  the  diode  ladder  circuit.  Analog 
switches,  placed  in  parallel  with  the  capacitors,  are  used  to  reset  the  circuit  after  each 
processed  signal. 

in  Fig.  3-8.  Each  of  the  three  rightmost  bread  boards  in  the  figure  contains  four  stages 
of  diodes,  series  resistors,  and  double  capacitors.  The  left  most  bread  board  con¬ 
tains  pulse-generation  circuitry,  and  the  remaining  bread  board  contains  the  voltage- 
controlled  current  source.  A  two-soliton  signal  generated  by  this  circuit  is  shown  on 
the  oscilloscope  traces  in  Fig.  3-9.  The  bottom  trace  in  the  figure  corresponds  to  the 
input  current  to  the  circuit,  and  the  remaining  traces,  from  bottom  to  top,  show  the 
current  through  the  second,  third  and  fourth  diodes  in  the  lattice. 

For  a  more  complex  example,  a  simple  waveform  consisting  of  three  component 
soliton  signals,  periodically  repeated,  was  used  to  drive  the  diode  ladder  circuit  with 
several  additional  stages.  Using  a  digital  oscilloscope  to  sample  real-time  circuit 
waveforms,  measurements  of  the  diode  currents  were  transferred  to  a  computer,  and 
then  plotted  online.  The  measured  currents  through  each  of  the  diodes  are  shown  in 
Fig.  3-10.  The  time  axis  of  the  figure  is  such  that  t  —  0  corresponds  to  the  beginning  of 
a  period.  The  largest  amplitude  soliton  measures  17mA,  with  a  pulse  width  of  82/us  at 
the  first  node.  As  measured  on  the  fifteenth  node,  the  amplitude  is  14mA  with  pulse 
width  of  86/us.  This  decay  in  the  soliton  amplitude  is  on  the  order  of  1%  per  node 
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Figure  3-8:  Hardware  implementation  of  the  diode  ladder  circuit.  The  first  column 
from  the  left  contains  the  pulse-generation  circuitry;  the  second  contains  the  voltage- 
controlled  current  source;  each  of  the  last  three  contains  four  stages  of  diodes,  series 
resistors  and  double  capacitors. 
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Figure  3-9:  Oscilloscope  traces  for  two  solitons  in  the  diode  ladder  circuit.  The  traces 
correspond  to  the  currents  in  the  first  four  diodes. 

and  may  have  several  causes  in  addition  to  deviations  of  the  circuit  components  from 
their  idealized  models.  Specifically,  as  stated  in  [32],  dissipative  effects  in  the  lattice 
are  contrary  to  the  conservative  nature  of  the  Toda  lattice,  and  will  necessarily  lead 
to  energy  loss.  Also,  as  shown  in  [50],  inter-node  parameter  fluctuations  can  lead  to 
dispersion,  causing  decay  as  well  as  additional  non-soliton  components.  This  leads  to 
a  change  in  the  soli  ton  parameter  /?,  resulting  in  a  decrease  in  soliton  velocity  as  they 
propagate  through  the  lattice.  In  the  figure,  there  is  also  a  small  spike  that  appears 
in  each  of  the  diode  currents  near  the  time  t  =  —1ms.  This  results  from  the  reset 
signal,  Vsy/ ,  propagating  down  the  lattice  and  resetting  adjacent  double  capacitors  at 
slightly  different  times. 


3.5  Circuit  Model  for  Discrete-KdV 

The  discrete-KdV  equation  (dKdV),  sometimes  referred  to  as  the  nonlinear  ladder 
equations  [1],  or  the  KM  system  (Kac  and  vanMoerbeke)  [67]  is  governed  by  the 
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Figure  3-10:  Diode  currents  measured  from  the  diode  ladder  circuit  in  operation.  The 
input  signal  consists  of  three  square  pulses  of  different  areas.  The  spike  that  appears 
in  the  figure  near  t  =  -1  ms  is  a  result  of  the  signal  that  resets  the  lattice. 
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equations, 


(3.20) 


un(t)  =  eUn~1^  -  e“n+1^. 

This  equation  is  first  order  in  time,  which  makes  the  dynamics  less  complex  than  the 
Toda  lattice,  yet  as  evidenced  in  Sec.  2.2,  the  inverse  scattering  framework  is  very 
similar.  Although  much  of  the  theory  for  this  nonlinear  system  has  been  developed 
with  the  theory  for  the  Toda  lattice,  the  discrete-KdV  equations  are  generally  ignored 
since  there  is  no  clear  physical  analog  of  these  equations.  However,  there  is  special  re¬ 
lationship  known  as  a  Backlund  transformation  which  provides  a  connection  between 
this  system  and  the  Toda  lattice  [67,  33,  34]. 

Given  our  success  with  the  Toda  lattice  using  a  ladder  comprising  diodes  and 
double  capacitors  and  the  similarity  between  (3.20)  and  (3.7),  we  first  consider  a 
ladder  of  diodes  with  shunt  capacitors.  A  similar  analysis  leads  to  the  following  set 
of  equations 

dvn{t)  _  h  ^e(u„_i (t)-Vn(t))/vt  _  e(vn(t)-vn+1(t))/vt'j  ,  (3.21) 

which  look  deceptively  similar  to  (3.20).  However,  there  is  no  apparent  means  of 
decoupling  the  node  voltage,  vn(t)  from  vn+\(t)  and  u„_i (t)  as  would  be  required. 
This  similarity  is  not  a  mere  coincidence  and  actually  leads  to  a  realization  of  the 
discrete-KdV  equation  using  two  Toda  circuits.  Following  [33]  and  [67],  let 

un  — y  —Rn,  t  —y-  — t,  (3.22) 


which  transforms  (3.20)  to 


Rn  = 


_  g  Rn+1 


(3.23) 


Letting  qn  =  Rn  +  Rn+i,  then  qn  satisfies 

qn  =  2e~9n  -  e~9n~2  -  egn+2.  (3.24) 
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Figure  3-11:  Illustration  of  the  relationship  between  two  adjacent  Toda  lattices, 
=  1,2  and  the  discrete-KdV  equation.  This  process  suggests  a  possible  im¬ 
plementation  of  the  dKdV  equation  using  two  adjacent  Toda  lattice  circuits. 

Taking  every  other  term,  i.e.  =  q2n,  <?2,n  =  <?2n+i,  we  have 

=  2e~Qi'n  -  e~gi’n~1  -  eQi-n+1 ,  *  =  1,2.  (3.25) 

Setting  qiin  =  -  ln(l  +  /i>n),  yields 

1  ~  2fi,n  +  fi,n+i  *  =  1,2,  (3.26) 

which  are  each  a  Toda  lattice  equation.  The  physical  interpretation  of  this  transfor¬ 
mation  is  the  following:  if  /i)n  and  /2,n  are  each  defined  within  a  different  Toda  lattice, 
then  Rn  defined  by  R2n  =  fi,n  ~  /2,n,  ^Wi  =  /2,n+i  -  fi,n  satisfies  the  discrete-KdV 
equation  in  the  form  (3.23).  This  process  is  illustrated  in  Fig.  3-11.  Although  a  dKdV 
circuit  could  be  so  constructed,  the  resulting  circuitry  would  be  twice  as  complex  as 
the  Toda  circuit.  A  much  simpler  implementation  can  be  found  by  maintaining  the 
aspects  of  the  diode  ladder  that  are  useful,  namely  the  exponential  current  relation¬ 
ship  of  the  diodes,  while  removing  the  aspects  which  are  troublesome,  viz.  the  ladder 
interconnections. 

Since  the  desired  equations  are  first  order,  we  use  capacitor  voltages  for  state 
variables,  i.e.  vn(t )  will  be  the  voltage  on  the  n-th  capacitor.  Rather  than  assembling 
the  capacitors  in  a  ladder  network,  we  consider  a  collection  of  nodes  with  nearest 
neighbor  coupling  as  shown  in  Fig.  3-12.  Each  node  is  maintains  a  node  voltage,  vn(t), 
and  also  maintains  a  voltage  that  is  proportional  to  eVn^\  which  can  be  accomplished 
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Figure  3-12:  Collection  of  nodes  for  the  discrete-KdV  circuit. 


with  a  voltage  follower  and  a  diode  as  shown  in  Fig.  3-12.  Since  the  voltage  follower 
mirrors  the  voltage  on  the  capacitor,  neglecting  the  voltage  drop  from  the  resistor, 
the  voltage  across  the  diode  is  approximately  vn(t).  Hence  the  current  through  the 
diode  is  in{t)  ~  Is(exp(vn(t)  /  vt)  —  1),  where  vt  is  the  thermal  voltage.  Since  each 
node  maintains  a  node  voltage  along  with  a  voltage  proportional  to  exp(n„(t)),  all 
that  remains  is  to  construct  a  current  source  that  is  proportional  to  the  difference 
in  the  exponential  reference  voltages  of  the  neighboring  nodes.  If  a  current  source 
is  used  to  drive  the  capacitor  as  shown  in  Fig.  3-12,  then  the  node  voltage,  vn(t),  is 
governed  by 

vn(t)  =  ^  _  evn+1(t)M)  .  (3.27) 

The  required  differential  voltage  controlled  current  source  can  be  the  same  as  the  one 
that  drives  the  Toda  ladder  circuit,  shown  in  Fig.  3-6.  Therefore,  for  a  collection  of 
nodes  such  as  those  in  Fig.  3-12,  the  node  capacitor  voltages  are  governed  by  the 
discrete-KdV  equation. 

The  time  scale  of  the  circuit  can  be  set  by  proper  choice  of  the  ratio  Is/C.  Specif¬ 
ically,  if  Is  =  avtC,  the  node  voltage  satisfies, 


dvn(r)/vt 

dr 


—  ^evn-i(r)/vt  _  evn+i(T)/vt'j  ^ 


(3.28) 


where  r  =  t/a.  Thus  vn(t)/vt  satisfies  the  discrete-KdV  equation  on  a  time-scale  t/a. 


56 


An  HSPICE  simulation  of  this  circuit  indeed  verifies  the  propagation  of  dKdV 
solitons.  Since  this  circuit  is  first  order,  the  state  of  the  system  is  completely  specified 
by  the  capacitor  voltages.  Rather  than  processing  continuous-time  signals  as  with 
the  Toda  lattice  system,  we  can  use  this  system  to  process  discrete-time  solitons 
as  specified  by  vn.  For  the  purposes  of  simulation,  we  consider  the  periodic  dKdV 
equation  by  setting  ura+i(t)  =  Vo(t)  and  initializing  the  system  with  the  discrete-time 
signal  corresponding  to  a  listing  of  node  capacitor  voltages.  We  can  place  a  multi- 
soliton  solution  in  the  circuit  using  inverse  scattering  techniques  to  construct  the 
initial  voltage  profile.  The  single  soliton  solution  to  the  dKdV  system  is  given  by 

m  ,  ( cosh (7  (n-2)-pt)  cosh (7  (n+l)-pt)\ 

Vn  n  \  cosh (7  (n  —  1)  —  f3t)  cosh(7  n  —  fl  t)  ) 

where  /3  =  sinh(27).  Shown  in  Fig.  3-13,  is  the  result  of  an  HSPICE  simulation 
of  the  circuit  with  30  nodes  in  a  loop  configuration,  each  with  2 nF  capacitors  and 
diodes  with  a  saturation  current  chosen  to  be  Is  =  vt  x  2nA.  Thus,  the  time  scale 
of  the  circuit  is  unity,  a  =  1.  The  initial  condition  was  set  such  that  a  soliton  with 
7x  =  2.5  was  placed  on  node  0  and  a  second  soliton  with  72  =  2  was  placed  on  node 
10.  As  with  the  diode  ladder  implementation  of  the  Toda  circuit,  this  circuit  model  is 
very  accurate,  and  can  support  a  wide  range  of  soliton  wavenumbers.  Discrete-time 
soliton  signals  could  be  used  as  an  initial  condition  to  the  dKdV  lattice  and  then  be 
processed  by  their  evolution  in  the  circuit. 

3.6  Further  Considerations 

In  this  chapter  we  have  illustrated  three  circuits  which  can  be  used  to  both  generate 
and  to  process  soliton  signals.  The  first  of  the  three,  developed  by  Hirota  and  Suzuki, 
demonstrated  the  feasibility  of  using  such  nonlinear  circuits  as  analogs  of  soliton 
systems.  Their  results  [29,  64,  65],  along  with  the  experimental  work  of  Toda  in  [68], 
Scott  in  [56,  57],  and  a  host  of  others  illustrate  the  potential  for  development  of  a 
variety  of  analog  hardware  for  soliton  processing  techniques. 
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Figure  3-13:  To  the  left,  the  normalized  node  capacitor  voltages,  vn(t)/vt  for  each 
node  is  shown  as  a  function  of  time.  To  the  right,  the  state  of  the  circuit  is  shown 
as  a  function  of  node  index  for  five  different  sample  times.  The  bottom  trace  in  the 
figure  corresponds  to  the  initial  condition. 

We  then  presented  a  new  diode  ladder  circuit  which  models  the  Toda  lattice  more 
accurately  than  the  nonlinear  LC  circuit  and  can  be  operated  over  a  larger  range 
of  soliton  amplitudes.  Increasing  the  range  of  operation  has  a  direct  impact  on  the 
complexity  of  the  processing  hardware  that  would  be  used,  for  example,  to  separate 
multi-soliton  signals.  Specifically,  the  number  of  nodes  required  (length  of  the  lattice) 
is  significantly  reduced.  A  real-time  implementation  of  this  circuit  was  constructed 
in  hardware  using  standard  components  which  appears  to  be  the  first  such  circuit  to 
experimentally  demonstrate  true  soliton  collisions.  It  remains  to  be  seen  whether  or 
not  the  diode  ladder  circuit  can  operate  without  the  need  for  periodic  resetting  of 
the  double  capacitor  sub-circuits.  The  resolution  of  this  issue  might  avail  the  diode 
ladder  circuit  to  a  greater  range  of  potential  processing  techniques  and  is  therefore 
an  interesting  avenue  for  future  work. 

Finally,  we  presented  a  new  circuit  which  implements  the  discrete-KdV  equation, 
for  which  no  other  circuit  models  have  previously  been  presented.  Experiments  using 
HSPICE  demonstrate  the  viability  of  the  circuit  model  for  a  variety  of  discrete-time 
soliton  processing  techniques.  As  shown,  the  dKdV  circuit  is  placed  in  a  loop  con¬ 
figuration  and  the  circuit  iterates  forward  from  an  initial  condition.  This  circuit  can 
also  process  continuous  signals  in  an  analogous  manner  to  the  Toda  lattice.  However, 
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the  problem  of  properly  terminating  the  dKdV  circuit  differs  from  its  Toda  lattice 
counterpart  due  to  the  more  complex  circuit  topology. 

The  circuits  developed  in  this  chapter  can  be  used  to  generate  multi-soliton  signals 
as  well  as  perform  a  variety  of  processing  operations.  For  example  multiple  solitons 
can  be  multiplexed  by  using  a  signal  with  solitons  arranged  in  increasing  amplitude 
as  input,  allowing  them  to  collide,  and  then  extracting  the  signal  from  the  circuit. 
Similarly  the  separation  of  multiple  overlapping  solitons  could  be  achieved  by  allowing 
them  to  propagate  at  different  velocities  and  again  extracting  the  signal  after  sepa¬ 
ration.  Each  of  these  otherwise  complex  nonlinear  tasks  can  be  completed  naturally 
by  the  dynamics  of  the  nonlinear  systems.  Such  systems  also  indicate  the  viability 
of  analog  hardware  implementations  of  a  large  class  of  nonlinear  systems  exhibiting 
soliton  behavior.  Additionally,  many  of  the  soliton  processing  techniques  could  be 
mapped  onto  the  soliton  or  soliton-like  behavior  that  has  been  exhibited  by  a  large 
class  of  cellular  automata  [1,  7,  53,  62,  63,  66],  for  which  many  fast  algorithms  [52] 
and  hardware  [46,  70,  71]  have  been  developed. 
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Chapter  4 


Communication  with  Soliton 
Signals 


Many  traditional  communication  systems  use  a  form  of  sinusoidal  carrier  modula¬ 
tion,  such  as  amplitude  modulation  (AM)  or  frequency  modulation  (FM)  to  transmit 
a  message-bearing  signal  over  a  physical  channel.  The  reliance  upon  sinusoidal  sig¬ 
nals  is  due  in  part  to  the  simplicity  with  which  such  signals  can  be  generated  and 
processed  using  linear  systems.  More  importantly,  information  contained  in  sinu¬ 
soidal  signals  with  different  frequencies  can  easily  be  separated  using  linear  systems 
or  Fourier  techniques.  This  has  led  to  the  development  of  standards  and  regulation 
of  the  sinusoidal  composition  of  signals  transmitted  electromagnetically  based  on  the 
allocation  of  separate  bands  of  frequencies  to  different  users. 

The  complex  dynamic  structure  of  soliton  signals  and  the  ease  with  which  these 
signals  can  be  both  generated  and  processed  with  analog  circuitry  renders  them  po¬ 
tentially  applicable  in  the  broad  context  of  communication  in  an  analogous  manner 
to  sinusoidal  signals.  We  have  also  seen  that  either  by  using  analog  circuitry  or 
the  inverse  scattering  transform,  the  individual  solitons  in  a  multi-soliton  signal  can 
also  be  easily  separated.  In  this  chapter  we  will  present  a  paradigm  for  modulation 
of  information  on  soliton  carriers  and  illustrate  the  potential  applicability  of  such 
a  framework  to  multi-user  communication.  Specifically,  we  consider  a  scenario  in 
which  solitons  are  used  as  carrier  waveforms  in  an  amplitude-  or  position-modulation 
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Figure  4-1:  A  soliton  carrier  signal  for  the  Toda  lattice. 


communication  system  using  the  soliton  systems  to  perform  signal  multiplexing  and 
demultiplexing.  In  this  context,  we  will  show  that  the  signal  amplitude  reduction 
that  occurs  during  soliton  interaction  may  provide  increased  energy  efficiency,  mak¬ 
ing  such  techniques  particularly  attractive  for  a  broad  range  communication  contexts 
over  power-limited  channels  including  portable  wireless  communication,  satellite  or 
deep-space  applications. 

4.1  Soliton  Modulation 

Although  the  mechanisms  for  the  generation  and  processing  of  soliton  signals  differs, 
many  standard  sinusoidal  carrier  modulation  techniques  could  be  implemented  using 
soliton  carriers.  In  this  chapter,  we  define  a  soliton  carrier  as  a  signal  which  is 
composed  of  a  periodically-repeated  single  soliton  solution  to  a  particular  nonlinear 
system.  For  example,  a  soliton  carrier  signal  for  the  Toda  lattice  is  shown  in  Fig.  4-1 
and  can  be  written, 

U{t)  =  ^  |  Y  /?2sech2  [a(n  -  A i)  -  /3t]  -  2 /?*/ j  ,  (4.1) 

where  a  and  f3  are  parameters  related  to  the  elliptic  functions  defined  in  Chap  2, 
Eq.  (2.11).  In  addition  to  the  Toda  lattice,  periodic  soliton  solutions  of  many  soliton 
systems  can  be  written  in  terms  of  Jacobian  elliptic  functions  as  in  Eq.  (2.10). 

As  a  soliton  carrier  such  as  that  for  the  Toda  lattice  in  Fig.  4-1  is  generated, 
a  simple  amplitude  modulation  scheme  could  be  devised  by  modulating  the  soliton 
parameter  /?,  since  the  amplitude  of  the  Toda  lattice  solitons  is  proportional  to  j32. 
Note  that  the  pulse-width  also  changes  with  /?,  giving  rise  to  a  scale  modulation 


Figure  4-2:  Modulating  the  relative  amplitude  or  position  of  soliton  carrier  signal  for 
the  Toda  lattice. 


rather  than  a  pure  amplitude  modulation.  Similarly,  an  analog  of  FM /PM  or  pulse- 
position  modulation  could  be  achieved  by  modulating  the  relative  position  of  each 
soliton  in  a  given  period,  as  shown  in  Fig.  4-2.  Although  the  mechanisms  developed 
for  sinusoidal  carrier  demodulation  would  also  have  to  be  adapted,  at  least  from  a 
signal  generation  viewpoint,  a  variety  of  sinusoidal  carrier  modulation  techniques  can 
be  directly  mapped  onto  soliton  carriers.  We  have  again  focused  on  the  Toda  lattice 
as  an  example,  however,  parameter  modulation  of  soliton  carriers  could  be  applied  to 
the  large  class  of  soliton  signals  that  can  be  both  generated  and  processed  by  solvable 
nonlinear  systems. 

4.2  Soliton  Multiplexing 

As  a  simple  extension,  the  soliton  modulation  techniques  described  in  Section  4.1  can 
be  generalized  to  include  multiple  solitons  in  each  period  and  accommodate  multiple 
information-bearing  signals,  as  shown  in  Fig.  4-3  for  a  four  soliton  example  using  the 
Toda  lattice  circuits  developed  in  Chap.  3.  In  the  figure,  a  signal  is  generated  as 
a  periodically-repeated  train  of  four  solitons  of  increasing  amplitude.  The  relative 
amplitudes  or  positions  of  each  of  the  component  solitons  could  be  independently 
modulated  about  their  nominal  values  to  accommodate  multiple  information  signals 
in  a  single  soliton  carrier. 

The  nominal  soliton  amplitudes  can  be  appropriately  chosen  so  that  as  this  signal 
is  processed  by  the  diode  ladder  circuit,  the  larger  amplitude  solitons  propagate  faster 
than  the  smaller  solitons,  and  each  of  the  solitons  can  become  nonlinearly  superim- 
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Figure  4-3:  Multiplexing  of  a  four  soliton  solution  to  the  Toda  lattice. 


posed  as  viewed  at  a  given  node  in  the  circuit.  From  an  input-output  perspective,  the 
diode  ladder  circuit  can  be  used  to  make  each  of  the  solitons  coincidental  in  time.  As 
indicated  in  the  figure,  this  packetized  soliton  carrier  could  then  be  transmitted  over 
a  communication  channel.  At  the  receiver,  the  multi-soliton  signal  can  be  processed 
with  an  identical  diode  ladder  circuit  which  is  naturally-suited  to  perform  the  non¬ 
linear  signal  separation  required  to  demultiplex  the  multiple  soliton  carriers.  As  the 
larger  amplitude  solitons  emerge  before  the  smaller,  after  a  given  number  of  nodes, 
the  original  multi-soliton  carrier  re-emerges  from  the  receiver  in  amplitude-reversed 
order.  At  this  point,  each  of  the  component  soliton  carriers  could  be  demodulated  to 
recover  the  individual  message  signals  it  contains.  Aside  from  a  packetization  of  the 
component  solitons,  we  will  see  that  multiplexing  the  soliton  carriers  in  this  fashion 
can  lead  to  an  increased  energy  efficiency  for  such  carrier  modulation  schemes. 

Since  the  Toda  lattice  equations  are  symmetric  with  respect  to  time  and  node 
index,  solitons  can  propagate  in  either  direction.  As  a  result,  a  single  diode  ladder 
implementation  could  be  used  as  both  a  modulator  and  demodulator  simultaneously, 
as  shown  in  Fig.  4-4.  From  the  figure  at  time  t  =  0,  a  signal  composed  of  three  solitons 
is  multiplexed  by  the  system,  by  propagating  the  separated  solitons  along  the  lattice 
from  node  0  to  node  35.  The  resulting  multiplexed  signal  can  be  extracted  from  node 
35  around  t  =  30.  Also  at  t  —  0,  a  received  signal  composed  of  three  multiplexed 
solitons  is  input  to  the  system  at  node  35.  These  solitons  are  induced  in  the  reverse 
direction,  allowing  them  to  separate  by  the  time  they  arrive  at  node  0  around  t  = 
30.  Since  the  forward  propagating  solitons  correspond  to  positive  eigenvalues  in 
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Figure  4-4:  Toda  lattice  response  to  an  outgoing  signal  at  node  0  and  to  a  received 
signal  at  node  35,  each  with  three  component  solitons.  The  solitons  that  are  input  at 
node  0  are  multiplexed  by  the  system,  as  viewed  on  node  35.  The  solitons  that  are 
input  to  node  35,  propagate  in  the  reverse  direction,  and  are  demultiplexed  by  the 
system,  as  viewed  on  node  0. 

the  inverse  scattering  transform  and  the  reverse  propagating  solitons  have  negative 
eigenvalues,  the  dynamics  of  the  two  signals  will  be  completely  decoupled.  The  ability 
to  propagate  solitons  in  both  directions  is  a  property  shared  by  many  soliton  systems 
with  dynamics  governed  by  differential  equations  that  are  second  order  in  time. 

Although  presented  in  the  context  of  the  Toda  lattice  equation,  the  soliton  mod¬ 
ulation  techniques  presented  are  applicable  to  a  variety  of  soliton  systems,  both  in 
discrete-  and  continuous-time.  For  example,  a  similar  technique  can  be  used  for  mod¬ 
ulation  of  the  discrete-time  solitons  generated  by  the  discrete-KdV  equation.  Both 
the  soliton  amplitude  and  relative  position  for  a  multi-soliton  solution  of  the  dKdV 
equation  could  be  independently  modulated  when  the  solitons  are  well-separated  in 
time.  This  multi-soliton  signal  could  then  be  used  as  the  initial  condition  to  a  dKdV 
circuit  developed  in  Chap.  3,  and  then  iterated  forward  in  time  until  the  component 
solitons  were  mutually  overlapping.  This  multiplexed  signal  could  be  transmitted 
to  a  receiver  where  identical  dKdV  circuitry  is  used  to  demultiplex  the  component 
solitons. 
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4.3  Spectrum  of  Toda  Lattice  Solitons 


To  draw  comparisons  with  sinusoidal  AM  and  FM  modulation,  we  consider  the  Fourier 
spectrum  of  a  periodic  single  soliton  carrier  for  the  Toda  lattice.  We  consider  a 
solution  of  the  form  of  Eq.  (4.1)  for  a  fixed  node, 

OO 

m  =  E  Aech 2mt  -  er))  -  w,  (4.2) 

£=—00 

where  T  is  the  period  of  one  message  frame.  The  Fourier  transform  of  a  single  soliton 
(or  one  frame)  is  given  by[13] 


/OO 

v(t)e~jutdt 

-OO 


mv 


sinh  (7ru/2P) 


27t/36(u). 


(4.3) 


In  the  small  amplitude  limit,  the  periodic  traveling  wave  solutions  to  the  Toda  lattice 
approach  sinusoids.  To  see  that  this  is  reasonable,  we  consider  the  Fourier  transform 
of  the  periodic  soliton,  Eq.  (4.2),  which,  ignoring  effects  due  to  time  aliasing,  is  given 
by 


7 xu 


F ^  t=S&o  sinh  (™/2p) 


6  ( u)  —  £u>c) 


(4.4) 


where  cuc  =  2-ir/T.  Since  the  Fourier  series  coefficients  fall  off  exponentially  in  fre¬ 
quency,  for  small  amplitudes  (0  <C  ouc), 


F{oj )  «  —  uc )  +  4tt0S(u  +  ue), 


(4.5) 


which  is  in  line  with  the  asymptotic  sinusoidal  form  of  f(t). 

We  now  consider  an  AM  modulation  with  a  single  soliton  waveform, 

fAM(t)  =  /(*)(  1  +  m(t)),  (4.6) 

where  m(t)  is  the  message  signal.  We  also  assume  that  the  modulation  depth  is 
small,  i.e.  the  bandwidth  of  m{t)  is  much  smaller  than  u >c.  To  illustrate  the  resulting 
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spectrum,  we  set  m(t)  =  mcos(u;mt),  which  yields 


Fam{u)  -  ^2 


FLU 


^  frr  ,/9/T)  ^  -  tUc)  +  -  fac  +  Um)  +  \s{u  -  £uc- 

e=-£&0 smh  K™/2(3)  2  2 


(4.7) 

Again,  making  the  small  amplitude  approximation  (/3  <C  luc),  the  positive  half  of  the 
spectrum  is  given  by 


F+(uj)  w  47r/36(u;  -  uic)  +  2f/35(lu  -  (luc  -  um))  +  2irf36(co  -  (u?c  -f  tum )),  (4.8) 

which  is  the  same  as  the  spectrum  for  a  sinusoidal  AM  modulation. 

In  analogy  to  FM,  we  consider  a  soliton  carrier  with  the  relative  position  of  the 
component  solitons  modulated, 

fFM(t)  =  f(t  +  77lCOs(wmt)),  (4.9) 


again  assuming  that  cum  <C  u>c.  This  expression  can  be  expanded  using  a  trigonometric 
identity,  to 


fFM(t)  =  ^2 


e=-oo  (cosh (f3ri)  cosh(/3m  cos(wmi))  +  sinh(/3Te)  sinh(/3m  cos(cumt)))y 


-2/3, 


(4.10) 

where  rt  =  t  —  IT .  Making  the  assumption  that  m  <C  T,  this  expression  can  be 
approximated  to  first  order  in  m  by 


fFM{t )  « 


(32 


cosh \&Ti)  +  2  cosh (Pti)  sinh(/3r^)m  cos(umt) 


-2/3 


y)  /32sech2(/3r£)  —  2/32sech2(/3 r^)  tanh(/3T£)racos(u;mi)  —  2/3j(4.11) 


&=— oo 


Noting  that 


— /32sech2(/3i)  =  —  2/33sech2(/3£)  tanh(/3t), 


(4.12) 
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we  obtain  the  approximate  Fourier  transform  of  /fat(^)  as 

p  (:A~  ™  i  Mv-um)2/ 2  +  um)2/2 

FM  sink  (ttu>/ 2 (3)  sinh  —  um)/2/3)  smh(n(u  +  um)/2f3)'  ' 

Again,  making  an  approximation  for  (3  <C  1,  and  noting  that  in  the  limit  of  small 
/?,  the  period  T  of  the  signal  must  approach  (3,  the  positive  half  of  the  spectrum  is 
approximated  by 

VFm,+(u)  ~  47t(36(uj  -  ue)  +j27rp5(u  -  (vc  -  u>m))  +  j2n  (35  (uu  -  (u>c  +  u>m)),  (4.14) 

which  has  the  same  magnitude  as  the  spectrum  for  a  sinusoidal  AM  modulation,  and 
is  the  same  spectrum  that  would  be  obtained  by  sinusoidal  FM  modulation. 

For  the  two-soliton  signal,  from  Chapter  2,  Eq.  (2.6),  a  single  period  of  the  signal 
is  of  the  form 

_  ffisech2(7?i)  +  /?fsech2(772)  +  Asech^T^sech2^) 

(cosh  (0/2)  +  sinh(0/2)  tanh(r/i)  tanh^))2 

with  rji  =  Pi(t  —  5i).  Thus,  for  positional  modulation,  Si  =  rrii{t),  in  addition  to 
the  spectral  components  of  the  individual  messages,  there  is  a  contribution  from 
the  product  of  the  individually  modulated  solitons.  However,  when  the  solitons  do 
not  overlap,  A  «  0  and  the  modulation  is  essentially  the  sum  of  the  individually 
modulated  waveforms.  As  the  solitons  begin  to  overlap,  the  contribution  from  the 
cross  terms  becomes  significant,  and  spectral  mixing  of  the  component  messages  will 
occur.  This  results  in  an  expansion  of  Fourier  bandwidth  of  the  multi-soliton  signal. 
Examining  the  Fourier  bandwidth  of  soliton  signals  is  of  interest  for  several  reasons. 
First,  such  analysis  offers  comparisons  to  the  spectral  analyses  of  standard  sinusoidal 
modulation  techniques.  Also,  our  engineering  intuition  often  relates  the  bandwidth 
of  a  signal  to  the  available  degrees  of  freedom  or  information  content  of  the  signal. 
Finally,  most  transmission  standards  for  electromagnetic  as  well  as  wired  propagation 
are  governed  by  an  allocation  of  a  finite  region  of  Fourier  spectrum  or  employ  channels 
which  are  frequency-selective,  making  such  analyses  critical  for  use  in  any  practical 


68 


transmission  context. 

However,  in  the  absence  of  regulation  on  the  Fourier  content  of  signals,  it  is  per¬ 
haps  more  appropriate  to  consider  the  nonlinear  spectrum,  as  defined  in  Chap.  2, 
when  analyzing  soliton  carrier  modulation  techniques.  The  nonlinear  spectrum  plays 
essentially  the  same  role  for  soliton  modulation  that  the  Fourier  spectrum  plays  for  si¬ 
nusoidal  modulation.  For  example,  in  sinusoidal  AM  or  FM,  multiple  carriers  occupy 
disjoint  regions  of  the  Fourier  spectrum,  enabling  signal  separation  via  LTI  process¬ 
ing.  Similarly,  in  the  soliton  carrier  modulation  techniques  outlined  in  this  chapter, 
multiple  carriers  would  occupy  disjoint  regions  of  the  nonlinear  spectrum,  or  different 
eigenvalues  of  the  linear  operator  L  in  the  inverse  scattering  transform.  This  enables 
separation  of  the  different  carrier  signals  through  evolution  of  the  nonlinear  equations 
or  through  an  inverse  scattering  analysis.  In  either  a  wired  context  or  another  situ¬ 
ation  where  all  of  the  signals  are  constructed  for  soliton  modulation,  although  there 
may  be  Fourier  spectral  broadening  of  the  component  signals,  there  is  no  expansion 
of  the  nonlinear  spectral  composition  of  multi-soliton  signals.  This  is  due  to  the  time- 
invariance  of  the  eigenvalues  in  the  inverse  scattering  transform  under  the  evolution 
of  the  nonlinear  dynamics. 

4.4  Low  Energy  Signaling 

In  Chapters  2  and  3,  we  observed  a  compaction  of  signal  amplitude  as  multiple  soli- 
tons  interact.  This  nonlinear  coupling  among  the  component  solitons  can  be  exploited 
to  yield  a  reduction  in  transmitted  signal  energy  for  soliton  carrier  modulation  tech¬ 
niques.  Since  the  nonlinear  coupling  of  multi-soliton  signals  is  system-specific,  we 
focus  our  attention  on  signal  amplitude  reduction  in  the  Toda  lattice.  Specifically,  it 
can  be  shown  [29]  that  for  a  solution  to  the  Toda  lattice  in  the  form  of  Eq.  (3.1), 

/OO 

vn(t)dt  =  constant,  (4.16) 

-OO 
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where  vn  is  the  voltage  on  the  n-th  nonlinear  capacitor, 


/oo 

in(t)dt  =  constant, 

-OO 


(4.17) 


where  in  is  the  current  through  the  n-th  inductor,  and  that 


/OO 

vn{t)in-i(t)dt  —  constant. 

-OO 


(4.18) 


As  a  consequence  of  these  conservation  laws,  we  see  that  if  two  solitons  are  well 
separated,  their  individual  contributions  in  Eqs.  (4.16)  and  (4.17)  are  added,  yet  as 
they  approach  one  another,  their  joint  amplitude  is  restricted  by  Eq.  (4.18).  For 
example,  if  two  solitons  are  well-separated  in  time  as  viewed  on  node  k,  and  mutually 
overlapping  as  viewed  on  node  A  we  obtain 

/oo  roo  roo 

{ve(t)ie-!{t))dt  =  /  vhk{ri)ii,k-i{ri)dt  +  /  v2,k{r2)i2,k-i(r2)dt 

-oo  J — oo  J— oo 

/oo 

(vi,*(r i)  +  v2^k(r2))(ii,k-i(ri)  +  i2,k-i{r2))dt,  (4.19) 

-OO 


where  vx^{ri)  is  the  voltage  contribution  from  the  first  soliton  when  the  solitons  are 
separated  as  viewed  on  node  k,  and  vt(t)  is  the  voltage  contribution  from  combined 
solitons  as  viewed  on  node  i.  Note  also  that  the  time  origin  of  each  signal  has  been 
shifted,  t  Ti  so  that  each  soliton  is  re-centered  about  t  =  0.  Due  to  the  conservation 
law  of  Eq.  (4.18),  the  joint  amplitude  of  the  solitons  when  combined  must  be  less  than 
the  sum  of  the  individual  amplitudes  when  separated  as  illustrated  in  Fig.  4-5.  A 
consequence  of  the  restriction  placed  on  the  amplitude  of  the  overlapping  solitons 
is  a  reduction  of  energy  in  the  transmitted  signal  for  the  modulation  techniques  of 
Section  4.2.  In  fact,  as  a  function  of  the  relative  separation  of  the  two  solitons,  the 
minimum  energy  of  the  transmitted  signal  is  obtained  precisely  at  the  point  of  overlap. 
This  can  be  shown  for  the  two-soliton  case  by  analysis  of  the  form  of  the  equation  for 
the  energy  in  the  voltage  waveform,  u(t), 


/OO 

v(t-,6u52)2dt, 

-OO 


(4.20) 
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Figure  4-5:  A  two-soliton  solution  is  depicted  in  the  Toda  lattice.  Each  horizontal 
trace  is  the  response  at  a  successive  node  in  the  lattice.  In  this  case,  the  two  soliton 
wavenumbers  are  p\  —  2  and  p2  =  1-3. 

where  v(t;  81,82)  is  given  in  Eq.  (4.15). 

Without  loss  of  generality,  we  assume  81  =  0,  and  seek  the  value  of  82  that 
minimizes  (4.20).  Differentiation  of  (4.20),  yields, 

dE_ 
d82 


/OO 

2  v(t) 

-OO 


dv(t) 

dS2 


dt. 


(4.21) 
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We  now  seek  the  value  of  82  that  makes  the  integral  (4.21)  equal  to  zero.  This  is 
accomplished  by  first  noting  that  when  the  two  solitons  have  the  same  phase,  i.e., 
<5x  —  S2  in  (4.15),  then  v(t)  is  an  even  function  of  time,  centered  about  5i  —  52.  In 
this  case,  setting  S2  =  <5i  =  0,  makes  v(t)  even.  If  it  can  be  shown  that  by  setting 
S2  =  0,  that  dv(t)/d52  is  an  odd  function,  then  the  integral  in  (4.21)  is  trivially  zero. 
This  indeed  turns  out  to  be  the  case, 


dv(t)  (2  /?isech2(/?i.  t)  +  2  /3fsech2(j72)  +  2  A  sech2(/?x  t)sech2(7?2)) 

°^2  (cosh(|)  +  sinh(|)  tanh(/3x  t )  tanh(772)) 

x  sinh(^)  tanh(/?x  t )  fl  —  tanh2^))  /?2 

Z 

[  2  /3|sech2(%)  tanh(?72)  +  2  A sech2(/?i  t)sech2(r]2)  ta,nh(r}2) /32  ^  ^ 

^cosh(|)  +  sinh(|)  tanh(/5x  t)  tanh(772)) 

Note  that  setting  S2  =  0  makes  each  of  the  terms  in  the  numerator  of  the  first  line 
of  (4.22)  an  even  function  of  time.  Setting  52  =  0  also  makes  the  denominator  of 
the  first  term  an  even  function.  This  term  is  then  multiplied  by  the  second  line  of 
(4.22),  which  is  a  constant,  sinh(0/2),  times  an  odd  function,  tanh^t),  times  an 
even  function.  Hence,  we  have  several  even  functions  multiplying  an  odd  function, 
making  the  entire  first  term  in  (4.22)  an  odd  function  of  time.  The  second  term  is 
also  seen  to  be  an  odd  function  by  similar  analysis.  As  a  result,  setting  S2  =  0  is  a 
stationary  point  of  (4.20).  To  check  that  this  is  a  minimum,  we  need  to  verify  that 

>  0.  (4.23) 

£2=0 


First  we  note  that 


d 2 
d52 


v2(t)dt 


s2= 0 


+  2  v(t) 


d?v(t) 


dbl 


dt. 


<s2=o 


(4.24) 


Since  dv(t)/d52\s2=o  is  real,  the  first  term  in  (4.24)  is  positive.  The  second  term  con¬ 
tains  v(t),  which  is  real,  and  positive,  and  d?v(t)/d6%\s2=o,  which  can  be  seen  to  also 
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be  positive, 


C fv(t )  _  4 /3|sech (r2)2  tanh(-r2)2-2 sech(r2)2 (l-taah(r2)2) 

^2  J2=0  (cosh(f  )+sinh(f )  tanh(ri)  tanh(T2)) 

4  A sech(ri )2sech(r2 )2  tanh(-r2)2/3f-2  A  sech(Tj )2sech(-r2)2  (l— tanh(r2)2 )/?2  ^ 

(cosh(|)+sinh(|)  tanh(ri)  tanh(r2)) 

(8/3|sech(T2)2tanh(r2)+8  Asech(7-i)2sech(r2)2tanh(r2)<32)  sinh(| )  tanh(ri)(l-tanh(r2)2)/32  ^ 

(cosh(f  )+sinh(f )  tanh(ri)  tanh(r2)) 

(6/3fsech(r1)24-6/32sech(T2)2+6  Asech(ri)2sech(r2)2)  sinh(|)2  tanh(ri)2(l-tanh(r2)2)2/3|  ^ 

(cosh(f)+sinh(f )  tanh(Ti)  tanh(r2)) 

(4  /3^sech(ri)2+4^2sech(r2)2+4  Asech(rx)2sech(r2)2)  sinh(|)  tanh(Tj)  tanh(r2)(l-tanh(T2)2)/32 

(cosh(f  )+sinh(f )  tanh(ri)  tanh(r2)) 

where  r,  =  /?;  t.  Since  each  term  in  (4.24)  is  positive,  the  integral  is  therefore  positive, 
and  62  =  Si  is  indeed  a  minimum. 

To  illustrate  the  energy  reduction,  a  two-soliton  solution  is  shown  in  Fig.  4-6  as 
a  function  of  both  time  and  mutual  separation,  <5i  —  <52-  As  shown  in  the  figure, 
when  the  two  solitons  are  well-separated  in  time,  |5i  —  S2\  3>  0,  the  two  component 
solitons  are  each  distinguishable  as  a  function  of  time  and  the  signal  appears  as  a 
linear  superposition  of  the  two.  In  this  case,  A  «  0  in  Eq.  (4.15).  However,  as  the 
solitons  come  to  interact,  ^  «  S2,  the  nonlinear  cross  term  in  Eq.  (4.15)  becomes 
significant,  and  the  combined  signal  amplitude  decreases. 

The  effect  of  the  energy  reduction  is  illustrated  in  Fig.  4-7  for  several  different 
values  of  the  parameter  f32  holding  fix  fixed,  fix  >  P2.  When  the  smaller  soliton  is 
roughly  the  same  amplitude  as  the  larger,  the  energy  decrease  is  on  the  order  of  10% 
of  the  maximum  signal  energy.  However  as  the  smaller  soliton  reduces  to  roughly 
half  the  amplitude  of  the  larger,  we  see  that  the  energy  reduction  is  almost  30%. 
Significant  energy  reduction  occurs  for  a  fairly  wide  range  of  separations,  indicating 
that  the  modulation  techniques  described  here  could  take  advantage  of  this  reduction. 
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Figure  4-6:  A  two-soliton  solution  to  the  Toda  lattice  as  a  function  of  both  time  and 
mutual  separation,  Si  —  <52. 

4.5  Gain  Normalization 

In  any  practical  communication  context,  a  modulation  system  must  be  able  to  combat 
the  presence  of  an  unknown  channel  gain  due  to  fluctuations  in  the  channel  character¬ 
istics.  This  is  a  potential  drawback  of  using  Toda  lattice  soliton  carrier  signals  since 
these  solitons  have  a  specific  relationship  between  the  amplitude  and  pulse-width 
which  must  be  preserved.  If  the  soliton  signal, 

s(t)  =  /32sech2(/3t)  (4.25) 

is  transmitted  through  a  channel  and  arrives  at  the  receiver  with  an  unknown  gain, 
r(t )  =  as(t),  then  the  soliton  dynamics  of  the  signal  s(t)  can  no  longer  be  observed 
from  processing  the  signal  r(t).  In  general,  the  signal  r(t)  will  give  rise  to  both  soliton 
and  non-soliton  components,  where  the  soliton  component  corresponds  to  a  different 
soliton  parameter  ^  ^  S3. 

Many  communication  systems  combat  gain  fluctuations  by  using  a  form  of  auto¬ 
matic  gain  control  (AGC)  to  dynamically  adjust  the  gain  of  a  pre-amplifier  in  the 
receiver,  r(t)  =  ay s(t).  To  demonstrate  the  feasibility  of  such  AGC  techniques  for 
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Figure  4-7:  Normalized  signal  energy  for  a  two-soliton  solution  to  the  Toda  lattice 
holding  pi  fixed  for  3  values  of  /?2-  The  signal  energy  is  normalized  by  the  maximum 
signal  energy  of  the  separated  solitons. 

the  soliton  modulation  systems  investigated  in  this  chapter,  we  will  explore  the  ef¬ 
fect  of  an  unknown  gain  on  a  single  soliton  as  it  is  processed  by  the  Toda  lattice. 
Since  the  effect  of  processing  the  soliton  will  correspond  to  a  simple  time-delay  only 
when  the  unknown  gain  has  been  removed,  cry  =  1,  an  AGC  system  might  exploit 
differences  between  the  input  and  the  processed  waveforms,  and  adjust  the  gain  7 
until  the  processed  waveform  is  a  pure  time-delay  of  the  input.  If  both  the  input  and 
processed  signals  are  made  zero-mean,  then  the  energy  in  the  difference  between  the 
input  and  the  output  when  the  largest  soliton  in  the  output  is  time  aligned  with  the 
input  corresponds  to 

min  J (( r(t )  -  mT)  -  ( rn(t  -  r)  -  mrn))2dt  =  (4.26) 

mm£r  +  ETn  -  2  J (r(t)  -  mT)(rn(t  -  r)  -  mrn)dt, 

where  mr  and  mrn  are  the  means  and  Er  and  Ern  are  the  energies  of  the  input  and 
output,  respectively.  When  the  gain  has  been  properly  adjusted,  cry  =  1,  this  error 
will  be  a  minimum.  If  each  signal  is  energy  normalized,  Er  —  Ern  =  1,  then  minimiz¬ 
ing  Eq.  (4.26)  amounts  to  maximizing  the  cross-correlation  between  the  normalized 
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Figure  4-8:  Normalized  cross-covariance  of  input  with  the  processed  signal  as  a  func¬ 
tion  of  the  unknown  gain,  a. 

signals.  The  cross-correlation  between  the  two  normalized  signals  can  be  calculated 
as  the  peak  of  the  cross-covariance  function, 

Cr, m(T)  =  / (r(*)  “  mT){rn{t  -  r)  -  mrn)dt.  (4.27) 

llirJl/rn  J 

Since  rn(t)  =  r(t  —  r)  only  when  cry  =  1,  then  by  the  Schwartz  inequality,  the  peak 
of  CV>rn(r)  <  1,  with  equality  only  achieved  when  when  07  =  1.  As  an  example, 
in  Fig.  4-8,  the  peak  of  the  cross-covariance  between  the  gain  adjusted  soliton  signal 
r(t )  =  ajs(t),  and  the  processed  signal  at  the  n-th  node  of  the  lattice  rn(t)  is  shown  as 
a  function  of  the  unknown  gain  cry.  For  this  example,  j3  =  sinh(2),  and  0.1  <  a  <  3. 
As  shown,  the  normalized  cross-covariance  has  a  unique  maximum  of  C  —  1,  when  the 
input  to  the  Toda  lattice  has  been  properly  rescaled,  cry  =  1.  This  gives  an  indication 
that  AGC  techniques  based  on  a  feedback  of  the  normalized  cross-covariance  between 
the  input  and  the  processed  signal  might  be  effective  in  combating  unknown  channel 
gain.  However,  a  variety  of  issues  including  how  such  AGC  might  be  performed  for 
multi-soliton  signals  or  modulated  multi-soliton  carriers  remains  an  interesting  area 
for  future  research. 
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4.6  Other  Potential  Applications 


We  have  explored  in  a  very  preliminary  manner  a  modulation  technique  using  soli- 
tons  as  carrier  waveforms  and  nonlinear  evolution  equations  as  tuned  transmitters 
and  receivers.  There  are  several  interesting  aspects  of  such  techniques  that  warrant 
further  study.  For  example,  noting  that  the  same  nonlinear  network  can  support  a 
spectrum  of  solitons  ranging  from  those  with  small  amplitudes  and  wide  pulse- widths, 
to  relatively  large  amplitudes  and  narrow  pulse- widths,  indicates  the  potential  for  op¬ 
erating  essentially  the  same  modulator-demodulator  networks  at  variable  data  rates 
depending  on  the  bandwidth  requirements.  In  the  case  of  on-off  keying  modulation, 
where  a  bit  could  be  indicated  by  the  presence  or  absence  of  a  soliton,  a  tradeoff 
may  be  made  between  the  data  rates  and  the  power  in  the  transmitted  signal.  When 
communication  bandwidth  requirements  are  low,  small  amplitude,  wide  solitons  may 
be  used  and  the  necessary  transmitted  power  is  low.  When  bandwidth  requirements 
increase,  so  does  the  requisite  output  power  as  narrower  larger-amplitude  solitons 
are  used.  This  same  analogy  holds  for  either  the  amplitude  or  position  modulation 
schemes  addressed  in  this  chapter. 

Such  signaling  techniques  may  also  prove  useful  in  the  context  of  recent  advances 
in  multi-resolution  signal  representations.  In  a  hierarchical  modulation  technique 
reminiscent  of  fractal  modulation  [77]  a  multi-resolution  representation  of  a  signal 
might  be  transmitted  such  that  each  scale  of  the  signal  were  encoded  in  a  soliton  of  a 
different  amplitude.  If  each  of  these  soliton  waveforms  are  combined  and  time-aligned 
according  to  the  nonlinear  superposition  of  the  network,  the  response  of  the  receiver 
network  to  this  signal  will  be  a  gradual  separation  of  each  of  the  component  soli¬ 
ton  waveforms.  In  this  manner,  compactly  represented  signals  might  be  constructed 
whereby  varying  amounts  of  processing,  or  equivalently,  longer  delay,  yield  signal 
representations  of  varying  fidelity.  As  the  received  waveform  is  processed,  the  in¬ 
formation  present  in  the  higher-amplitude,  faster  solitons  emerges  quickly  from  the 
composite  waveform  and  may  be  decoded.  As  the  message  is  passed  further  down 
the  chain,  the  information  in  the  next  set  of  solitons  may  be  decoded,  and  so  on. 
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This  type  of  modulation  may  be  useful  in  a  variety  of  communications  or  broadcast 
contexts  which  contain  a  large  number  of  receivers  of  variable  processing  power.  This 
may  also  be  useful  in  a  signaling  context  with  messages  of  variable  priority.  In  either 
case,  successively  more  processing  power,  more  accurate  hardware,  or  more  process¬ 
ing  delay  may  be  required  to  extract  the  messages  encoded  on  the  smallest  amplitude 
solitons.  This  could  correspond  to  either  users  of  lower  priority  in  a  priority-based 
multiple  access  communications  scenario,  or  to  an  incremental  fidelity  enhancement 
as  envisioned  for  progressive  transmission  or  embedded  coding  systems. 

An  essential  ingredient  in  the  modulation  concepts  explored  in  this  chapter  is  the 
nonlinear  superposition  of  the  component  soliton  signals  for  combined  packetization 
and  energy  reduction  of  the  multi-soliton  signal.  As  we  shall  see  in  Chap.  6  it  is  during 
this  nonlinear  superposition  that  the  susceptibility  to  errors  in  the  multi-soliton  signal 
may  be  traded  off  between  the  component  solitons.  In  fact,  over  certain  ranges  of 
parameter  values,  in  addition  to  a  net  reduction  in  the  transmitted  signal  energy, 
there  can  be  net  increase  in  fidelity  of  the  recovered  component  messages. 

A  technique  for  modulation  of  information  on  soliton  carriers  was  also  proposed 
by  Hirota,  Suzuki  and  Yoshikawa  in  [64]  and  [65].  In  their  work,  an  amplitude  and 
phase  modulation  of  a  two-soliton  solution  to  the  Toda  lattice  were  presented  as  a 
technique  for  private  communication.  Although  their  signal  generation  and  process¬ 
ing  methods  relied  on  an  inexact  phenomenon  known  as  recurrence,  the  modulation 
paradigm  they  presented  is  essentially  a  two-soliton  version  of  the  carrier  modulation 
paradigm  presented  in  this  chapter.  Rather  than  generating  multi-soliton  signals  in 
the  lattice  directly,  a  Toda  lattice  circuit  was  driven  with  a  sinusoidal  signal.  For 
certain  frequencies,  a  two-soliton-like  solution  is  induced  into  the  line.  Each  of  the 
component  solitons  were  separately  modulated,  and  an  appropriate  length  of  the  Toda 
chain  was  used  to  allow  the  solitons  to  recombine.  This  signal  was  then  considered 
the  transmitted  signal,  to  be  separated  by  another  Toda  circuit  receiver.  This  overall 
process  is  illustrated  in  Fig.  4-9. 

These  methods  were  proposed  as  techniques  for  secure  communication,  the  jus¬ 
tification  being  that  the  nonlinear  effects  of  the  soliton  superposition  would  cause 
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Figure  4-9:  Schematic  diagram  of  AM-like  modulation  of  Hirota  and  Suzuki.  Redrawn 
based  on  [64,  65] 


enough  spectral  mixing  of  the  component  message  waveforms  that  decoding  without 
the  correct  nonlinear  circuit  would  be  difficult.  However,  in  their  experiments  the 
message  waveforms  were  pure  sinusoids  and  as  they  point  out  from  observations  of 
the  signal  spectrum,  these  messages  could  be  recovered  by  using  a  narrow  band  filter 
near  the  carrier  frequency. 

Before  such  methods  could  be  considered  for  private  or  LPI  (low  probability  of 
intercept)  communications  in  practice,  as  is  the  case  for  direct-sequence  spread  spec¬ 
trum  systems,  an  additional  level  of  security  would  have  to  be  employed.  Typically  in 
direct  sequence  spread  spectrum,  the  binary-valued  symbol  sequence  is  multiplied  by 
a  pseudo-random  binary-valued  sequence  known  only  to  the  transmitter  and  receiver. 
In  the  soliton  position  modulation  example,  this  could  correspond  to  a  pseudo-random 
perturbation  to  the  relative  phasing  of  the  component  solitons.  Again,  the  efficacy  of 
such  techniques  would  have  to  be  explored  in  detail. 

The  issue  of  the  relative  phasing  of  the  soliton  carrier  signals  brings  to  light 
another  interesting  analogy  with  more  traditional  modulation  techniques.  An  issue  of 
interest  in  modern  multi-user  communication  systems  is  the  relative  performance  of 
time-division-multiple-access  (TDMA)  versus  code-division-multiple-access  (CDMA) 
systems  and  the  corresponding  complexity  of  such  systems.  If  the  relative  phases 
of  the  component  solitons  in  the  soliton  PPM  system  are  chosen  such  that  they 
overlap  temporally,  then  at  least  in  the  linear  limit,  a  form  of  CDMA  is  approached. 
Conversely,  if  the  relative  phasing  is  chosen  such  that  the  solitons  do  not  overlap, 
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then  a  TDMA  system  is  employed.  If  the  phasing  is  chosen  such  that  only  the 
highest  amplitude  soliton  is  separate  from  the  composite,  a  hybrid  TDMA/CDMA- 
like  system  is  realized.  In  addition,  a  variety  of  digital  communications  techniques 
could  be  implemented  using  soliton  carrier  signals.  However,  the  potential  benefits 
of  such  techniques  remain  an  area  of  future  study. 
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Chapter  5 


Noise  Dynamics  in  Soliton  Systems 


We  have  seen  that  soliton  solutions  to  nonlinear  evolution  equations  provide  an  in¬ 
triguing  and  potentially  viable  alternative  to  sinusoidal  signals  as  carrier  waveforms 
in  the  broad  context  of  communication.  In  order  to  address  the  efficacy  of  the  mod¬ 
ulation  techniques  presented  in  Chapter  4,  accurate  models  are  needed  for  the  effects 
of  random  fluctuations  on  the  dynamics  of  soliton  systems.  Such  disturbances  could 
take  the  form  of  additive  noise  or  interference  or  convolutional  corruption  incurred 
during  terrestrial  or  wired  transmission,  circuit  thermal  noise,  or  modeling  errors  due 
to  system  deviation  from  the  idealized  soliton  dynamics.  A  fundamental  property  of 
solitons  that  has  made  them  appeal  to  such  a  broad  area  of  science  and  engineering 
is  that  they  are  stable  in  the  presence  of  a  variety  of  disturbances.  From  their  early 
discovery  in  the  Scottish  canal,  it  was  the  inherent  stability  of  one  such  solitary  wave 
that  first  attracted  John  Scott-Russell  and  influenced  his  study  of  the  phenomena 
for  many  years  [57].  For  several  decades,  however,  solitary  waves  were  considered 
to  be  fluke  solutions  to  a  special  class  of  nonlinear  evolution  equations  and  it  was 
conjectured  that  two  solitary  waves  would  annihilate  one  another  due  to  nonlinear 
interaction.  It  was  after  early  computer  experiments  by  Perring  and  Skyrme  [57] 
with  the  sine-Gordon  equation,  and  later  numerical  studies  of  KdV  by  Zabusky  and 
Kruskal  that  revealed  the  contrary,  prompting  Zabusky  and  Kruskal  to  term  these 
stable  solitary  waves  solitons  since  they  “pass  through  one  another  without  losing 
their  identity”  [78]. 
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In  the  literature,  several  measures  of  the  stability  or  robustness  of  such  nonlin¬ 
ear  systems  and  their  soliton  solutions  have  been  investigated.  In  addition  to  early 
numerical  work,  studies  like  [51]  have  empirically  investigated  the  stability  of  soli¬ 
ton  solutions  in  the  presence  of  additive  corruption.  Lindgren  and  Buratti  [43]  have 
studied  analytically  the  stability  of  solitons  in  the  sine-Gordon  equation  through  a 
linearization  about  a  known  soliton  solution.  Some  success  has  also  been  achieved  for 
such  a  study  of  one-dimensional  nonlinear  lattices  by  Flytzanis  et  al.  [17].  Many  forms 
of  perturbation  theory  and  approximate  linear  analysis  have  also  been  applied  to  the 
nonlinear  Schrodinger  equation,  demonstrating  the  viability  of  proposed  telecommu¬ 
nications  systems  as  well  [40].  In  [5],  Benjamin  demonstrates  the  Lyapunov  stability 
of  a  single  soliton  in  the  fully  nonlinear  KdV  equation. 

There  has  also  been  increasing  interest  in  the  solvability  of  soliton  equations  in 
the  presence  of  additive  noise.  This  area  of  the  literature  concerns  systems  such  as 
the  “stochastic  KdV  equation”  which  is  a  rather  restrictive  setting  in  which  additive 
noise  is  additive  and  is  a  function  of  time,  while  remaining  a  constant  function  of 
space, 

Uf  4“  QlLUx  4"  'U'xxx  —  n(t).  (5.1) 

This  system  can  be  shown  to  possess  an  exact  soliton  solution  with  a  phase  drift  that 
is  given  by  a  Wiener  process,  when  n(t)  is  a  stationary  white  Gaussian  process  [44,  74]. 

With  the  development  of  the  inverse  scattering  framework  and  the  discovery  that 
many  soliton  systems  were  conservative  Hamiltonian  systems,  many  of  the  questions 
regarding  the  stability  of  soliton  solutions  are  readily  answered.  For  example,  since 
the  eigenvalues  of  the  associated  linear  operator  remain  unchanged  under  the  evolu¬ 
tion  of  the  dynamics,  then  any  solitons  that  are  initially  present  in  a  system  must 
remain  present  for  all  time,  regardless  of  their  interactions.  Similarly,  the  dynam¬ 
ics  of  any  non-soliton  components  that  are  present  in  the  system  are  orthogonal  to 
the  dynamics  of  the  solitons.  However,  in  the  communication  scenario  discussed  in 
Chapter  4,  soliton  waveforms  are  generated  and  then  propagated  over  a  noisy  chan¬ 
nel.  During  transmission,  these  waveforms  are  susceptible  to  additive  corruption 
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from  the  channel.  When  the  waveform  is  received  and  processed,  the  inverse  scat¬ 
tering  framework  can  provide  useful  information  about  the  soliton  and  noise  content 
of  the  received  waveform.  Since  the  time  evolution  of  the  dynamics  of  such  systems 
leaves  the  eigenvalues  of  a  linear  operator  unchanged,  investigations  into  their  stabil¬ 
ity  involves  a  study  of  the  spectra  of  random  linear  operators.  Although  there  has 
been  recent  work  in  this  area  [8,  54],  little  that  has  been  published  applies  directly  to 
our  communication  context. 

In  this  chapter,  we  will  assume  that  soliton  signals  generated  in  the  communication 
contexts  of  Chapter  4  have  been  transmitted  over  an  additive  white  Gaussian  noise 
channel.  We  can  then  consider  the  effects  of  additive  corruption  on  the  processing 
of  soliton  signals  with  their  nonlinear  evolution  equations.  Two  general  approaches 
are  taken  to  this  problem.  The  first  primarily  deals  with  linearized  models  and 
investigates  the  dynamic  behavior  of  the  noise  component  of  signals  comprising  an 
information  bearing  soliton  signal  and  additive  noise.  The  second  approach  is  taken  in 
the  framework  of  inverse  scattering  and  is  based  on  some  results  from  random  matrix 
theory.  Although  the  analysis  techniques  developed  in  this  chapter  are  applicable 
to  a  large  class  of  soliton  systems,  we  focus  our  attention  on  the  Toda  lattice  as  a 
representative  example. 


5.1  Toda  Lattice  Small  Signal  Model 

If  a  signal  that  is  processed  in  a  Toda  lattice  receiver  contains  only  a  small  amplitude 
noise  component,  then  the  dynamics  of  the  receiver  can  be  approximated  by  a  small 
signal  model.  Starting  with  the  nonlinear  transmission  line  model, 

A  1„(1  +  V„(t))  =  -  2V„(t)  +  14+i(t)),  (5.2) 
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and  using  the  approximation,  ln(l  +  x)  «  x,  we  see  that  the  lattice  equations  can  be 
approximately  described  by  the  linear  lattice  equations 

=  -£(Vi-i (t)  -  2 V„(t)  +  K+1(i)),  (5.3) 

when  the  amplitude  of  Vn(t)  is  appropriately  small.  Since  this  model  is  linear,  we 
may  decompose  solutions  into  harmonic  components  of  the  form 


Vn(t)  =  V+ej{kn~wt)  +  y_ertfcn+wt). 


(5.4) 


From  Eqs.  (5.3)  and  (5.4)  we  see  that  the  frequency  of  a  single  forward  propagating 
solution  must  satisfy  the  dispersion  relation 


-u)2  = 


1 

LC 


(e~jk  -2  +  ejk), 


(5.5) 


which  reduces  to 


u>  — 


2sin(k/2) 

Vlc 


Therefore  the  lattice  is  dispersive,  with  frequency-dependent  velocity, 


c(k)  = 


2sin(fc/2) 
ky/LC  ’ 


(5.6) 


(5.7) 


or 


CM  =  r— — 

‘1  cin  L 


u 


28m-1  ( 'uy/LC/2 ))' 
Note  that  we  can  also  write  the  dispersion  relation  as 


(5.8) 


k  =  2  sin  1  ( 'uVLC/2 )  , 


(5.9) 


from  which  k  is  only  real  if 


M  < 


Vlc’ 


(5.10) 
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for  which  there  are  propagating  waves.  When  u>  is  outside  this  region,  the  wavenum¬ 
ber,  k,  is  complex,  corresponding  to  evanescent  waves  of  the  form 


uj  =  —j=  cosh(Im(fc)/2),  Re(A;)  =  7r. 
v  LC 

These  solutions  decay  as  they  pass  through  the  lattice, 


(5.11) 


| y  j  _  |j/^|e~2cosh  1(wvTC/2)n 


(5.12) 


which  for  u>  2  j\jLC  corresponds  to 


v"=^(^y 


(5.13) 


If  we  consider  processing  signals  with  an  infinite  linear  lattice  and  obtain  an  input- 
output  relationship  where  a  signal  is  input  at  the  zeroth  node  and  the  output  is  taken 
as  the  voltage  on  the  iV-th  node,  we  see  that  the  input-output  frequency  response  of 
the  system  can  be  given  by 


HN(jv)  = 


g— 2j sin- 1  (u) VXC/2) N ^  |w|  <  2/^/ZC 

e{jn-2cosh~1(aVLC/2)]N  eJge 


(5.14) 


As  shown  in  Fig.  5-1,  the  lattice  behaves  as  a  low  pass  filter,  and  for  N  1,  ap¬ 
proaches 


2  _  J  1,  M  <  u>c  =  2/y/LC 


0,  else. 


(5.15) 


5.2  Linearized  Model 


Our  small  signal  model  indicates  that  in  the  absence  of  solitons  in  the  received  signal, 
small  amplitude  noise  will  be  processed  by  a  low  pass  filter.  If  the  received  signal 
also  contains  solitons,  then  the  small  signal  model  of  Eq.  (5.3)  will  no  longer  hold.  A 
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Figure  5-1:  The  log  magnitude  of  the  frequency  response  from  the  input  node  to  the 
iV-th  node  as  a  function  of  normalized  frequency.  As  indicated,  the  response  rapidly 
drops  off  as  a  function  of  N  for  u  >  uc. 


linear  small  signal  model  can  still  be  used  if  we  linearize  Eq.  (5.2)  about  the  known 
soliton  signal  of  the  form 

(5.16) 

Assuming  that  the  solution  contains  a  single  soliton  in  small  amplitude  noise,  Vn(t)  = 
Sn(t)  +  u„(t),  we  can  write  Eq.  (5.2)  as 

<P  1 

—  In  (1  +  Sn  +  vn)  =  ~  2Sn  +  Sn+ 1  +  un_i  -  2vn  +  u„+i).  (5.17) 


Sn(t )  =  fl2sech2  I  pn 


Qt 

Vlc 


After  factoring  the  argument  of  the  logarithm  and  cancelling  terms  from  both  sides 
that  correspond  to  the  known  soliton  solution,  we  have  an  exact  equation  that  is 
satisfied  by  the  non-soliton  component, 

dfi  ln  +  1  +  Sn{t))  =  LC{Vn~l{t)  ~  2vn(t)  +  Vn+l(t)),  (5.18) 
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Figure  5-2:  Receiver  model  comprising  a  low  pass  filter  followed  by  a  Toda  lattice 
circuit. 


which  can  be  viewed  as  the  fully  nonlinear  model  with  a  time-varying  parameter, 
(l  +  Sn(t)).  As  a  result,  over  short  time  scales  relative  to  Sn(t),  we  would  expect  this 
model  to  behave  in  a  similar  manner  to  the  small  signal  model  of  Eq.  (5.3).  With 
vn(t)  (1  +  Sn(t)),  we  obtain 

^TTsjf)  “  “  2v"(t)  +  (519) 

When  the  contribution  from  the  soliton  is  small,  Eq.  (5.19)  reduces  to  the  linear 
system  of  Eq.  (5.3).  We  would  therefore  expect  that  both  before  and  after  a  soliton  has 
passed  through  the  lattice,  the  system  essentially  low  pass  filters  the  noise.  However, 
as  the  soliton  is  processed,  there  will  be  a  time- varying  component  to  the  filter. 


5.3  Simulation  of  the  Lattice  in  Noise 

To  confirm  the  intuition  developed  through  small  signal  analyses,  we  have  simulated 
the  fully  nonlinear  dynamics.  As  indicated  in  Chap.  3,  we  will  work  with  a  finite- 
length  lattice  which  is  terminated  with  its  linearized  impedance.  We  then  focus  on 
the  dynamics  of  the  small  amplitude  noise  component  in  the  response  of  the  lattice 
to  a  signal  comprising  a  single  soliton  in  white  Gaussian  noise  with  noise  power  JV0. 

Our  primary  interest  in  this  chapter  is  to  characterize  the  effects  of  additive  noise 
in  a  receiver  for  a  potential  soliton  modulation  system.  Since  the  bandwidth  limita¬ 
tions  of  the  receiver  for  any  of  the  communications  scenarios  discussed  in  Chap.  4  will 
restrict  the  possible  range  of  soliton  parameters,  without  loss  of  generality,  we  may 
assume  that  the  receiver  comprises  a  low  pass  filter  followed  by  a  Toda  lattice  circuit 
as  shown  in  Fig.  5-2.  We  also  assume  that  the  bandwidth,  27r/A,  of  the  low  pass 
filter  in  Fig.  5-2  is  wide  enough  to  pass  the  soliton  component  of  r(t)  completely.  The 
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input  to  the  Toda  lattice  circuit,  Vq  (t),  then  comprises  the  soliton  signal  in  low  pass 
Gaussian  noise.  Simulations  are  performed  using  a  Runge-Kutta  integration  routine 
with  a  fixed  step  size,  A.  To  model  the  effects  of  the  noise,  an  i.i.d.  Gaussian  random 
sequence,  w(kA )  ~  N( 0,  cr^),  was  added  to  the  samples  of  the  input  sequence  V0(kA) 
resulting  in  an  effective  white  noise  power  of  N0  =  A o^. 

The  circuit  equations  governing  the  resistance-terminated  lattice  are  given  in 
Eq.  (5.2)  for  n  <  N.  At  the  termination,  n  =  N,  we  have 


VN(t)  = 


vk-m  -  v„(t)  vwtn  n . ,,  ,  i 'm 

- LC - ~RC~ j  (1  +  lV(t>)  +  T7W 


(5.20) 


Writing  Eq.  (5.2)  and  (5.20)  as  2 N  first-order  differential  equations,  we  obtain 


MV 

wn(t ) 

wN(t) 


Wn(t),  1  <  n  <  N  (5.21) 

+ 1  i  < « <  n 


(Vn-l(t)  -  Vw(t) 

v  LC 


RC  j 


(1  +  Viv(£))  + 


i  +  vN(ty 


(5.22) 


From  our  linearized  analyses,  we  anticipate  that  the  response  of  the  lattice  to  a 
soliton  in  small  amplitude  Gaussian  noise  will  essentially  comprise  the  unperturbed 
soliton  with  an  additional  small  amplitude  low  pass  Gaussian  component.  In  Fig.  5-3 
we  show  the  response  of  the  lattice  to  a  single  soliton  at  20dB  signal-to-noise  ratio, 
where  the  SNR  is  calculated  using 


SNR  =  10  log 


(5.23) 


As  expected,  the  response  to  the  lattice  is  essentially  the  unperturbed  soliton 
with  an  additional  low  pass  perturbation.  The  noise  component  of  the  response  on 
the  third  node  of  the  lattice  is  indicated  in  Fig.  5-4.  The  spectrum  of  the  noise 
remains  essentially  flat  over  the  bandwidth  of  the  soliton  and  is  attenuated  out  of 
band. 
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SNR  =  20  dB 


Figure  5-3:  Response  to  a  single  soliton  with  j3  =  sinh(l)  in  20  dB  Gaussian  noise. 
The  spectrum  of  the  noise  process  is  flat  out  to  half  the  sample-rate  of  the  integration 
routine.  The  corresponding  in-band  SNR  is  approximately  24  dB. 

As  the  SNR  decreases,  and  correspondingly  vn(t)/(  1  4-  Sn(t))  becomes  significant 
relative  to  unity,  the  linear  approximation  no  longer  applies  and  the  noise  dynamics 
must  be  described  with  the  fully  nonlinear  equations.  As  shown  in  Fig.  5-5,  the 
response  to  the  noise  term,  vn(t),  contains  small-amplitude  solitons  as  well  as  non- 
soliton  components.  The  behavior  of  solitons  induced  by  noise  is  better  handled  in 
the  framework  of  inverse  scattering,  which  is  the  topic  of  Sec.  5.6. 

5.4  Noise  Correlation 

The  statistical  correlation  of  the  system  response  to  the  noise  component  can  also 
be  estimated  from  our  linear  analyses.  From  Sec.  5.1,  the  small  signal  model  for  the 
nonlinear  lattice  approximately  satisfies  the  linear  equations, 

Q p  ^ 

^n(t)  =  -^(w„_i(i)  -  2 vn{t)  +  Vn+1(t)),  (5.24) 
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Figure  5-4:  Noise  response  to  a  single  soliton  in  20  dB  Gaussian  noise  as  viewed  from 
the  third  node  in  the  lattice. 

which  has  a  magnitude-squared  frequency  response  at  the  n-th  node,  n  3>  1,  of 


1,  H  <  2/VZC 

0,  else. 


(5.25) 


Therefore,  vn{t)  is  zero  mean  and  has  an  auto-correlation  function  given  by 


Rn,n(r)  =  E{vn{t)vn(t  +  r)}  «  N0 


sin  (u;cr) 


7TT 


(5.26) 


and  a  variance  a%n  fa  N0uc/ir,  for  1. 

Although  the  autocorrelation  of  the  noise  at  each  node  is  only  affected  by  the 
magnitude  response  of  Eq.  (5.14),  the  cross-correlation  between  nodes  is  also  affected 
by  the  phase  response.  The  cross-correlation  between  nodes  m  and  n  is  given  by 


Rm,n(T)  =  hm{r)  *  K{-t)  =  hm(r)  *  hm(-r)  *  hn_m(-r),  (5.27) 

where  hm(r)  is  the  inverse  Fourier  transform  of  Hm(juj)  in  Eq.  (5.14).  Since  hm{r)  % 
hm{—T)  approaches  the  impulse  response  of  an  ideal  low  pass  filter  for  m  »  1,  we 
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SNR  =  10  dB 


Figure  5-5:  Response  to  a  single  soliton  with  /?  =  sinh(l)  in  10  dB  Gaussian  noise, 
have 

Rm,niT)  ~  N0Sm{u}cT)  *  h„_m(r).  (5.28) 

7TT 

In  Fig.  5-6,  Rm^n(j)  is  shown  for  n  >  m  »  1.  Note  that  for  u  small  in  Eq.  (5.14), 
sin-1(u;/u;c)  «  cu/cuc,  and  the  lattice  looks  like  a  pure  delay  of  a  =  2 (n  —  m)/ujc, 
corresponding  to  Rm,n(T )  ~  sin(cuc(T  —  a))/(7r(r  —  a)).  This  approximation  is  only 
valid  in  the  low  frequency  limit  and  corresponds  to  the  diagonal  translation  of  the 
largest  lobe  of  Rm,n{r)  in  Fig.  5-6. 

For  small  amplitude  noise,  the  correlation  structure  can  be  examined  through  the 
linear  lattice,  which  acts  as  a  dispersive  low  pass  filter.  A  corresponding  analysis  of 
the  nonlinear  system  in  the  presence  of  solitons  becomes  prohibitive  in  closed  form. 
However  we  can  explore  the  analyses  numerically  by  linearizing  the  dynamics  of  the 
system  about  the  known  soliton  trajectory. 

To  examine  the  correlation  structure  in  the  presence  of  soliton  components,  we 
use  the  state  space  framework  of  linear  dynamic  systems.  The  state  space  model 
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Figure  5-6:  Cross-correlation,  /^^(r),  between  the  m-th  and  the  n-th  node  voltages 
in  the  linearized  lattice. 

comprises  a  linear  system  of  finite  dimension,  N,  with  state  vector 

x(t)  =  [x0(t),...,a;v-i(i)]T 

and  dynamics  in  state  space  form, 

x(t)  =  A(t)x(t)  +  b(t)u(t),  (5.29) 

where  A(t)  is  an  N  x  N  state  transition  matrix,  b(t)  is  an  N  x  1  vector,  and  u(t)  is  a 
scalar  input.  We  consider  u(t)  to  be  a  zero-mean,  white  Gaussian  noise  process  with 
noise  power  a\. 

For  a  linear  system  of  the  form  (5.29),  the  state  covariance  matrix, 

P(t)  =  E  {(x(i)  -  £{x(t)})(x(i)  -  E{x(i)})T}  ,  (5.30) 

satisfies  the  following  differential  equation  [19] 

P{t)  =  A(t)P(t)  +  P(t)A(t)T  +  b{t)a2ub(t)r.  (5.31) 
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To  limit  the  number  of  state  variables  in  x,  we  again  terminate  the  nonlinear 
lattice  with  its  linearized  impedance.  If  we  assume  that  the  input  to  the  nonlinear 
lattice  is  of  the  form, 

Vin(t)  =  V0°(t)+u(t),  (5.32) 

where  u(t)  is  a  small  amplitude  white  Gaussian  noise  process,  and  V^t)  corresponds 
to  a  known  soliton  input,  we  may  linearize  the  dynamics  about  the  known  response 
of  the  system. 

By  seeking  a  response  Vn(t)  =  V°(t)  +  vn(t),  and  Wn(t)  =  Vn(t)  =  W„(t)  +  wn{t), 
where  V^(t)  and  W®(t)  are  the  known  responses  to  the  input  V^t),  and  vn  (t)  and 
wn(t)  are  the  responses  to  the  small  amplitude  noise  component,  u(t),  we  obtain 


V°  +  Vn  = 
W%  +  Wn  = 


+ 


K  +  «£, 

W  +  Wnf 
l  +  V°+vn 
1  +  V®  +  Vn 
LC 


1  <n<N 
1  <  n  <  N 


(kn+i  2Vn  +  +  Un+i  2vn  +  vn _ 


i)> 


/  Vn-  i  —  Vn  +  vn- i  ~  vn 
V  LC 


WN  +  wN\ 
RC  ) 


(1  +  WN  +  wn ) 


(WN  +  wn)2 
1  +  VN  +  vn  ’ 


(5.33) 


where  =  V0° (t)  and  v0  =  u(t).  Cancelling  terms  that  correspond  to  the  known 
input  and  response  and  terms  higher  than  first  order,  we  obtain 


Vn  =  t£>®,  1  <  n  <  N 

XV  -  2VW°  +  14-1  (X)2  2(1  +  VJ) 

LC 


Wn  = 


+ 

wN  = 

+ 


1  +  X 

LC 


(?Wl  +  Uji-l)  + 


(i  +  v»y  lc 

2W° 


i  +  vs 


Wr 


1  <  n  <  N 


VN-1  -  VS  W% 


LC  RC 
vn- i  ~  vn  wn  ~ 


vN 


2W% 


1  “JV  .  0 

lc  rc\  (1 + Vw)  +  rry? 


(5.34) 
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Eqs.  (5.34)  can  be  written  in  the  form, 


v(f) 

w(t) 


—  A(t) 


v(t) 

w(t) 


+  b  (t)u(t), 


(5.35) 


where  v(t)  =  [m(t), . . .  ,Ujv(i)]T,  and  w (t)  =  [wi (t), . . . ,  w^(t)]T. 

From  our  earlier  linearized  analyses,  the  linear  time-varying  small  signal  model 
can  be  viewed  over  short  time  scales  as  a  linear  time-invariant  chain,  with  a  slowly 
varying  parameter.  The  resulting  input-output  transfer  function  can  be  viewed  as  a 
low  pass  filter  with  time  varying  cutoff  frequency  equal  to  u ;c  when  a  soliton  is  far  from 
the  node,  and  to  w0y  1  +  V®  as  a  soliton  passes  through.  Thus,  we  would  expect  the 
variance  of  the  node  voltage  to  rise  from  a  nominal  value  as  a  soliton  passes  through. 

We  numerically  integrate  the  corresponding  Riccati  equation,  (5.31),  for  the  node 
covariance  and  in  Fig.  5-7,  the  resulting  variance  of  the  noise  component  on  each  node 
is  shown.  In  this  example,  the  input  to  the  lattice  was  a  periodically  repeated  single 
soliton  with  an  initial  SNR  of  30  dB.  Since  the  lattice  was  assumed  initially  at  rest, 
there  is  a  startup  transient,  as  well  as  an  initial  spatial  transient  at  the  beginning 
of  the  lattice,  after  which  we  see  that  the  variance  of  the  noise  is  amplified  from  the 
nominal  variance  as  each  soliton  passes  through,  confirming  our  earlier  intuition. 


Figure  5-7:  The  variance  of  each  node  voltage  as  a  function  of  time. 
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5.5  Noise  Dynamics  for  the  Diode  Ladder  Circuit 


The  prior  analyses  developed  in  this  chapter  apply  to  the  Toda  lattice  when  the  noise 
component  of  the  signal  can  be  considered  small  in  comparison  to  the  remaining 
arguments  of  the  logarithm  in 

^'a(1+W]= ik{K-i{t)  - 2K(() + wt))-  (5-36) 

When  Vn(t )  is  small  as  compared  with  Vo,  then  Eq.  (5.36)  behaves  like  a  linear  LC 
ladder.  However,  for  the  diode  ladder  circuit  which  satisfies, 

+  =^(V,W-2/»(t) +  /„+,(*)),  (5.37) 

for  a  similar  small  signal  analysis,  In(t )  would  have  to  be  small  in  comparison  to 
the  saturation  current,  Is.  This  would  either  require  diodes  with  an  unusually  large 
saturation  current  or  very  small  signal  levels. 

For  a  solution  containing  a  soliton  signal,  J„(t),  and  a  small  amplitude  noise  signal, 
in(t),  an  exact  expression  for  the  small  amplitude  component  is  given  by 

dP ln  (X  +  In(if+/S)  =  “  2ln ^  +  ?n+1  (5-38^ 

The  linearization  that  results  from  the  assumption  that  the  current  in(t )  is  small 
in  comparison  to  the  saturation  current  is  tantamount  to  replacing  the  diodes  with 
their  equivalent  linearized  resistance,  Req  =  vt/Is,  which  is  on  the  order  of  MO.  The 
resulting  small  signal  model  has  a  low  pass  characteristic  with  a  cutoff  frequency  of 
uc  =  a /Req  «  8kHz  for  the  range  of  circuit  parameters  used  in  Chap.  3.  Since  the 
soliton  pulse- widths  considered  there  were  on  the  order  of  fxs  for  mA  amplitudes,  the 
bandwidth  of  the  small  signal  model  is  extremely  narrow  in  comparison  to  that  of  the 
soliton. 

Observations  of  our  circuit  implementation  of  the  diode  ladder  circuit  seem  to 
indicate  that  this  bandwidth  is  too  narrow  to  explain  the  level  of  higher-frequency 
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circuit  noise  present.  This  is  partially  explained  by  the  change  in  the  cutoff  frequency 
as  solitons  are  processed  through  a  node.  That  is,  over  regions  where  the  soliton  com¬ 
ponent  is  significant,  the  equivalent  resistance  of  the  diode  becomes  Req  =  vt/In(t), 
which  is  on  the  order  of  25 Cl  for  a  1mA  soliton.  The  effect  this  has  on  the  linearized 
lattice  is  to  make  the  lattice  effectively  an  all  pass  filter  in  the  vicinity  of  propagating 
solitons. 

As  a  practical  matter,  we  note  that  there  appears  to  be  a  small  amount  of  diode 
leakage  current  present  in  the  circuit  implementation  and  will  explore  the  effect  of 
a  small  bias  current  on  the  dynamics  of  both  the  soliton  components  and  the  small 
amplitude  perturbation.  For  a  solution  containing  a  soliton,  In(t ),  a  small  amplitude 
component,  in(t),  and  a  small  bias  current,  R,  the  resulting  system  dynamics  are 

— 2  In  (l  +  — — =  — {In- 1  —  2 In  +  In+ 1  +  in- 1  —  2 in  +  in+ 1),  (5.39) 

(tv  \  lg  / 

which  reduces  to 

d^ln  (X  +  ljt)  +  l)  =  ~  2ln ®  +  tn+1^-  (5'40) 

When  the  noise  component  in(t)  is  small  as  compared  with  R,  and  away  from  the 
peaks  of  the  soliton  signal,  In(t)  <  R,  the  dynamics  further  reduce  to 

~  ^(»n-l(t)  -  2 in(t)  +  R+l  (*)),  (5.41) 

where,  the  diodes  are  replaced  by  their  linearization  about  the  bias  current,  Req  = 
vt/Ib,  leading  to  an  increase  in  the  bandwidth  of  the  effective  low  pass  filter. 

In  summary,  due  of  the  scaling  of  the  diode  ladder  circuit,  in  order  for  the  linear 
analyses  to  hold,  the  noise  must  be  small  in  comparison  to  the  diode  saturation 
current,  Is.  When  a  soliton  is  present,  if  the  noise  is  small  compared  with  the  soliton, 
then  a  linear  model  can  hold  as  with  the  LC  ladder.  When  there  is  no  soliton,  if  there 
is  a  small  bias  current  that  is  larger  than  the  noise,  this  can  also  lead  to  a  simple 
linear  model.  When  there  is  neither  a  bias  nor  a  soliton  present,  if  the  noise  is  not 
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small  as  compared  with  the  diode  saturation  current,  then  the  noise  satisfies  the  fully 
nonlinear  system.  The  resulting  disturbance  is  better  described  in  terms  of  inverse 
scattering  and  leads  to  the  problem  of  determining  the  spectrum  of  random  linear 
operators,  or  random  matrices. 


5.6  Inverse  Scattering-Based  Noise  Modeling 

The  inverse  scattering  transform  provides  a  particularly  useful  mechanism  for  ex¬ 
ploring  the  long  term  behavior  of  soliton  systems.  In  a  similar  manner  to  the  use 
of  the  Fourier  transform  for  describing  the  ability  of  linear  processors  to  extract  a 
signal  from  a  stationary  random  background,  the  nonlinear  spectrum  of  a  received 
soliton  signal  in  noise  can  effectively  characterize  the  ability  of  the  nonlinear  system 
to  extract  or  process  the  component  solitons.  In  this  section,  we  focus  on  the  ef¬ 
fects  of  random  perturbations  on  the  dynamics  of  solitons  in  the  Toda  lattice  from 
the  viewpoint  of  inverse  scattering.  A  discussion  of  this  topic  for  soliton  systems 
in  general  would  require  a  level  of  mathematical  background  that  is  largely  outside 
the  scope  of  this  thesis.  However,  some  interesting  results  can  be  obtained  for  the 
Toda  lattice  through  an  application  of  some  results  from  matrix  theory,  highlighting 
some  potential  directions  for  future  research  in  the  broad  context  of  soliton  signal 
processing. 

As  seen  in  Sec.  2.2,  the  dynamics  of  the  Toda  lattice  may  be  described  by  the 
evolution  of  the  matrix 


an-2  (t) 

i(t) 

On-l(i) 

m  = 

O'Tl— 1  (t) 

K(t) 

un(t) 

dn(t) 

bn+1  (t)  Q“n- 1-1  (^) 

(5.42) 


whose  eigenvalues  outside  the  range  |A|  <  1  give  rise  to  soliton  behavior.  By  consid¬ 
ering  the  effects  of  small  amplitude  perturbations  to  the  sequences  an(t )  and  bn(t)  on 
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the  eigenvalues  of  L(t),  we  can  observe  the  effects  on  the  soliton  dynamics  through  the 
eigenvalues  corresponding  to  solitons.  There  is  a  large  body  of  literature  on  the  per¬ 
turbation  of  eigenvalues  of  linear  operators  and  matrices,  e.g.  [6,  37,  73,  76].  We  will 
largely  use  results  from  [73]  for  randomly  perturbed  symmetric  matrices  with  simple 
eigenvalues,  since  both  the  matrix  L(t)  and  its  perturbation  will  be  symmetric. 
Following  [73],  we  write  the  N  x  N  matrix  L  in  the  form 

L  =  L0  +  D,  (5.43) 

where  L0  is  the  unperturbed  symmetric  matrix,  and  D  is  the  symmetric  random 
perturbation.  It  is  further  assumed  that  each  element  of  the  perturbation  matrix  is 
almost-surely  bounded,  such  that 

(  N  V/2 

|£>l=lE4J  <e  a-sr  (5-44) 

where  e  >  0  is  sufficiently  small.  By  expanding  the  eigenvalues  of  the  matrix  L  in 
a  convergent  series  and  matching  terms  in  the  expansion,  vomScheidt  and  Purkert 
derive  expressions  for  the  first  four  terms  in  the  expansion, 

OO 

Xg  =  Ms  ~  E  Xgk,  (5.45) 

k= 1 

where  ng  is  the  g-th.  eigenvalue  of  Lq  and  Xgk  is  of  &-th  order  in  the  elements  of  D. 
The  expectations  and  correlations  of  the  eigenvalues  can  also  be  derived  up  to  fourth 
order  in  the  elements  of  D.  To  second  order,  the  eigenvalues  are  given  by 

Affl  =  ~dgg, 

N  d  d 

Xg2  =  E  — ,  (5.46) 

where  Hig  =  fii  —  ng,  and  dig  are  the  elements  of  the  matrix  D  defined  by, 


D  =  CtDC , 


(5.47) 
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and  C  is  a  matrix  that  diagonalizes  L, 


CtL0C  =  diag(/xi, . . . ,  hn)-  (5.48) 

To  second  order,  the  means  of  the  eigenvalues  are  given  by 

£{A„}  =  ft,-  £  (5.49) 

i=l,i£g  rig 

indicating  that  the  eigenvalues  of  L  are  asymptotically  unbiased  estimates  of  the 
eigenvalues  of  L0-  To  first  order, 

Ap  ~  /ip  —  dgg,  (5.50) 


and  dgg  is  a  linear  combination  of  the  elements  of  D, 

N 

dgg  -  ^  ^  CgrCgsdrs.  (5.51) 

r=l,s=l 

Therefore,  if  the  elements  of  D  are  jointly  Gaussian,  then  to  first  order,  the  eigenvalues 
of  L  will  be  jointly  Gaussian,  distributed  about  the  eigenvalues  of  Lq. 

Note  that  the  second  order  perturbation  to  the  eigenvalues  depends  not  only  on 
the  perturbation  matrix,  Z),  but  also  on  the  proximity  to  the  other  eigenvalues.  Also 
apparent  from  (5.49)  is  that  the  smallest  and  largest  eigenvalues  are  biased,  that 
is  £7{Ai}  <  Hi  and  E{ A#}  >  hn-  This  is  of  particular  interest,  since  the  largest 
eigenvalue  typically  corresponds  to  a  soliton. 

The  variance  of  the  eigenvalues  are  in  general  given  by 


(5.52) 


to  third  order.  When  there  is  only  small  amplitude  noise  in  the  processed  signal,  we 
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can  evaluate  (5.52)  to  second  order.  If  the  perturbation  matrix  is  given  by 


an-2 

&n- 2  fin—1  &n— 1 


D  = 


&n— 1  Pn  &n 


®-n  Pn+l  ^n+l 


(5.53) 


where  pn  and  an  are  independent,  identically-distributed  (iid)  random  sequences  with 
variances  and  erg,  then  the  variances  of  the  eigenvalues  reduce  to 


Var(Ap) «  + 


(5.54) 


to  second  order.  This  indicates  that  the  eigenvalues  of  L  are  consistent  estimates  of 
the  eigenvalues  of  Lq. 

Although  we  have  only  given  a  cursory  investigation  of  this  area,  the  inverse  scat¬ 
tering  framework  appears  to  provide  unique  perspective  on  the  effects  of  perturbations 
on  the  long-term  behavior  of  the  soliton  solutions.  If  this  method  is  to  be  applied 
to  the  driven  Toda  lattice  as  discussed  in  Chap.  3,  then  the  sequences  an  and  Pn 
would  no  longer  be  i.i.d.,  although  to  first  order,  they  would  remain  jointly  Gaussian. 
However,  for  the  communication  example  using  the  dKdV  equation,  the  transmitted 
signal  could  correspond  to  the  sequence  an  and  this  analysis  would  apply  directly. 
Although  the  methods  developed  in  [73]  assume  a  bounded  perturbation,  D ,  it  can 
be  shown  that  such  methods  remain  accurate  for  Gaussian  perturbations  of  small 
variance. 

In  this  chapter,  we  have  examined  some  of  the  effects  of  small  amplitude  noise  on 
the  processing  of  soliton  signals  with  the  Toda  lattice.  When  processing  small  ampli¬ 
tude  noise  alone,  the  lattice  can  be  viewed  as  a  dispersive  low  pass  filter.  Therefore, 
the  response  of  the  system  will  be  low  pass  Gaussian  noise  at  each  node.  When 
solitons  are  processed  in  noise,  our  linearized  analyses  indicate  that  the  solitons  will 
be  essentially  unperturbed,  and  the  noise  will  remain  Gaussian  and  low  pass.  These 
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conclusions  are  validated  through  numerical  integration  of  the  fully  nonlinear  system 
and  through  a  state  space  approach.  Through  numerical  solution,  we  see  that  the 
variance  of  the  noise  at  each  node  is  essentially  the  variance  that  would  be  observed 
at  the  output  of  the  low  pass  filter  with  a  time-varying  cutoff  frequency.  Useful 
insight  can  also  be  obtained  from  the  nonlinear  spectrum  of  the  inverse  scattering 
approach.  We  see  that  when  processing  small  amplitude  noise  alone,  the  noise  only 
excites  eigenvalues  corresponding  to  non-soliton  components.  When  solitons  are  also 
processed,  the  noise  induces  a  small  Gaussian  perturbation  to  the  soliton  eigenvalues 
as  well. 
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Chapter  6 


Estimation  of  Soliton  Signals 


In  the  communication  techniques  suggested  in  Chapter  4,  the  parameters  of  a  multi- 
soliton  carrier  are  modulated  with  message-bearing  signals  and  the  carrier  is  then 
processed  with  the  corresponding  nonlinear  evolution  equation.  A  potential  advantage 
to  transmission  of  this  packetized  soliton  carrier  is  a  net  reduction  in  the  transmitted 
signal  energy.  However  during  transmission,  the  multi-soliton  carrier  signal  can  be 
subjected  to  distortions  due  to  propagation,  which  we  have  assumed  can  be  modeled 
as  additive  white  Gaussian  noise,  (AWGN).  In  Chapter  5,  we  characterized  some  of  the 
effects  of  additive  channel  distortion  on  the  generation  and  processing  of  soliton  signals 
using  their  associated  nonlinear  evolution  equations.  In  this  chapter,  we  investigate 
the  fidelity  of  the  demodulated  message  signals  from  the  perspective  of  the  ability  of 
a  receiver  to  estimate  the  parameters  of  a  multi-soliton  carrier. 

Specifically,  if  a  message  is  encoded  in  a  soliton  carrier  by  modulating  the  relative 
spacing  of  the  component  solitons  in  each  period,  then  it  is  important  to  understand 
the  problem  of  time-of-arrival  estimation  both  in  terms  of  algorithm  performance  as 
well  as  theoretical  performance  bounds.  More  generally,  the  efficacy  of  any  system 
employing  soliton  signals  will  be  fundamentally  influenced  by  the  ability  to  detect, 
estimate,  or  otherwise  process  the  signals  themselves.  The  nonlinear  interaction  of 
the  component  waveforms  make  the  signal  dynamics  unique  to  soliton  systems,  and 
it  is  therefore  important  to  begin  to  understand  how  such  signals  behave  under  the 
influence  of  traditional  signal  processing  algorithms. 
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In  this  chapter,  we  consider  the  problems  of  estimating  the  scaling  parameters 
and  the  relative  positions  of  component  solitons  of  multi-soliton  solutions  and  once 
again  focus  on  the  Toda  lattice  as  an  example.  For  each  of  these  problems,  we 
derive  Cramer-Rao  lower  bounds  for  the  estimation  error  variance  through  which 
several  properties  of  multi-soliton  signals  can  be  observed.  Specifically,  we  will  see 
that  although  the  net  transmitted  energy  in  a  multi-soliton  signal  can  be  reduced 
through  nonlinear  interaction,  the  estimation  performance  for  the  parameters  of  the 
component  solitons  can  also  be  enhanced.  However,  at  the  receiver  there  are  inherent 
difficulties  in  parameter  estimation  imposed  by  this  nonlinear  coupling.  We  will 
demonstrate  that  the  Toda  lattice  can  act  as  a  tuned  receiver  for  the  component 
solitons,  naturally  decoupling  them  so  that  the  parameters  of  each  soliton  can  be 
independently  estimated.  Based  on  this  strategy,  we  develop  asymptotically  optimal 
algorithms  for  maximum  likelihood  parameter  estimation.  We  also  extend  the  analogy 
of  the  inverse  scattering  transform  as  an  analog  of  the  Fourier  transform  for  linear 
techniques,  by  developing  a  maximum  likelihood  estimation  algorithm  based  on  the 
nonlinear  spectrum  of  the  received  signal. 

6.1  Soliton  Parameter  Estimation:  Bounds 

In  order  to  understand  the  effects  of  additive  channel  corruption  on  some  of  the 
modulation  techniques  discussed  in  Chap.  4,  we  consider  some  of  the  basic  signal 
processing  issues  that  must  be  addressed  by  any  receiver.  Specifically,  in  this  sec¬ 
tion  we  will  determine  bounds  on  the  performance  of  any  unbiased  estimator  for  the 
problem  of  estimating  both  the  scale  parameters,  /%,  and  the  relative  positions,  of 
multi-soliton  signals  in  stationary  AWGN.  In  our  simplified  channel  model,  the  re¬ 
ceived  signal  r(t)  contains  a  multi-soliton  signal  s(t)  in  an  additive  noise  background 
n(t)  with  noise  power,  N0. 

A  bound  on  the  variance  of  an  estimate  of  the  parameter  f3  may  be  useful  in 
determining  the  demodulation  performance  of  an  AM-like  modulation  or  PAM,  where 
the  component  soliton  wavenumbers  are  slightly  amplitude  modulated  by  a  message- 
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bearing  waveform.  The  fidelity  of  the  recovered  message  will  be  directly  affected  by 
the  ability  of  a  receiver  to  estimate  the  soliton  amplitude,  (3 2 .  When  s(t)  contains  a 
single  soliton  for  the  Toda  lattice, 


s(t)  —  /32sech2  (/?£), 


(6.1) 


the  variance  of  any  unbiased  estimator  of  (3  must  satisfy  the  Cramer-Rao  lower 
bound  (CRB), 


Var  (fl) 


(6.2) 


where  the  observation  interval  is  assumed  to  be  tt  <  t  <  tf.  For  the  infinite  observa¬ 
tion  interval,  —  oo  <  t  <  oo,  the  CRB  (6.2)  is  given  by 


Var 


^  No  _  N0 
-  (|  +  if)/?  ~  3.544/)' 


(6.3) 


Note  that  the  CRB  is  a  decreasing  function  of  (3.  Although  the  width  of  s(t)  increases 
as  (3  — >  0,  the  amplitude  decreases  proportional  to  (32,  and  s(t)  -4  0.  Therefore,  in  a 
fixed  level  of  noise,  estimation  of  (3  will  eventually  become  impossible,  and  the  CRB 
becomes  unbounded.  On  the  other  hand,  as  (3  increases,  the  amplitude  increases  as 
/ 3 2  and  in  the  limit,  (3  can  be  estimated  arbitrarily  well  and  the  CRB  vanishes. 

The  proof  of  the  bound  (6.3)  is  based  on  the  evaluation  of  the  following  integral, 


(6.4) 


For  the  single  soliton  this  becomes, 


r  ( 9/?2sech2(/3t)\2 

L{ — w — j dt 

Expanding  the  integrand  yields, 


(6.5) 


j™  (4/?2sech 3{/3t)  -  8/33sech4(/3t)t  +  4/?4sech4(/?t)  tanh {f3t)t2)  dt.  (6.6) 
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The  first  two  integrands  evaluate  to 


4/?2  (3/?)  ’  and  8/?3  (3/52)  ’ 

respectively.  It  is  the  third  integrand  in  (6.6)  that  causes  difficulty, 


(6.7) 


16  r°°  tdt 
30  J-00  e2t  +  1 


(6.8) 


4(1  +  5e2t  -  5e4t  +  15e6t)t2  -  4(e2t  -  3e4t  -  3e6<  +  e8t)t  +  2(e2t  +  e4t  -  e6t  -  e8t) 


15  (e2t  +  l)b 

(6.9) 

The  second  term  on  the  right  is  easily  seen  to  vanish.  The  integral  (6.9)  evaluates 
to  [22] 

1;5™  =  £  <610> 


A  slightly  different  bound  may  be  useful  in  determining  the  demodulation  perfor¬ 
mance  of  an  FM-like  modulation  or  PPM,  where  the  soliton  position,  or  time-delay, 
is  slightly  modulated  by  a  message-bearing  waveform.  The  fidelity  of  the  recovered 
message  waveform  will  be  directly  affected  by  the  ability  of  a  receiver  to  estimate  the 
soliton  position.  When  the  signal  s(t )  contains  a  single  soliton, 


s(t)  =  /32sech.2(f3(t  —  <5)), 


(6.11) 


where  S  is  the  relative  position  of  the  soliton  in  a  period  of  the  carrier,  the  CRB  for 
8  is  given  by  [72] 


(6.12) 


which  evaluates  to 


Var  (S) 


> 


N0 


N0 


ft!  4/56sech4(/3(t  -  5))  tanh2(/?(t  -  8))dt  (|f )  p ' 


(6.13) 
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As  a  comparison,  for  estimating  the  time  of  arrival  of  the  raised  cosine  pulse, 


s(t)  =  /32  (1  +  cos  (27r/?(t  -  5))) ,  (6-14) 

the  CRB  for  this  more  traditional  pulse  position  modulation  would  be 

Var(i)>^,  (6.15) 

which  has  the  same  dependence  on  signal  amplitude  as  (6.13). 


6.2  Multi-soliton  Parameter  Estimation:  Bounds 


When  the  received  signal  comprises  a  multi-soliton  waveform  where  the  multiple  soli- 
tons  are  not  time  resolved  and  are  therefore  combined  according  to  the  nonlinear 
superposition  of  the  system,  the  estimation  problem  becomes  more  difficult.  Ap¬ 
propriately,  the  bounds  for  estimating  the  parameters  of  such  signals  must  also  be 
sensitive  to  the  relative  positions  of  the  component  solitons.  The  ability  to  estimate 
the  soliton  parameters  based  on  observations  of  the  multi-soliton  received  waveform 
will  give  direct  insight  to  the  potential  performance  of  any  demodulation  algorithm. 

We  will  focus  our  attention  on  the  two-soliton  solution  to  the  Toda  lattice, 

s(t)  _  A2sech2(r?i)  +  ftsech2(7fe)  +  Asech2(T7i)sech2(7?2) 

(cosh(^/2)  +  sinh(<ji/2)  tanh(7?i)  tanh^))2 

where 

A  =  sinh(<£/2)  ((/?2  +  sinh(0/2)  +  2ft &  cosh(<ft/2))  ,  (6.17) 

and  when  both  solitons  travel  in  the  same  direction  and  ft  >  ft, 


4>  —  In 


/  sinh((p!  —  p2)/2)\ 
Vsinh((pa  +p2)l2))  ’ 


(6.18) 


where  ft  =  sinh(pj),  and  rji  =  ft(t  -  ft). 

From  the  modulation  techniques  discussed  in  Chapter  4,  we  are  generally  inter- 
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ested  in  estimating  the  parameters  of  the  multi-soliton  carrier  for  an  unknown  relative 
spacing  among  the  solitons  present  in  the  carrier  signal.  Either  the  relative  spacing 
of  the  solitons  has  been  modulated,  and  is  therefore  unknown,  or  the  parameters  0i 
and  02  are  slightly  modulated  and  the  induced  phase  shift  in  the  received  solitons, 
<f>,  is  unknown.  In  either  case,  the  CRB  for  jointly  estimating  the  parameters  of  a 
multi-soliton  signal  from  observations  of  r(t)  can  be  obtained  numerically  by  forming 
the  Fisher  information  matrix,  7(0),  where  ©  =  [SihPiPzV , 


['(©)k;  = 


1  rlf  ( ds(r ;  ©)  ds{r\ 0)  \ 

Nju  V  de>  d&~) 


dt, 


(6.19) 


where  ©{  is  the  z-th  element  of  0.  The  resulting  bound  on  the  estimation  variance 
for  parameter  i  is  given  by  [72] 


Var(0,)>  [/-‘(e)]...  (6.20) 

For  large  separations,  5  =  Si  —  52,  the  CRB  for  estimating  the  parameters  of  any  of 
the  component  solitons  will  be  unaffected  by  the  parameters  of  the  other  solitons.  As 
shown  in  Fig.  6-1,  when  the  component  solitons  are  well  separated,  the  CRB  for  either 
Pi  or  02  approaches  the  CRB  for  estimation  of  a  single  soliton  with  that  parameter 
value  in  the  same  level  of  noise.  The  bounds  for  estimating  0i,  and  02  are  shown  in 
Fig.  6-1  as  a  function  of  the  relative  separation,  <5. 

Note  that  both  of  the  bounds  are  reduced  by  the  nonlinear  superposition,  indi¬ 
cating  that  the  potential  performance  of  the  receiver  is  enhanced  by  the  nonlinear 
superposition.  However,  if  we  let  the  parameter  difference  02  —  Pi  increase,  we  notice 
a  different  character  to  the  bounds.  Specifically,  we  maintain  0i  —  sinh(2),  and  let 
0i  =  sinh(1.25)  and  plot  in  Fig.  6-2,  the  same  bounds  as  Fig.  6-1. 

Note  that  now,  the  performance  of  the  larger  soliton  is  inhibited  by  the  nonlinear 
superposition,  while  the  smaller  soliton  is  still  enhanced.  In  fact,  the  CRB  for  the 
smaller  soliton  becomes  lower  than  that  for  the  larger  soliton  near  the  range  5  =  0. 
The  reason  for  this  has  to  do  with  the  relative  sensitivity  of  the  signal  s(t )  to  each  of 
the  parameters  0i  and  02- 


108 


Figure  6-1:  The  Cramer- Rao  lower  bound  for  estimating  81  —  sinh(2)  and  82  = 
sinh(1.75)  with  all  parameters  unknown  in  AWGN  with  N0  =  1.  The  bounds  are 
shown  as  a  function  of  the  relative  separation,  5  =  —  S2.  The  CRB  for  estimating 

fli  and  82  of  a  single  soliton  with  the  same  parameter  value  is  indicated  with  ‘o’  and 
‘x’  marks,  respectively. 


The  ability  to  simultaneously  enhance  estimation  performance  while  decreasing 
signal  energy  is  an  inherently  nonlinear  phenomena.  From  Eq.  (6.2),  we  see  that  the 
CRB  for  linear  estimation  in  AWGN  is  inversely  proportional  to  the  signal  energy, 


Var  (/?) 


> 


Np 

ft!  s2(t)dt 


SNR"1. 


(6.21) 


However,  for  nonlinear  parameter  estimation,  the  CRB  is  dependent  upon  the  energy 
in  the  derivative  of  the  signal  with  respect  to  the  parameter,  rather  than  signal  energy. 

The  bounds  for  estimating  the  times  of  arrival  of  the  two  component  solitons  can 
also  be  seen  to  agree  with  the  single  soliton  bounds,  as  shown  in  Fig.  6-3.  Note  that 
while  the  estimation  performance  of  the  larger  of  the  two  solitons  is  inhibited  by 
nonlinear  superposition,  the  estimation  of  the  smaller  of  the  two  solitons  is  enhanced. 

For  the  two-soliton  signal,  the  nonlinear  interaction  of  the  component  solitons  can 
therefore  enhance  the  parameter  estimation  performance  or  permit  a  tradeoff  in  per¬ 
formance  between  the  component  solitons.  This  potential  performance  enhancement 
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Figure  6-2:  The  Cramer-Rao  lower  bound  for  estimating  Pi  =  sinh(2)  and  P2  = 
sinh(1.25)  with  all  parameters  unknown  in  AWGN  with  Nq  =  1.  The  bounds  are 
shown  as  a  function  of  the  relative  separation,  6  —  Si  —  62-  The  CRB  for  estimating 
Pi  and  P2  of  a  single  soliton  with  the  same  parameter  value  is  indicated  with  ‘o’  and 
‘x’  marks,  respectively. 

is  achieved  in  combination  with  a  net  energy  reduction  in  the  multi-soliton  signal. 
The  combination  of  these  properties  may  make  the  superimposed  solitons  particularly 
attractive  for  a  variety  of  communications  contexts  by  providing  an  overall  signal  to 
noise  ratio  enhancement. 


6.3  Estimation  Algorithms 

6.3.1  General  Approach 

In  this  section  we  will  present  and  analyze  the  performance  of  several  algorithms  for 
estimating  the  parameters  of  soliton  signals.  Although  the  algorithms  presented  do 
not  amount  to  a  demodulation  of  any  of  the  communication  techniques  suggested 
in  Chap.  4,  the  bounds  presented  in  the  previous  sections  should  give  an  idea  of 
the  best  achievable  performance  of  any  such  demodulation  hardware.  Similarly,  the 
performance  of  the  algorithms  presented  here  may  indicate  the  degree  to  which  such 
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Figure  6-3:  The  Cramer-Rao  bounds  for  estimating  the  time  of  arrival  for  each  soliton 
in  a  two-soliton  signal  with  (3\  =  sinh(2)  and  /?2  =  sinh(1.25)  in  AWGN  with  N0  =  1. 
The  asymptotic  values  of  each  of  the  bounds  agree  with  the  CRB  for  estimating  the 
time  of  arrival  of  a  single  soliton  with  the  same  parameter  value  as  indicated  with  ‘o’ 
and  ‘x’  marks. 

bounds  may  be  approached. 

As  an  example,  we  will  focus  on  the  diode  ladder  circuit  implementation  of  the 
Toda  lattice  equations, 

cP 

—  ln(l  +  in(t))  =  -  2 in(t)  +  (6.22) 

where  in(t)  is  the  current  through  the  n-th  diode,  z0 (f)  =  and  for  simplicity, 

the  parameters  of  all  circuit  elements  have  been  normalized  to  unity. 

In  order  to  motivate  the  general  approach,  we  start  by  considering  the  problem  of 
estimating  the  position,  5,  of  a  single  soliton  solution  to  (6.22), 

s(t ;  8)  =  /?2sech2(/?(i  -  (5)),  (6.23) 

with  the  parameter  (3  known.  This  is  a  classical  time-of-arrival  estimation  problem. 
For  observations  r(t)  =  s(t)  +  n(t),  where  n(t)  is  a  stationary  white  Gaussian  process, 
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the  maximum  likelihood  estimate  is  given  by  the  value  of  the  parameter  8  which 
maximizes  the  likelihood  of  the  observation,  or  equivalently  minimizes  the  exponent 
of  the  likelihood  function  yielding  the  expression 

5  =  argmin  f  f  (r(t)  -  s(t;r))2dt,  (6.24) 

Jti 

where  s(t;r)  is  the  set  of  “replica  signals”  corresponding  to  s(t  —  r)  for  all  possible 
values  of  the  parameter  <5.  Since  the  replica  signals  all  have  the  same  energy,  we  can 
represent  the  minimization  in  (6.24)  as  a  maximization  of  the  correlation, 

~  /■*/ 

8  =  argmin/  r(t)s(t  —  r)dt.  (6.25) 

T  Jti 

It  is  well-known  that  an  efficient  way  to  perform  the  correlation  (6.25)  with  all  of  the 
replica  signals  s(t  —  r)  over  the  range  (5mill  <  r  <  8max,  is  through  convolution  with 
a  matched  filter  followed  by  a  peak-detector  [72].  The  estimate  8  is  then  determined 
from  the  time-location  of  the  peak. 

When  the  signal  r(t )  contains  a  multi-soliton  signal,  s(t;{3,8 ),  where  we  wish  to 
estimate  the  parameter  vector  8,  the  estimation  problem  becomes  more  involved. 
If  the  component  solitons  are  well  separated  in  time,  then  the  maximum  likelihood 
estimator  for  the  positions  of  each  of  the  component  solitons  would  again  involve  a 
matched  filter  processor  followed  by  a  peak-detector  for  each  soliton. 

If  the  component  solitons  are  not  well-separated  and  are  therefore  nonlinearly 
combined,  the  estimation  problems  are  tightly  coupled  and  should  not  be  performed 
independently.  For  observations  in  AWGN,  the  maximum  likelihood  processor  would 
involve  a  minimization  of  the  energy  in  the  difference  between  the  observed  signal 
r(t )  and  a  replica  signal  s(t;[3,8)  over  the  parameter  space, 

a  rtf 

8  =  argmin  /  (r(t)  -  s(t] /3,0))2dt.  (6.26) 

2.  Jt{ 

Note  that,  since  the  set  of  signals  s(t:  fl,  6)  are  not  all  of  equal  energy,  the  minimization 
cannot  be  represented  as  a  correlation.  Further,  since  the  signals  are  not  all  related 
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by  a  time-shift,  there  is  no  simple  way  to  reproduce  each  member  of  the  set. 

The  estimation  problems  can  be  decoupled  by  preprocessing  the  signal  r(t)  with 
the  Toda  lattice.  By  setting  im(t)  =  r(t),  that  is  the  current  through  the  first  diode 
in  the  diode  ladder  circuit,  then  as  the  signal  propagates  through  the  lattice,  the 
component  solitons  will  naturally  separate. 

Defining  the  signal  and  noise  components  as  viewed  on  the  k- th  node  in  the  lattice 
as  ik(t )  by  sk(t)  and  nk(t),  respectively,  i.e., 

h(t)  =  sk(t)  +  nk(t),  (6.27) 

where  n0(t)  is  the  stationary  white  Gaussian  noise  process  n(t),  in  Chapter  5,  we  saw 
that  in  the  high  SNR  limit,  nk(t )  will  be  low  pass  and  Gaussian.  In  this  limit,  the 
ML  estimator  for  the  positions,  5,,  can  again  be  formulated  using  matched  filters  for 
each  of  the  component  solitons.  Since  the  lattice  equations  are  invertible,  at  least  in 
principle  through  inverse  scattering,  then  the  ML  estimate  of  the  parameter  5  based 
on  r(t)  must  be  the  same  as  the  estimate  based  on  in(t)  =  T(r(f)),  for  any  invertible 
transformation  T(-).  If  the  component  solitons  are  well-separated  as  viewed  on  the 
IV-th  node  of  the  lattice,  ijv(f),  then  an  ML  estimate  based  on  observations  of  ijv(£) 
will  reduce  to  the  aggregate  of  ML  estimates  for  each  of  the  separated  component 
solitons  in  low  pass  Gaussian  noise.  For  soliton  position  estimation,  this  amounts 
to  a  bank  of  matched  filters.  We  can  view  this  estimation  procedure  as  a  form  of 
nonlinear  matched  filtering,  whereby  first,  dynamics  matched  to  the  soliton  signals 
are  used  to  perform  the  necessary  signal  separation,  and  then  filters  matched  to  the 
separated  signals  are  used  to  estimate  their  arrival  time. 

To  examine  the  empirical  performance  of  this  approach,  we  consider  a  soliton 
signal  in  AWGN  with  noise  power  Nq.  Since  the  bandwidth  limitations  of  the  receiver 
will  restrict  the  possible  range  of  parameters,  and  to  simplify  our  simulations,  we 
assume  that  the  receiver  comprises  a  low  pass  filter  followed  by  a  Toda  lattice  circuit 
as  shown  in  Fig.  6-4.  We  also  assume  that  the  bandwidth,  2n/A,  of  the  low  pass  filter 
in  Fig.  6-4,  is  wide  enough  to  pass  the  soliton  components  of  r(t )  completely.  The 
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Figure  6-4:  Toda  lattice  receiver  model. 


input  to  the  Toda  lattice  circuit,  iin(f),  then  comprises  the  soliton  signal  in  low  pass 
Gaussian  noise.  Simulations  are  performed  using  a  Runge-Kutta  integration  routine 
with  a  fixed  step  size,  A.  To  model  the  effects  of  the  noise,  an  i.i.d.  Gaussian  random 
sequence,  w(kA)  ~  iV(0,  a£),  is  added  to  the  samples  of  the  input  sequence  iin(kA) 
resulting  in  an  effective  white  noise  power  of  No  =  A a2.  Equivalently,  to  simulate 
the  effect  of  a  white  noise  component  in  r(t )  with  noise  power  N0,  and  a  bandwidth 
of  the  low  pass  filter  of  2ojc,  the  Runge-Kutta  time  step  should  be  set  to  A  =  n/uc 
and  the  variance  of  the  Gaussian  i.i.d.  sequence  should  be  a 2  =  N0/A. 


6.3.2  Position  Estimation 

We  will  focus  our  attention  on  the  two-soliton  signal  (6.16).  If  the  component  solitons 
are  well-separated  as  viewed  on  the  N -th  node  of  the  Toda  lattice,  the  signal  appears 
to  be  a  linear  superposition  of  two  solitons, 

iN(t)  pa  /3lsech2(/3i(t  —  Si)  - piN  -  <j)/2)  (6.28) 

+  /?fsech2(/?2(t  -  fa)  -  PiN  +  <j>/ 2), 


where  (j>/2  is  the  time-shift  incurred  due  to  the  nonlinear  interaction.  Matched  filters 
can  now  be  used  to  estimate  the  time  of  the  arrival  of  each  soliton  at  the  iV-th  node. 
We  formulate  the  estimate 


PiN  +  <j>/2\ 

A  ) 


p2N  -  <j>/2\ 
"fa  )’ 


(6.29) 


where  taN  l  is  the  time  of  arrival  of  the  i-th  soliton  an  node  N.  The  performance 
of  this  algorithm  for  a  two-soliton  signal  with  /?  =  [sinh(2),sinh(1.5)]  is  shown  in 
Fig.  6-5.  Note  that  although  the  error  variance  of  each  estimate  appears  to  be  a 
constant  multiple  of  the  CRB,  the  estimation  error  variance  approaches  the  CRB  in 
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Figure  6-5:  The  CRJBs  for  ft  and  <52  are  shown  with  solid  and  dashed  lines,  while 
the  estimation  error  results  of  100  Monte-Carlo  trials  are  indicated  with  ‘o’  and  ‘  x  ’ 
marks,  respectively. 

an  absolute  sense  as  jV0  — i ►  0. 


6.3.3  Velocity  Estimation 

Since  the  velocity  of  a  soliton  in  the  Toda  lattice  is  related  to  the  scale  parameter,  ft 
and  the  velocity  can  be  estimated  based  on  the  arrival  times  of  the  soliton,  we  may 
also  form  an  estimate  of  through  velocity  estimation.  From  the  set  of  arrival  times, 
taN  i  we  have  a  set  of  equations  for  each  of  the  unknown  parameters,  ft, 


PiN  +  0/2  =  ^t\N,  p2N  -  0/2  =  (32tlN,  (6.30) 


which  are  nonlinear  in  the  parameter  ft,  since  p*  =  sinh  x(ft). 
parameter  p*,  we  obtain, 


sinh  (pj) 
Pi 


Solving  for  the 
(6.31) 


where  tf  N  are  the  time-difference  of  arrivals  of  the  i-th  soliton  between  nodes  N  and 
N  +  1.  Equation  (6.31)  can  be  solved  numerically  using  Newton’s  method  and  for 
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large  p,  the  solution  can  also  be  expressed  in  terms  of  Lambert’s  omega  function,  pt  = 
— fi(— tf  N)  [9].  We  may  then  formulate  the  estimate  /%  =  sinh(p,).  This  technique  is 
similar  to  “velocity-filtering”  commonly  used  in  seismic  signal  processing  [55].  Note 
that  in  estimating  /%,  the  arrival  times  taN  -  must  be  selected  based  on  peak  detection 
of  the  signal  alone  without  the  aid  of  the  matched  filter  which  would  require 
knowledge  of  the  unknown  parameter  pi. 

A  scalar  closed-form  sufficient  statistic  for  the  estimation  of  p  based  on  observa¬ 
tions  of  a  single  soliton  in  Gaussian  noise  does  not  exist  due  to  the  highly  nonlinear 
manner  in  which  the  parameter  appears  both  in  the  time  and  the  amplitude  scales  of 
the  signal  which  preclude  factorization  of  the  likelihood  function  for  r(t).  Although 
an  ML  solution  could  be  found  by  numerical  maximization  of  the  likelihood  function, 
such  a  solution  would  be  computationally  intensive.  We  therefore  look  at  the  perfor¬ 
mance  of  the  ad-hoc  estimate  based  on  the  most  salient  characteristics  of  the  soliton 
through  the  lattice,  namely  the  scale-dependent  velocity. 

The  results  of  100  Monte  Carlo  trials  for  the  velocity-based  estimation  algorithm 
using  a  Toda  lattice  with  N  =  10  nodes  are  shown  in  Fig.  6-6.  The  parameter  values 
estimated  were  Pi  =  sinh(2)  and  p2  =  sinh(1.5)  with  both  solitons  time-aligned. 
Therefore  the  CRB  for  the  smaller  soliton  is  actually  lower  than  that  for  the  larger 
soliton. 

6.3.4  Estimation  Based  on  Inverse  Scattering 

Although  an  ML  estimate  of  the  parameter  P  cannot  be  formulated  directly  based 
on  observations  of  the  received  signal  r(t),  we  can  use  the  inverse  scattering  frame¬ 
work  along  with  some  of  the  results  from  Chapter  5  to  demonstrate  a  procedure  for 
constructing  the  ML  estimate  in  the  high  SNR  limit.  We  have  argued  that  ML  pa¬ 
rameter  estimation  based  on  a  received  multi-soliton  signal,  r(t),  is  equivalent  to  ML 
estimation  based  on  the  signal  as  observed  on  the  iV-th  node  in  the  Toda  lattice, 
ijv(t).  We  may  also  consider  the  transformation 

L(t)  =  T{r(t)},  (6.32) 
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Figure  6-6:  The  estimation  error  variance  for  the  velocity-based  algorithm.  The  CRB 
for  pi  and  02  is  shown  with  a  solid  and  dashed  line,  and  the  estimation  error  variances 
are  indicated  by  the  points  labeled  ‘o’  and  ‘  x  ’  respectively. 

where  L{t)  is  the  symmetric  matrix  from  the  inverse  scattering  transform.  This  trans¬ 
formation  T{-}  is  also  invertible  in  principle  through  inverse  scattering.  Therefore,  an 
ML  estimate  based  on  the  signal  L(t)  must  be  the  same  as  an  ML  estimate  based  on 
r(t).  We  therefore  seek  to  form  an  estimate  of  the  parameters  of  the  signal  r(t)  by  per¬ 
forming  the  estimation  in  the  nonlinear  spectral  domain.  This  can  be  accomplished 
as  shown  in  Fig.  6-7  where  the  Toda  lattice  can  be  viewed  as  a  nonlinear  filter  bank 
which  projects  the  signal  r(t)  onto  the  spectral  components  of  L(t).  This  use  of  the 
inverse  scattering  transform  is  analogous  to  performing  frequency  estimation  with  the 
Fourier  transform. 

From  Chapter  2  we  saw  that  if  vn(t )  evolves  according  to  the  Toda  lattice, 


-(Vn(t)-Vn-l(t)) 


e~(Vn+l(t)-Vn(t))\ 


(6.33) 
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Figure  6-7:  Procedure  for  computing  L(t)  by  processing  the  signal  r(t)  with  the  Toda 
lattice. 


then  the  eigenvalues  of  the  matrix, 


<^n—2  (t)  bn—l(.t) 

an-i(t) 

m  = 

an-i{t) 

bn(t)  On(^) 

®n(^)  ^n+lW  ®n+l(^) 

(6.34) 


are  time-invariant,  where,  an(t )  =  ^Vn^~Vn+1^)!2 ,  and  bn  =  vn(t)/2.  Note  that  vn(t) 
could  correspond  to  the  voltages  on  the  double  capacitors  in  the  diode  ladder  circuit. 
Further,  the  eigenvalues  of  L(t)  for  which  |A,|  >  1  correspond  to  soliton  solutions, 
with  ^  —  sinh(cosh-1(Aj))  =  J Af  —  1.  We  also  saw  in  Chapter  5,  that  the  eigenvalues 
of  L(t)  are,  to  first  order,  jointly  Gaussian  and  distributed  about  the  true  eigenvalues 
corresponding  to  the  original  multi-soliton  signal,  s(t).  Therefore,  estimation  of  the 
parameters  Pi  from  the  eigenvalues  of  L(t)  as  described  above  constitutes  a  maximum 
likelihood  approach  in  the  high  SNR  limit. 

The  parameter  estimation  algorithm  now  amounts  to  an  estimation  of  the  eigenval¬ 
ues  of  L(t).  Note  that  since  L(t)  is  tridiagonal,  very  efficient  techniques  for  eigenvalue 
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estimation  may  be  used  [21].  The  estimate  of  the  parameter  j3  is  then  found  by  the 
relation  fa  =  yjXj  —  1,  where  j A*  ]  >  1,  and  the  sign  of  fa  can  be  recovered  from  the 
sign  of  A;.  Clearly  if  there  is  a  pre-specified  number  of  solitons,  k,  present  in  the 
signal,  then  the  k  largest  eigenvalues  would  be  used  for  the  estimation.  If  the  number 
k  were  unknown,  then  a  simultaneous  detection  and  estimation  algorithm  would  be 
required. 

In  order  to  perform  eigenvalue  estimation,  we  must  once  again  deal  with  the  finite 
length  of  the  Toda  lattice.  This  can  be  resolved  by  either  using  the  periodic  Jacoby 
matrix  L(t)  that  results  from  making  the  periodic  assumption  a^(t)  —  a0(t), 


h(t)  afat)  0  •••  aN{t) 

afat)  b2(t)  a2(t ) 


L  Ojv(t) 


0 


aN-2{t ) 
aN-2(t)  bw-fat) 

ajv-i(t)  bff(t) 


(6.35) 


or  by  simply  truncating  the  matrix.  It  can  be  shown  that  for  the  periodic  Toda 
lattice,  the  eigenvalues  of  the  periodic  Jacoby  matrix  are  time-invariant  [68].  In 
practice,  truncating  the  matrix  around  the  region  of  support  of  the  solitons  has  little 
effect  on  the  eigenvalues.  Note  the  similarity  once  again  to  the  frequency  estimation 
problem,  whereby  the  Fourier  transform  cannot  be  computed  explicitly,  however  the 
discrete  Fourier  Transform,  (DFT),  is  often  used  which  makes  similar  assumptions 
about  the  region  of  support  of  the  underlying  signal. 

An  example  of  the  joint  estimation  of  the  parameters  of  a  two-soliton  signal  is 
shown  in  Fig.  6-8.  The  estimation  error  variance  decreases  with  the  noise  power 
at  the  same  exponential  rate  as  the  CRB.  The  normalized  mean  estimation  error, 
(fa  —  fa)/ fa  for  each  estimate  fa  are  shown  in  Fig.  6-9,  indicating  that  the  estimates 
are  empirically  unbiased. 

To  verify  that  the  performance  of  the  estimation  algorithm  has  the  same  depen¬ 
dence  on  the  relative  separation  of  solitons  as  indicated  in  Sec  6.2,  the  estimation 
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Figure  6-8:  The  estimation  error  variance  for  the  inverse  scattering-based  estimates 
of  (3i  —  sinh(2).  /?2  =  sinh(1.5).  The  bounds  for  and  fe  are  indicated  with  solid 
and  dashed  lines  respectively.  The  estimation  results  for  100  Monte  Carlo  trials  with 
a  diode  lattice  of  N  =  10  nodes  for  and  ft2  are  indicated  by  the  points  labeled  ‘o’ 
and  ‘x’  respectively. 

error  variance  is  also  indicated  in  Fig.  6-10  versus  the  relative  separation,  <5.  In  the 
figure,  the  mean-squared  parameter  estimation  error  for  each  of  the  parameters  /% 
are  shown  along  with  their  corresponding  CRB.  At  least  empirically,  we  see  that  the 
fidelity  of  the  parameter  estimates  are  indeed  enhanced  by  their  nonlinear  interac¬ 
tion,  even  though  this  corresponds  to  a  signal  with  lower  energy,  and  therefore  lower 
observational  SNR. 


6.4  Summary 

In  this  chapter,  we  have  explored  some  of  the  properties  of  multi-soliton  signals  from 
the  viewpoint  of  parameter  estimation.  In  fact,  some  of  the  rich  structure  of  these 
signals  is  apparent  from  the  behavior  of  the  Cramer-Rao  bounds  for  estimating  their 
parameters.  In  Chap.  4,  we  observed  that  the  energy  in  a  two-soliton  signal  is  re¬ 
duced  when  the  component  solitons  are  nonlinearly  superimposed.  In  this  chapter, 
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Figure  6-9:  Normalized  mean  parameter  estimation  error  for  the  inverse  scattering- 
based  estimation. 

we  observed  that  the  nonlinear  interaction  of  solitons  can  also  improve  the  estimation 
performance  of  each.  For  certain  parameter  ranges,  the  coupling  enables  the  larger 
soliton  components  to  enhance  the  estimation  performance  of  the  smaller  components. 

For  multi-soliton  parameter  estimation,  we  noted  that  the  coupling  of  component 
solitons  can  make  traditional  techniques  prohibitively  complex.  However,  using  the 
nonlinear  evolution  equations  as  a  preprocessor  to  perform  signal  separation,  we  were 
able  to  perform  maximum  likelihood  time-delay  estimation  and  an  associated  velocity 
estimation. 

We  also  demonstrated  an  algorithm  for  multi-soliton  parameter  estimation  based 
on  the  inverse  scattering  transform.  In  the  high  SNR  limit,  the  discrete  eigenvalues 
in  the  nonlinear  spectrum  can  be  used  to  formulate  maximum  likelihood  estimates 
of  the  scaling  parameters  for  Toda  lattice  solitons.  This  approach  is  analogous  to  a 
filterbank  which  measures  the  content  of  the  received  signal  in  the  nonlinear  spectrum 
of  the  associated  linear  operator.  Parameter  estimation  based  on  inverse  scattering 
techniques  also  has  many  direct  analogies  to  estimation  based  on  Fourier  methods  for 
linear  systems. 
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Figure  6-10:  The  estimation  error  variance  for  estimating  /?i  =  sinh(2),  #2  =  sinh(1.5) 
are  indicated  with  the  points  labeled  ‘o’  and  ‘  x  ’  and  the  CRBs  for  each  are  indicated 
with  solid  and  dashed  lines,  respectively. 

Perhaps  there  is  something  to  be  gained  from  a  more  thorough  investigation  of  the 
inverse  method  in  looking  for  a  sufficient  statistic  for  soliton  parameter  estimation. 
The  lack  of  a  sufficient  statistic  for  estimation  of  the  parameter  /?  stems  from  the 
highly  nonlinear  manner  in  which  the  parameter  appears  both  in  the  time-scale  and 
the  amplitude  scale  of  the  signal 

sn(t;  (3)  =  /?2sech2(sinh-1(/?)n  —  jit).  (6.36) 

An  investigation  into  exploiting  the  scale-invariance  of  such  signals  and  the  potential 
for  multi-scale  processing  of  soliton  signals  also  represents  an  interesting  area  for 
future  work. 
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Chapter  7 


Detection  of  Soliton  Signals 

7.1  General  Approach 

The  problem  of  detecting  a  single  soliton  or  multiple  non-overlapping  solitons  in 
AWGN  falls  within  the  theory  of  classical  detection.  The  general  approach  to  such 
problems  involves  first  mapping  the  observation  waveform,  r(t)  onto  a  convenient 
and  suitably  compact  decision  space  [72].  The  detection  problem  then  reduces  to 
a  partitioning  of  the  decision  space  into  decision  regions  according  to  a  particular 
criteria  of  optimality.  Consider  the  following  binary  hypothesis  problem, 

H0  :  r(t)  =  n(t), 

Hi  :  r(t)  =  VEs(t)  +  n(t), 

where  the  received  waveform  comprises  the  soliton  signal  y/Es^t)  in  white  Gaus¬ 
sian  observational  noise  n(t)  under  hypothesis  1,  or  just  the  observational  noise  under 
the  null  hypothesis,  for  later  convenience  we  have  normalized  s(t)  to  have  unit  en¬ 
ergy.  The  first  task,  reduction  of  r(t )  into  a  convenient  decision  space,  is  accomplished 
through  the  sufficient  statistic, 

r  —  f  f  r(t)s(t),  (7.1) 

Jti 
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which  can  be  efficiently  computed  with  a  matched  filter. 

Once  the  received  signal  has  been  transformed  into  the  decision  space,  we  have 
the  simple  binary  hypothesis  test 


H0  :  r  =  n, 

Hi  :  r  =  VE  +  n, 


where  under  H0  :  r  ~  7V(0,  No),  and  under  H i  :  r  ~  N(\/E,  jV0),  where  N0  is  the 
noise  power  of  the  white  noise  process,  n(t).  The  optimal  receiver  for  a  variety  of 
measures,  including  the  Bayes  or  Neyman-Pearson  criteria,  involves  a  likelihood  ratio 
test,  where  the  likelihood  function  is  given  by 


AM  =  Pr\Hi(R\Hi) 

[)  prlHo(R\H0y 


(7.2) 


resulting  in  the  test 


Hi 


f‘ r(t)s(t)dt  >  1 n»+ (7.3) 
Ho 

When  the  signal  r(t )  contains  a  multi-soliton  signal  where  the  component  solitons 
are  not  resolved,  the  detection  problem  becomes  more  involved.  Specifically,  consider 
a  signal  comprising  a  two-soliton  solution  to  the  Toda  lattice,  where  we  wish  to  make 
one  of  four  decisions, 


H0  : 

neither  soliton  is  present, 

Hi  : 

only  soliton  one  is  present, 

H2  •' 

only  soliton  two  is  present, 

Hu  : 

both  solitons  are  present. 

If  the  relative  positions  of  the  component  solitons  are  known  a  priori,  then  the  de- 
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tection  problem  reduces  to  deciding  which  among  4  possible  signals  is  present, 


Ho  : 

:  r(t)  =  n(t), 

Hi  : 

:  r(t)  =  Si(t)  -f-n(i), 

H2  : 

:  r(t)  =  s2{t) +n(t), 

Hu  : 

■  r(t)  -  Si2{t)  4- 

n{t), 

where  si(t),s2(t),  and  Si2(£)  are  soliton  one,  soliton  two  and  the  multi-soliton  sig¬ 
nals  respectively.  Once  again,  this  problem  can  be  solved  with  standard  Gaussian 
detection  theory. 

If  the  relative  positions  of  the  solitons  are  unknown,  as  would  be  the  case  for  a 
modulated  soliton  carrier,  then  the  signal  s  12(f)  will  vary  significantly  as  a  function 
of  the  relative  separation.  Similarly,  if  the  signals  are  to  be  transmitted  over  a  soliton 
channel  where  different  users  occupy  adjacent  soliton  wavenumbers,  any  detection  at 
the  receiver  would  have  to  be  performed  with  the  possibility  of  another  soliton  com¬ 
ponent  present  at  an  unknown  position.  We  therefore  obtain  a  composite  hypothesis 
testing  problem,  whereby  under  each  hypothesis,  we  have 


H0  :  r(t)  =  n(t), 

Hi  :  r(t)  =  Si(t;  5X)  4-  n(f), 
H2  :  r(t)  =  s2(t‘82) +n(t), 
Hi 2  :  r(t)  =  Si2(t;S)  +  n(t), 


where  8  =  [61,  <52]T-  The  general  problem  of  detection  with  an  unknown  parameter,  8, 
can  be  handled  in  a  number  of  ways.  For  example,  if  the  parameter  can  be  modeled 
as  random  and  the  distribution  for  the  parameter  were  known,  pa(£),  along  with  the 
distributions  pr\s,H{R\8,  Hi)  for  each  hypothesis,  then  the  Bayes  or  Neyman-Pearson 
criteria  can  be  used  to  formulate  a  likelihood  ratio  test  based  on  the  likelihood  func- 
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tions, 


*  /  \  =  I  dSps_(S)Pr\s,H(R\^Hi) 

U  SdSPi(®PrHH(R\&,HS)  1  j 

for  i  7^  j.  Unfortunately,  even  when  the  distribution  for  the  parameter  5  is  known, 
the  likelihood  ratios  given  in  (7.4)  cannot  be  found  in  closed  form  for  even  the  single 
soliton  detection  problem. 

Another  approach  that  is  commonly  used  in  radar  processing  [72,  24]  applies 
when  the  distribution  of  S  does  not  vary  rapidly  over  range  of  possible  values  while 
the  likelihood  function  (7.4)  has  a  sharp  peak  as  a  function  of  S.  In  this  case,  the 
major  contribution  to  the  integral  in  the  averaged  likelihood  function  (7.4)  is  due  to 
the  region  around  the  value  of  <5  for  which  A(r;  S)  is  maximum,  and  therefore  this 
value  of  the  likelihood  function  is  used  as  if  the  maximizing  value,  Jml,  were  the 
actual  value.  Since  this  amounts  to  a  generalization  of  the  standard  likelihood  ratio 
test,  it  is  often  called  a  “generalized  likelihood  ratio  test”  (GLRT). 

If  we  plan  to  employ  a  GLRT  for  the  multi-soliton  detection  problem,  we  are  again 
faced  with  the  need  for  an  ML  estimate  of  the  position,  |ML.  A  standard  approach  to 
such  problems  would  involve  turning  the  current  problem  into  one  with  hypotheses 
Ho,  Hi,  and  H2  as  before,  and  an  additional  M  hypotheses  —  one  for  each  value  of  the 
parameter  5  sampled  over  a  range  of  possible  values.  Additionally,  the  complexity  of 
the  detection  problem  increases  exponentially  with  the  number  of  component  solitons, 
Ns,  resulting  in  a  hypothesis  testing  problem  with  (M  +  1)^  hypotheses. 

However,  as  with  the  estimation  problems  in  Chap.  6,  the  detection  problems  can 
be  decoupled  by  preprocessing  the  signal  r(t)  with  the  Toda  lattice.  If  the  component 
solitons  separate  as  viewed  on  the  N-th  node  in  the  lattice,  then  the  detection  problem 
can  be  more  simply  formulated  using  ijv(t).  The  invertibility  of  the  lattice  equations 
implies  that  a  Bayes  optimal  decision  based  on  r(t)  must  be  the  same  as  that  based 
on  zjy(t).  Since  the  Bayes  optimal  decision  can  be  performed  based  on  the  likelihood 
function  A(r(t)),  and  A(i^(t))  =  A (T{r(t)j)  =  A (r(t)),  the  optimal  decisions  based 
on  r(t)  and  ijv(t)  must  be  the  same  for  any  invertible  transformation  T{-}.  Although 
we  will  be  using  a  GLRT,  where  the  value  of  SML  is  used  for  the  unknown  positions 
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of  the  multi-soliton  signal,  since  the  ML  estimates  based  on  r(t)  and  must  also 
be  the  same,  as  shown  in  Chap.  6,  the  detection  performance  of  a  GLRT  using  those 
estimates  must  also  be  the  same.  Since  at  high  SNR,  the  noise  component  of  the 
signal  iN{t)  can  be  assumed  low  pass  and  Gaussian,  as  in  Chap.  6,  the  GLRT  can  be 
performed  by  pre-processing  r(t)  with  the  Toda  lattice  equations  followed  by  matched 
filter  processing. 

7.2  Simulations 

To  illustrate  the  algorithm,  we  consider  the  hypothesis  test  between  Ho  and  Hu, 
where  the  separation  of  the  two  solitons,  —  62,  varies  randomly  in  the  interval 
[— 1//?2,  l/ft]-  The  detection  processor  comprises  a  Toda  lattice  of  N  =  20  nodes, 
with  the  detection  performed  based  on  the  signal  iio(t).  To  implement  the  GLRT, 
we  search  over  a  fixed  time  interval  about  the  expected  arrival  time  for  each  soliton. 
In  this  manner  we  obtain  a  sequence  of  1000  Monte  Carlo  values  of  the  processor 
output  for  each  hypothesis.  A  set  of  Monte  Carlo  runs  has  been  completed  for  each 
of  3  different  levels  of  the  noise  power,  N0. 

The  receiver  operating  characteristic  (ROC)  for  the  soliton  with  /32  =  sinh(1.5)  is 
shown  in  Fig.  7-1,  where  the  probability  of  detection,  Pd,  for  this  hypothesis  test  is 
shown  as  a  function  of  the  false  alarm  probability,  Pp.  For  comparison,  we  also  show 
the  ROC  that  would  result  from  a  detection  of  the  soliton  alone,  at  the  same  noise  level 
and  with  the  time-of-arrival  known.  The  detection  index,  d  =  ^Je/N0,  is  indicated 
for  each  case,  where  E  is  the  energy  in  the  component  soliton.  The  corresponding 
results  for  the  larger  soliton  are  qualitatively  similar,  although  the  detection  indices 
for  that  soliton  alone,  with  =  sinh(2),  are  5.6,  4,  and  3.3,  respectively.  Therefore 
the  detection  probabilities  are  considerably  higher  for  a  fixed  probability  of  false 
alarm.  Note  that  the  detection  performance  for  the  smaller  soliton  is  well  modeled  by 
the  theoretical  performance  for  detection  of  the  smaller  soliton  alone.  This  implies, 
at  least  empirically,  that  the  ability  to  detect  the  component  solitons  in  a  multi- 
soliton  signal  appears  to  be  unaffected  by  the  nonlinear  coupling  with  other  solitons. 
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Further,  although  the  unknown  relative  separation  results  in  significant  waveform 
uncertainty  and  would  require  a  prohibitively  complex  receiver  for  standard  detection 
techniques,  Bayes  optimal  performance  can  still  be  achieved  with  a  minimal  increase 
in  complexity. 


P2=sinh(1.5) 


Figure  7-1:  A  set  of  empirically  generated  ROCs  are  shown  for  the  detection  of  the 
smaller  soliton  from  a  two-soliton  signal.  For  each  of  the  three  noise  levels,  the  ROC 
for  detection  of  the  smaller  soliton  alone  is  also  indicated  along  with  the  corresponding 
detection  index,  d. 

7.3  Summary  and  Further  Considerations 

In  this  chapter  we  address  the  problem  of  multi-soliton  detection  in  an  AWGN  envi¬ 
ronment.  Our  approach  is  to  exploit  the  soliton  evolution  equations  as  a  preproces¬ 
sor,  separating  the  multi-soliton  observation  into  component  solitons.  Focusing  on 
the  Toda  lattice  as  an  example,  at  high  signal-to-noise  ratios,  the  noise  component 
in  the  receiver  remains  low  pass  and  Gaussian  and  is  decoupled  from  the  solitons. 
This  allows  each  of  the  component  solitons  to  be  detected  with  standard  matched 
filter  processing  using  a  GLRT.  Empirically  generated  ROCs  seem  to  indicate  that 
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the  detection  performance  for  the  component  solitons  is  comparable  to  that  for  single 
soliton  observations. 

One  outstanding  issue  is  the  determination  of  appropriate  bounds  on  the  theoret¬ 
ical  performance  of  the  GLRT  for  the  detection  of  soliton  and  multi-soliton  signals. 
Both  when  the  Toda  lattice  is  used  for  signal  separation  of  a  multi-soliton  signal  and 
when  a  single  soliton  is  observed  with  unknown  arrival  time,  the  probability  of  detec¬ 
tion  and  the  false  alarm  rate  have  to  be  computed  empirically,  due  to  the  complexity 
of  the  required  calculations.  The  false  alarm  rate  of  a  GLRT  for  a  single  soliton 
detection  with  unknown  arrival  time  is  determined  by  the  level  crossing  problem, 
whereby  we  seek  the  probability  that  the  output  of  the  matched  filter  will  cross  a 
given  threshold  over  a  fixed  time  interval.  A  description  of  this  problem  is  given  in 
[24].  In  general,  determination  of  the  false  alarm  rate  is  difficult,  however  there  are  a 
few  well-studied  cases.  For  example,  based  on  Siegert’s  work  with  Markov  processes, 
Helstrom  shows  that  for  a  first  order  Gauss-Markov  process  with  auto-covariance, 
Rxx(r)  =  i?(0)e~alTl,  for  large  time  intervals,  aT  »  1,  the  probability  of  false  alarm 
is  given  approximately  by  [24] 

P/«  *  1  -  e~vT ,  V  =  ■ --l-  ,  (7.5) 

y/**R{  0) 

where  it  is  assumed  that  the  threshold,  rj  ^R{ 0).  When  a  matched  filter  is  used  for 
detecting  a  single  soliton  of  unknown  arrival  time,  our  matched  filter  can  be  roughly 
approximated  by  a  filter  with  impulse  response,  h(t)  =  /32e~20^.  The  corresponding 
matched  filter  output  for  a  white  noise  input  will  have  auto-covariance  of  the  form 
/?4(2  +  |r|)e-2^Tl  and  for  (3  »  1  will  behave  much  like  a  first  order  Gauss  Markov 
process.  In  this  manner,  we  might  obtain  a  better  estimate  of  the  expected  false  alarm 
rate  for  single  and  multi-soliton  detection  with  unknown  position.  For  modest  signal 
to  noise  ratios,  the  detection  probability  is  well-modeled  by  the  detection  probability 
of  a  processor  with  the  arrival  time  known. 
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Chapter  8 


Conclusions  and  Directions  for 
Future  Work 


In  this  thesis,  we  have  developed  a  framework  for  exploring  the  generation  and  pro¬ 
cessing  of  soliton  signals  and  introduced  some  potentially  important  applications  of 
solitons  in  the  broad  context  of  communications. 

We  have  taken  the  viewpoint  of  using  solitons  as  carrier  signals  for  transmission 
over  linear,  rather  than  nonlinear  channels.  The  nonlinear  evolution  equations  can 
then  be  viewed  as  specialized  processors  of  this  class  of  signals,  which  are  naturally 
suited  to  performing  a  number  of  complex  signal  processing  tasks.  For  example, 
these  systems  can  efficiently  generate  soliton  signals  and  can  perform  the  nonlinear 
signal  separation  of  multi-soliton  carriers  necessary  for  multiplexing  and  demultiplex¬ 
ing  multiple  users  in  a  potential  soliton  wireless  communications  context. 

Focusing  specifically  on  two  soliton  systems,  the  Toda  lattice  and  the  discrete- 
KdV  equation,  we  develop  new  electrical  analogs  for  the  generation  and  processing  of 
soliton  signals.  Although  analog  circuit  models  have  been  previously  developed  for  a 
variety  of  nonlinear  wave  equations  in  general,  and  for  the  Toda  lattice  in  particular, 
our  diode  ladder  implementation  is  the  first  direct  electrical  analog  of  this  soliton 
system.  Further,  this  appears  to  be  the  first  circuit  model  of  the  Toda  lattice  which  is 
sufficiently  accurate  to  demonstrate  true  soliton  interactions.  The  diode  ladder  circuit 
was  implemented  in  hardware  using  standard  components  and  provides  a  platform  for 
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the  further  development  and  testing  of  real-time  soliton  processing  techniques.  We 
have  also  developed  a  new  circuit  model  for  the  discrete-KdV  equation;  a  nonlinear 
system  which  was  largely  ignored  for  having  no  prior  electrical  or  mechanical  analog. 
The  discrete-KdV  circuit  provides  a  framework  for  processing  discrete-time  soliton 
signals. 

Based  on  these  nonlinear  circuit  models,  we  have  suggested  and  explored  a  number 
of  communications  techniques  whereby  message-bearing  signals  are  used  to  modulate 
the  parameters  of  soliton  carriers.  The  nonlinear  circuits  are  used  to  generate  and 
to  multiplex  multiple  message-signals  onto  a  multi-soliton  carrier.  We  have  demon¬ 
strated  that  the  nonlinear  interaction  of  multiple  solitons  can  be  exploited  as  a  means 
of  reducing  the  transmitted  signal  energy  in  a  multi-user  communication  context. 

To  assess  the  efficacy  and  the  robustness  of  these  communication  techniques  in 
the  presence  of  background  noise,  we  have  analyzed  the  effects  of  small  amplitude 
disturbances  on  the  processing  of  soliton  signals  in  the  Toda  lattice  and  characterized 
the  statistics  of  the  noise  as  it  is  processed.  We  have  shown  that  in  a  high  SNR  AWGN 
background,  the  dynamics  of  soliton  signals  are  practically  unperturbed.  Also,  the 
noise  component  of  the  received  signal  remains  essentially  low  pass  and  Gaussian. 
Through  an  analysis  of  the  nonlinear  spectrum  from  the  inverse  scattering  transform, 
we  also  have  shown  that  at  high  SNR,  virtually  all  of  the  background  noise  maps  onto 
the  non-soliton  (continuum)  component  of  the  nonlinear  spectrum,  and  the  soliton 
spectrum  experiences  only  a  small  Gaussian  perturbation.  This  result  furthers  the 
analogy  of  the  inverse  scattering  transform  to  the  Fourier  transform  for  linear  systems. 

We  were  also  able  to  develop  maximum  likelihood  parameter  estimation  algorithms 
for  soliton  signals  corrupted  by  white  Gaussian  noise.  Calculation  of  the  Cramer-Rao 
bounds  for  the  estimation  error  variance  of  soliton  parameters  reveals  the  potential  for 
estimation  enhancement  due  to  multiple  soliton  interaction.  By  using  the  nonlinear 
circuits  as  preprocessors  which  perform  the  necessary  signal  separation,  the  joint 
estimation  of  multiple  soliton  parameters  can  be  achieved  through  standard  Gaussian 
techniques.  We  have  also  developed  a  maximum  likelihood  parameter  estimation 
algorithm  in  the  nonlinear  spectral  domain  of  the  inverse  scattering  transform.  Both 
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the  theoretical  bounds  and  the  empirical  algorithm  performance  demonstrate  that 
the  nonlinear  coupling  of  multiple  solitons  can  actually  enhance  parameter  estimation 
performance.  This  indicates  the  potential  for  wireless  multi-soliton  communication 
techniques  to  simultaneously  reduce  the  transmitted  signal  energy  while  enhancing 
communication  performance. 

By  a  similar  technique,  soliton  circuits  can  be  used  to  decompose  a  multi-soliton 
detection  problem  with  (M  +  1)^  hypotheses  into  ^-separate  binary  hypothesis 
tests  without  sacrificing  Bayes  optimality  in  high  SNR  Gaussian  environments.  Of 
particular  interest  is  the  case  of  multi-soliton  position  uncertainty  which  can  gives  rise 
to  dramatic  changes  in  the  received  multi-soliton  signal  and  therefore  would  have  a 
significant  impact  on  the  structure  of  a  traditional  detector.  However,  since  the  soliton 
circuits  decompose  the  multi-soliton  detection  problem  into  single  soliton  detections, 
our  detection  framework  is  unchanged. 

8.1  Future  Directions 

In  this  thesis,  we  began  an  exploration  of  a  number  of  new  and  interesting  applications 
of  soliton  signals  and  systems  from  a  signal  processing  viewpoint.  As  a  result,  there 
are  a  number  outstanding  issues  relating  to  many  of  the  topics  considered.  Although 
many  areas  of  future  study  have  been  identified  in  the  body  of  the  thesis,  we  mention 
several  particularly  intriguing  directions  here. 

Much  of  the  mathematical  foundation  for  soliton  dynamics  lies  within  the  frame¬ 
work  of  solvable  nonlinear  systems.  In  Sec.  2.2  we  have  attempted  to  present  some 
of  the  essential  ideas  behind  this  class  of  nonlinear  systems  and  relate  them  to  linear 
system  theory.  Although  many  of  the  results  in  this  thesis  pertain  to  the  detection 
and  estimation  of  parameters  of  received  soliton  signals,  there  are  many  interesting 
signal  processing  questions  regarding  the  dynamics  these  integrable  nonlinear  systems 
themselves.  Specifically,  it  may  be  possible  to  address  the  problems  of  recursive  esti¬ 
mation,  e.g.,  prediction,  extrapolation,  filtering,  and  smoothing  of  such  systems  with 
the  aid  of  the  inverse  scattering.  Traditional  barriers  to  such  investigations  for  nonlin- 
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ear  systems  lie  in  both  the  state  propagation  and  the  propagation  of  state  probability 
densities.  For  integrable  Hamiltonian  systems,  the  state  propagation  can  be  solved 
explicitly  via  inverse  scattering,  while  the  propagation  of  the  state  density  remains  a 
difficult  problem. 

Another  potentially  interesting  area  of  future  research  involves  extending  the  re¬ 
sults  of  integrable  system  theory  to  include  totally  discrete  (discrete-time,  finite  di¬ 
mensional)  systems.  Some  work  in  this  area  has  led  to  discretized  versions  of  the 
inverse  scattering  theory  for  the  Toda  lattice  [20,  26,  28]  the  Korteweg-de  Vries 
equation[25],  as  well  as  others  in  [1,  27].  Such  systems  may  lead  to  applications 
using  discrete  soliton  signals  for  digital  communication. 

In  Chap.  3  we  illustrate  how  analog  circuits  could  be  used  in  the  context  of  the 
communication  ideas  shown  in  Chap.  4.  These  results  along  with  those  of  Hirota  and 
Suzuki  [29,  64,  65],  Toda  [68],  Scott  [56,  57],  and  others  illustrate  the  potential  for 
development  of  a  variety  of  analog  hardware  for  soliton  modulation  techniques.  How¬ 
ever,  there  are  several  outstanding  issues  regarding  the  detailed  implementation  of 
the  communications  techniques  discussed  in  Chap.  4  using  such  circuits.  Specifically, 
the  demodulation  of  individual  soliton  carriers  still  must  be  addressed  along  with  a 
more  complete  investigation  into  the  effects  of  realistic  non-additive  channel  distur¬ 
bances.  There  may  also  be  benefit  to  further  exploration  of  the  soliton-like  behavior 
of  certain  cellular  automata  [1,  7,  53,  66]  and  the  possible  use  of  such  systems  for 
digital  communication  or  computation  [32,  62,  63]. 

The  nonlinear  circuit  models  developed  for  the  Toda  lattice  demonstrate  an  in¬ 
teresting  similarity  to  the  models  used  by  Keshner  in  his  analysis  of  1//  noise  [38]. 
Specifically,  a  small  signal  model  for  the  diode  ladder  circuit  yields  a  driving  point 
impedance  inversely  proportional  to  frequency,  and  thus  would  produce  a  “double 
Brownian”  motion  process  in  response  to  a  white  noise  input.  Such  an  analogy  leads 
directly  to  the  observation  that  at  least  for  the  small  signal  model,  noise  at  the  in¬ 
put  of  the  nonlinear  circuits  will  lead  to  1/ /-type  disturbances.  Another  interesting 
link  between  the  Iff  generators  and  soliton  circuits,  is  that,  like  the  family  of  1/f 
random  processes,  the  family  of  soliton  solutions  to  such  equations  also  have  a  multi- 
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scale  component.  This  apparent  multi-scale  invariance  of  the  soliton  dynamics  along 
with  the  similarity  in  the  circuit  models  to  those  which  generate  1//  noise,  may  in¬ 
dicate  the  presence  of  a  fundamental  nonlinear  multi-scale  behavior  inherent  in  the 
dynamics. 

Although  we  have  explored  some  basic  strategies  for  using  soliton  signals  for  com¬ 
munication,  there  remain  a  number  of  unanswered  questions  regarding  the  efficacy 
of  such  techniques  in  comparison  to  standard  methods.  For  example,  a  quantitative 
investigation  of  the  spectral  efficiency,  the  effects  of  inter-symbol  interference,  the 
error  probability,  or  decoder  complexity  for  such  methods  in  comparison  to  standard 
techniques  ought  to  be  undertaken.  Finally,  we  note  that  solitons  represent  a  rich 
class  of  signals  that  are  potentially  applicable  to  a  variety  of  signal  processing  and 
communications  contexts  beyond  those  mentioned  in  this  thesis.  Indeed,  perhaps 
the  most  exciting  and  potentially  rewarding  direction  for  future  research  lies  in  the 
discovery  of  new  and  promising  applications. 
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